Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,487,342 C→T 100% intergenic (+234/‑28) tsgA → / → nirB putative transporter/nitrite reductase, large subunit, NAD(P)H‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,487,3420CT96.4% 721.2 / 12.7 225intergenic (+234/‑28)tsgA/nirBputative transporter/nitrite reductase, large subunit, NAD(P)H‑binding
Reads supporting (aligned to +/- strand):  ref base C (4/4);  new base T (111/105);  total (116/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 231 total aligned reads displayed.

TCTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGA  >  CP009273/3487198‑3487485
                                                                                                                                                |                                                                                                                                               
tCTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAggtggt                                                                                                                                              >  2:3160072/1‑148 (MQ=255)
 cTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGACTCGTTAAGGTAGGTGGTa                                                                                                                                             >  1:3461290/1‑148 (MQ=255)
 cTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTa                                                                                                                                             >  2:2820002/1‑148 (MQ=255)
  tATTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTaa                                                                                                                                            >  1:2543362/1‑148 (MQ=255)
   aTTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAAt                                                                                                                                           >  1:316412/1‑148 (MQ=255)
   aTTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAAt                                                                                                                                           >  2:461387/1‑148 (MQ=255)
      ttttGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATaga                                                                                                                                        <  2:3385851/148‑1 (MQ=255)
        ttGCACAAACATGAAATATCCGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATagaaa                                                                                                                                      <  2:1429722/148‑1 (MQ=255)
        ttGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATagaaa                                                                                                                                      <  2:1965177/148‑1 (MQ=255)
        gtGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATagaaa                                                                                                                                      <  1:2459029/147‑1 (MQ=255)
             caAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATagaaaagaaa                                                                                                                                 <  1:587830/148‑1 (MQ=255)
              aaaCATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAAt                                                                                                                                >  2:2325919/1‑148 (MQ=255)
                aCATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCg                                                                                                                              >  1:1256629/1‑148 (MQ=255)
                 cATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGa                                                                                                                             >  1:1217575/1‑148 (MQ=255)
                  aTGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAg                                                                                                                            <  1:2362546/148‑1 (MQ=255)
                   tGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAgg                                                                                                                           <  2:3318410/148‑1 (MQ=255)
                    gAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGc                                                                                                                          >  1:2864225/1‑148 (MQ=255)
                           cAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATg                                                                                                                   >  2:1262915/1‑148 (MQ=255)
                           cAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATg                                                                                                                   <  2:682106/148‑1 (MQ=255)
                             gACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAg                                                                                                                 >  1:1579291/1‑148 (MQ=255)
                                aaTTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCaa                                                                                                              <  1:245275/148‑1 (MQ=255)
                                 aTTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCaaa                                                                                                             >  1:2949427/1‑148 (MQ=255)
                                 aTTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCaaa                                                                                                             >  1:1464457/1‑148 (MQ=255)
                                  ttCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAGGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAg                                                                                                            >  2:1301623/1‑148 (MQ=255)
                                     cGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCa                                                                                                         >  2:134401/1‑148 (MQ=255)
                                       tGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGa                                                                                                       >  1:585377/1‑148 (MQ=255)
                                                aaaaTTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAAtt                                                                                              <  1:1526684/148‑1 (MQ=255)
                                                   aTTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATc                                                                                           <  1:1025293/148‑1 (MQ=255)
                                                        taCAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTaa                                                                                      >  2:576732/1‑148 (MQ=255)
                                                           aaaTCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACgg                                                                                   >  1:152879/1‑148 (MQ=255)
                                                             aTCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTa                                                                                 >  2:964692/1‑148 (MQ=255)
                                                             aTCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGGAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTa                                                                                 >  2:832644/1‑148 (MQ=255)
                                                               cAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATg                                                                               >  2:2147835/1‑148 (MQ=255)
                                                                 gCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGt                                                                             >  1:1089792/1‑148 (MQ=255)
                                                                 gCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGt                                                                             >  1:3049578/1‑148 (MQ=255)
                                                                 gCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGCTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGt                                                                             <  1:2103203/148‑1 (MQ=255)
                                                                 gCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGCTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGt                                                                             <  1:2103046/148‑1 (MQ=255)
                                                                  cAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTc                                                                            >  2:3270574/1‑148 (MQ=255)
                                                                     tgtaCCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGc                                                                         <  2:3088214/146‑1 (MQ=255)
                                                                     tataCCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGc                                                                         >  1:1390470/1‑148 (MQ=255)
                                                                      ataCCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGcc                                                                        <  1:2202365/148‑1 (MQ=255)
                                                                       taCCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCa                                                                       <  2:2086017/148‑1 (MQ=255)
                                                                        aCCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCAt                                                                      >  1:238758/1‑148 (MQ=255)
                                                                        aCCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCAt                                                                      >  2:3005451/1‑148 (MQ=255)
                                                                        aCCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCAt                                                                      >  1:2665202/1‑148 (MQ=255)
                                                                          ccATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCg                                                                    <  2:3293541/148‑1 (MQ=255)
                                                                              tAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCttt                                                                >  1:2663465/1‑148 (MQ=255)
                                                                              tAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCttt                                                                >  2:781543/1‑148 (MQ=255)
                                                                                 ggAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATc                                                             >  1:1177134/1‑148 (MQ=255)
                                                                                 ggAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATc                                                             >  1:2954385/1‑148 (MQ=255)
                                                                                 gcagaatataAAGGTGAATTTGATTTACATCAATAAGAGGGGTTGATGAATCTTTAAGCTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAACAAACTCAGAAACGCAATAATCGGTAACGGTATGGTCGGCCATCGCTTTATc                                                             >  2:1005365/6‑148 (MQ=255)
                                                                                   aGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGa                                                           >  1:2701718/1‑148 (MQ=255)
                                                                                      atataAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGa                                                        <  1:2751458/148‑1 (MQ=255)
                                                                                       tataAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGAt                                                       <  2:2654749/148‑1 (MQ=255)
                                                                                        ataAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGAtc                                                      <  2:2890029/148‑1 (MQ=255)
                                                                                          aaaGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGAtctt                                                    <  1:1770782/148‑1 (MQ=255)
                                                                                            aGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGAtcttct                                                  >  1:352025/1‑148 (MQ=255)
                                                                                            aGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGAtcttct                                                  <  2:2649276/148‑1 (MQ=255)
                                                                                            aGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGAtcttct                                                  >  2:808147/1‑148 (MQ=255)
                                                                                             ggTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATcttctt                                                 >  2:498646/1‑148 (MQ=255)
                                                                                             ggTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATcttctt                                                 >  2:146940/1‑148 (MQ=255)
                                                                                             ggTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATcttctt                                                 >  2:1364594/1‑148 (MQ=255)
                                                                                             ggTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATctcctt                                                 >  2:332399/1‑148 (MQ=255)
                                                                                              gTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTg                                                <  2:3466348/148‑1 (MQ=255)
                                                                                               tGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGa                                               <  1:2404225/148‑1 (MQ=255)
                                                                                                  aTTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATaa                                            <  2:2228878/148‑1 (MQ=255)
                                                                                                  aTTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATaa                                            <  2:3445268/148‑1 (MQ=255)
                                                                                                  aTTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATaa                                            >  2:1311688/1‑148 (MQ=255)
                                                                                                   tttGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATaaa                                           <  2:1142861/148‑1 (MQ=255)
                                                                                                   tttGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATaaa                                           <  2:2921445/148‑1 (MQ=255)
                                                                                                   tttGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATaaa                                           <  2:2505756/148‑1 (MQ=255)
                                                                                                    ttGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAAt                                          >  2:35261/1‑148 (MQ=255)
                                                                                                     tGATTTACATCAATCAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATc                                         <  2:1487510/148‑1 (MQ=255)
                                                                                                       aTTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGTCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTg                                       <  2:1633178/148‑1 (MQ=255)
                                                                                                       aTTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTg                                       <  2:2543362/148‑1 (MQ=255)
                                                                                                       aTTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTg                                       <  2:1549994/148‑1 (MQ=255)
                                                                                                        tttACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGa                                      <  1:2147835/148‑1 (MQ=255)
                                                                                                          tACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATAGAGGCAAATATGAGCAAAGTCAGACTCGCAATTATCGTTATCGGTATGGTCGGCCATCTCTTTTTCGAAGATCTTCTTGATAAATCTGATg                                    <  2:941747/148‑1 (MQ=255)
                                                                                                               cAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGcc                               >  1:3279494/1‑148 (MQ=255)
                                                                                                                  tAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAAc                            <  2:775293/148‑1 (MQ=255)
                                                                                                                  tAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAAc                            >  1:2914987/1‑148 (MQ=255)
                                                                                                                  tAAGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAAc                            >  1:3258430/1‑148 (MQ=255)
                                                                                                                   aaGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACt                           >  2:599276/1‑148 (MQ=255)
                                                                                                                   aaGCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACt                           >  2:2720355/1‑148 (MQ=255)
                                                                                                                     gCGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACttt                         >  2:1031263/1‑148 (MQ=255)
                                                                                                                      cGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTg                        <  1:2484555/148‑1 (MQ=255)
                                                                                                                      cGGGGTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTg                        >  2:1842038/1‑148 (MQ=255)
                                                                                                                        gggTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGat                      >  1:937822/1‑148 (MQ=255)
                                                                                                                         ggTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGata                     >  2:2855403/1‑148 (MQ=255)
                                                                                                                          gTTGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGatat                    <  2:2412452/148‑1 (MQ=255)
                                                                                                                           ttGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATAtt                   >  2:1926917/1‑148 (MQ=255)
                                                                                                                           ttGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATAtt                   >  2:1062299/1‑148 (MQ=255)
                                                                                                                            tGCTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTa                  <  2:2429074/148‑1 (MQ=255)
                                                                                                                              cTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTAcc                <  1:48331/148‑1 (MQ=255)
                                                                                                                              cTGAATCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTAcc                <  1:1092836/148‑1 (MQ=255)
                                                                                                                                  aTCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGttt            >  2:1791120/1‑148 (MQ=255)
                                                                                                                                  aTCGTTAAGGTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGttt            >  2:1894112/1‑148 (MQ=255)
                                                                                                                                         aggtaggtGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTga     <  1:668930/148‑1 (MQ=255)
                                                                                                                                           gTAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTgaag   >  1:1905377/1‑148 (MQ=255)
                                                                                                                                            tAGGTGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTgaaga  <  2:2582650/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                               
TCTATTTTTTGCACAAACATGAAATATCAGACAATTCCGTGACTTAAGAAAATTTATACAAATCAGCAATATACCCATTAAGGAGTATATAAAGGTGAATTTGATTTACATCAATAAGCGGGGTTGCTGAATCGTTAAGGTAGGCGGTAATAGAAAAGAAATCGAGGCAAAAATGAGCAAAGTCAGACTCGCAATTATCGGTAACGGTATGGTCGGCCATCGCTTTATCGAAGATCTTCTTGATAAATCTGATGCGGCCAACTTTGATATTACCGTTTTCTGTGAAGA  >  CP009273/3487198‑3487485

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: