Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,684,575 G→A 5.2% S311S (TCC→TCT)  glyA ← serine hydroxymethyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,684,5750GA 5.2% 253.4 / 2.3 77S311S (TCC→TCT) glyAserine hydroxymethyltransferase
Reads supporting (aligned to +/- strand):  ref base G (34/39);  new base A (2/2);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

GTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACC  >  NC_000913/2684432‑2684721
                                                                                                                                               |                                                                                                                                                  
gTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAAc                                                                                                                                               >  8:52352/1‑149 (MQ=255)
gTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAAc                                                                                                                                               >  1:297663/1‑149 (MQ=255)
gTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAAc                                                                                                                                               >  8:45946/1‑149 (MQ=255)
     aaaCGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACtt                                                                                                                                          <  6:247784/149‑1 (MQ=255)
     aaaCGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACtt                                                                                                                                          <  7:399963/149‑1 (MQ=255)
      aaCGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACttt                                                                                                                                         <  6:359614/149‑1 (MQ=255)
        cGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGt                                                                                                                                       <  7:128934/149‑1 (MQ=255)
         gggCTCTTCGGATCGTTCGGTACGCTGTTGTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTa                                                                                                                                      <  7:358328/149‑1 (MQ=255)
                     tCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACATGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAg                                                                                                                          <  8:41952/149‑1 (MQ=255)
                     tCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAg                                                                                                                          <  1:197891/149‑1 (MQ=255)
                     tCGTTAGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGTCGCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCATCCAGGAACAGGTGGTTATCAGTGCCGACGGAAACCTCTTTGTAGCGGCGCTCgt                                                                                                                           <  8:286273/148‑2 (MQ=255)
                      cGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAg                                                                                                                          <  2:297663/148‑1 (MQ=255)
                        ttCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGATTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGaa                                                                                                                       >  6:154363/1‑149 (MQ=255)
                                  tGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTAACGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGTTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTAc                                                                                                              >  6:270282/1‑148 (MQ=255)
                                    tttttGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTAcc                                                                                                             <  2:227221/147‑1 (MQ=255)
                                      tttGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCg                                                                                                         >  7:110031/1‑149 (MQ=255)
                                            gACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTa                                                                                                    <  8:95518/148‑1 (MQ=255)
                                                    tGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGtttt                                                                                             >  8:96453/1‑147 (MQ=255)
                                                    tGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTa                                                                                           <  7:304039/149‑1 (MQ=255)
                                                    ggTTAGCACGGCCCAGAGCGGCGTCTGCGTCTTTACCGGGCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTa                                                                                           <  4:69117/148‑1 (MQ=255)
                                                     gTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAg                                                                                          >  1:204076/1‑149 (MQ=255)
                                                         gCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGa                                                                                       <  8:236622/148‑1 (MQ=255)
                                                         gCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGAc                                                                                      <  7:52352/149‑1 (MQ=255)
                                                         gCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGAc                                                                                      >  8:301250/1‑149 (MQ=255)
                                                         gCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGAc                                                                                      <  7:45946/149‑1 (MQ=255)
                                                         gCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGGCAGGTTTTTATCAACAAGATCAAACAGGATCAGGTGGTTATCATGGCAGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGa                                                                                       <  5:129064/148‑1 (MQ=255)
                                                           aCGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACct                                                                                    <  7:96453/149‑1 (MQ=255)
                                                             ggCCCAGAGCGGCGTCTGCATCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTATCCGCGCTCGATGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACctgc                                                                                  <  6:124173/149‑1 (MQ=255)
                                                                          gTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGttt                                                                     >  8:222832/1‑149 (MQ=255)
                                                                               cttcttTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACttt                                                                                                                                         <  3:79935/76‑1 (MQ=255)
                                                                               cttcttTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACttt                                                                                                                                         >  4:79935/1‑76 (MQ=255)
                                                                               cttcttTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTg                                                                   <  7:301250/146‑1 (MQ=255)
                                                                                 tcttTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGAGACCTGCTGCTGGTAAGTTTTGAACt                                                               >  2:137675/1‑148 (MQ=255)
                                                                                  cttTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTc                                                              <  7:272593/148‑1 (MQ=255)
                                                                                  cttTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTc                                                              <  1:137675/148‑1 (MQ=255)
                                                                                    ttACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGAcc                                                                                     <  1:131986/123‑1 (MQ=255)
                                                                                    ttACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGAcc                                                                                     >  2:131986/1‑123 (MQ=255)
                                                                                        cGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTc                                                        >  5:158342/1‑148 (MQ=255)
                                                                                          gTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCAt                                                     >  1:78610/1‑149 (MQ=255)
                                                                                          gTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCAt                                                     >  5:351916/1‑149 (MQ=255)
                                                                                           tCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATc                                                    >  7:140906/1‑149 (MQ=255)
                                                                                             aGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGc                                                  >  6:27097/1‑149 (MQ=255)
                                                                                               gTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGttt                                                                     <  3:112894/128‑1 (MQ=255)
                                                                                               gTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGttt                                                                     >  4:112894/1‑128 (MQ=255)
                                                                                               gTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGc                                                          >  4:90069/1‑139 (MQ=255)
                                                                                               gTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGc                                                          <  3:90069/139‑1 (MQ=255)
                                                                                                   ttATCAACCAGATCAACCAGGACCAGGTGGGTTTCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTCCCCTCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCttt                                            >  2:214962/1‑149 (MQ=255)
                                                                                                    tATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCttt                                            <  8:110031/148‑1 (MQ=255)
                                                                                                     atcaaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCttt                                            >  8:17777/1‑147 (MQ=255)
                                                                                                     atcaaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTACTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTc                                           <  8:140906/148‑1 (MQ=255)
                                                                                                      tcaaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCag                                         <  2:204076/149‑1 (MQ=255)
                                                                                                      tcaaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCa                                          <  3:136903/148‑1 (MQ=255)
                                                                                                        aaccagatcaaccagTAACAGGTGGTTATCAGTGCCGCAGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCaga                                        <  4:349797/148‑1 (MQ=255)
                                                                                                        aaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCa                                                      <  4:100828/134‑1 (MQ=255)
                                                                                                        aaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCa                                                      >  3:100828/1‑134 (MQ=255)
                                                                                                        aaccagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCagag                                       >  6:189161/1‑149 (MQ=255)
                                                                                                         accagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCagag                                       >  4:165075/1‑148 (MQ=255)
                                                                                                         accagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCagag                                       >  4:302994/1‑148 (MQ=255)
                                                                                                         accagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCagag                                       >  8:324487/1‑148 (MQ=255)
                                                                                                           cagatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCaa                                    <  2:78610/149‑1 (MQ=255)
                                                                                                            agatcaaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAAc                                   <  3:302994/149‑1 (MQ=255)
                                                                                                               tcaaccagGAACAGGTGGTTATCAGTGCCGCCAGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCagagagatag                                 >  1:239116/1‑142 (MQ=255)
                                                                                                               tcaaccagGAACAGGTGGTTATCAGTGCCGCCAGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCaga                                        <  2:245176/141‑1 (MQ=255)
                                                                                                               tcaaccagGAACAGGTGGTTATCAGTGCCGCCAGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCaga                                        <  2:239116/141‑1 (MQ=255)
                                                                                                               tcaaccagGAACAGGTGGTTATCAGTGCCGCCAGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCaga                                        >  1:245176/1‑141 (MQ=255)
                                                                                                                caaccagGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCt                               >  5:400038/1‑149 (MQ=255)
                                                                                                                        aaCAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGc                       <  5:189161/149‑1 (MQ=255)
                                                                                                                                gtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGt               <  7:222832/149‑1 (MQ=255)
                                                                                                                                gtTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGt               <  6:400038/149‑1 (MQ=255)
                                                                                                                                 ttATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTg              >  5:228498/1‑149 (MQ=255)
                                                                                                                                 ttATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTg              >  3:64715/1‑149 (MQ=255)
                                                                                                                                   aTCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGc             >  6:401043/1‑148 (MQ=255)
                                                                                                                                    tCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCAt           >  3:33241/1‑149 (MQ=255)
                                                                                                                                      aGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCttt                                            >  3:258970/1‑114 (MQ=255)
                                                                                                                                      aGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCttt                                            <  4:258970/114‑1 (MQ=255)
                                                                                                                                        tGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAAc       >  2:16105/1‑149 (MQ=255)
                                                                                                                                         gccgccGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAAc       <  5:270282/148‑1 (MQ=255)
                                                                                                                                           cgccgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCATCTTTCAGAGCAACCGATTAACCGGCGATTACGTGCATCAACGGa    >  8:90409/1‑149 (MQ=255)
                                                                                                                                             ccgGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGAcc  >  1:56367/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
GTCACAAACGGGCTCTTCGGATCGTTCGGTACGCTGTTTTTGTTGACGGTGATGTTAGCACGGCCCAGAGCGGCGTCTGCTTCTTTACCGGTCAGGTTTTTATCAACCAGATCAACCAGGAACAGGTGGTTATCAGTGCCGCCGGAAACCACTTTGTAGCCGCGCTCGAGGAACACTTCTACCATCGCTTTAGCGTTTTTAGCGACCTGCTGCTGGTAAGTTTTGAACTCAGGCTCCATCGCTTCTTTCAGAGCAACCGCTTTACCGGCGATTACGTGCATCAACGGACC  >  NC_000913/2684432‑2684721

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑