Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,531,029 G→T 100% R180S (CGC→AGC)  narW ← nitrate reductase 2 (NRZ), delta subunit (assembly subunit)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,531,0290GT96.4% 668.8 / 13.3 233R180S (CGC→AGC) narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base G (5/4);  new base T (105/119);  total (110/123)
Fisher's exact test for biased strand distribution p-value = 7.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 233 total aligned reads displayed.

CCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGAT  >  CP009273/1530885‑1531173
                                                                                                                                                |                                                                                                                                                
ccGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGgt                                                                                                                                               <  2:165358/148‑1 (MQ=255)
  gCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGgtgt                                                                                                                                             >  2:2010199/1‑148 (MQ=255)
    aCTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCa                                                                                                                                           >  1:1093939/1‑148 (MQ=255)
     cTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATAATTATTAATCCGCGACCTGTCACACCCCGTGGCCTTAACATAAATAAACTTCACCAGTTACAATTCCCACCCAGCCGCAACCGACTGGCTGGTGGCAt                                                                                                                                          >  1:2103766/1‑148 (MQ=255)
      tGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATc                                                                                                                                         >  1:3202758/1‑148 (MQ=255)
       gATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATcg                                                                                                                                        >  2:316871/1‑148 (MQ=255)
        aTGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATcgc                                                                                                                                       >  2:760962/1‑148 (MQ=255)
          gTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCgcgc                                                                                                                                     >  2:2086392/1‑148 (MQ=255)
            cGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGctc                                                                                                                                   <  2:569338/148‑1 (MQ=255)
            cGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGctc                                                                                                                                   <  2:3504871/148‑1 (MQ=255)
            cGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGctc                                                                                                                                   >  1:1928664/1‑148 (MQ=255)
            cGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGctc                                                                                                                                   <  2:853158/148‑1 (MQ=255)
             gACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGctct                                                                                                                                  <  1:2811610/148‑1 (MQ=255)
              aCATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCtt                                                                                                                                 >  1:1738208/1‑148 (MQ=255)
                  atTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGAATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGcg                                                                                                                             >  2:1918804/1‑147 (MQ=255)
                   ttGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTg                                                                                                                            <  1:799099/148‑1 (MQ=255)
                   ttGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTg                                                                                                                            >  1:2245926/1‑148 (MQ=255)
                    tGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGt                                                                                                                           <  1:3469451/148‑1 (MQ=255)
                       ggcgCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTa                                                                                                                        <  1:117997/148‑1 (MQ=255)
                       ggcgCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTa                                                                                                                        >  1:2579620/1‑148 (MQ=255)
                         cgctACATCCTGGCTAAAGCGTCGCTGATATTGATTAATCGGCGAGCTGTCACACGCCGTGTCGTTATCTTCAATAAACTTATACTGTTCCTCTTCCAACAATGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACt                                                                                                                      <  2:222975/148‑1 (MQ=255)
                          gcgACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACtt                                                                                                                     <  1:557016/148‑1 (MQ=255)
                           cgACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTg                                                                                                                    <  2:2360808/148‑1 (MQ=255)
                            gACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGt                                                                                                                   <  1:2483925/148‑1 (MQ=255)
                                tCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtg                                                                                                               <  2:11301/148‑1 (MQ=255)
                                   tGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGa                                                                                                            <  1:1852601/148‑1 (MQ=255)
                                   tGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGa                                                                                                            <  1:1379614/148‑1 (MQ=255)
                                      cTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGacac                                                                                                         >  2:1045172/1‑148 (MQ=255)
                                        aaaGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTg                                                                                                       <  2:2285052/148‑1 (MQ=255)
                                         aaGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGt                                                                                                      >  2:920468/1‑148 (MQ=255)
                                         aaGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGt                                                                                                      >  1:2720449/1‑148 (MQ=255)
                                           gCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCa                                                                                                    >  2:2756739/1‑148 (MQ=255)
                                              tCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACtt                                                                                                 >  2:1462568/1‑148 (MQ=255)
                                              tCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACtt                                                                                                 >  2:3229455/1‑148 (MQ=255)
                                                 cTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGaa                                                                                              <  1:1986476/148‑1 (MQ=255)
                                                 cTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGaa                                                                                              <  1:3392024/148‑1 (MQ=255)
                                                    atcttGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAgg                                                                                           <  2:1139180/145‑1 (MQ=255)
                                                    atatTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAgg                                                                                           <  1:3149298/148‑1 (MQ=255)
                                                     tatTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAggg                                                                                          <  1:2679969/148‑1 (MQ=255)
                                                     tatTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAggg                                                                                          <  1:461254/148‑1 (MQ=255)
                                                      atTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGa                                                                                         >  1:3318344/1‑148 (MQ=255)
                                                        tGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGGTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGc                                                                                       <  2:1899490/148‑1 (MQ=255)
                                                        tGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGCTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGc                                                                                       <  1:1820366/148‑1 (MQ=255)
                                                         gATTAAGCGGCTAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTGCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCGGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCt                                                                                      <  2:2103766/148‑1 (MQ=255)
                                                          aTTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCtt                                                                                     >  1:2595181/1‑148 (MQ=255)
                                                            tAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                   <  2:180279/148‑1 (MQ=255)
                                                             aaGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTccc                                                                                  <  2:3441212/148‑1 (MQ=255)
                                                              aGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCg                                                                                 >  2:2779362/1‑148 (MQ=255)
                                                                cggcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGcc                                                                               >  1:2298857/1‑148 (MQ=255)
                                                                   cgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGt                                                                            >  2:803991/1‑148 (MQ=255)
                                                                            acacGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCaga                                                                   <  1:3020002/148‑1 (MQ=255)
                                                                               cGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCa                                                                >  1:1354695/1‑148 (MQ=255)
                                                                                    tGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCaaa                                                           >  2:133672/1‑148 (MQ=255)
                                                                                     ggCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAaac                                                          >  1:2661367/1‑148 (MQ=255)
                                                                                         ttATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAcg                                                      >  2:2652719/1‑148 (MQ=255)
                                                                                          tATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAcgc                                                     >  1:1684884/1‑148 (MQ=255)
                                                                                          tATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAcgc                                                     >  1:517809/1‑148 (MQ=255)
                                                                                          tATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGATGCAGCAGAGCATCAAACAAcgc                                                     >  1:1055165/1‑148 (MQ=255)
                                                                                           aTCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAcgcg                                                    <  2:1578147/148‑1 (MQ=255)
                                                                                             cTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTa                                                  <  1:3082731/148‑1 (MQ=255)
                                                                                              ttCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTAc                                                 >  2:2901870/1‑148 (MQ=255)
                                                                                              ttCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTAc                                                 >  2:2902142/1‑148 (MQ=255)
                                                                                               tCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTAcc                                                <  2:1928664/148‑1 (MQ=255)
                                                                                                cAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCa                                               >  1:3457195/1‑148 (MQ=255)
                                                                                                  aTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACg                                             <  1:127136/148‑1 (MQ=255)
                                                                                                   tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                            <  1:2035512/148‑1 (MQ=255)
                                                                                                   tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                            >  2:261209/1‑148 (MQ=255)
                                                                                                   tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                            <  1:1530914/148‑1 (MQ=255)
                                                                                                   tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                            <  1:2036191/148‑1 (MQ=255)
                                                                                                     aaCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTg                                          <  2:668283/148‑1 (MQ=255)
                                                                                                      aCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGc                                         <  2:3498240/148‑1 (MQ=255)
                                                                                                       cTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGcc                                        <  1:2614920/148‑1 (MQ=255)
                                                                                                           aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgc                                    >  2:1544007/1‑148 (MQ=255)
                                                                                                            ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgcg                                   <  1:1434898/148‑1 (MQ=255)
                                                                                                            ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTAAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgcg                                   <  1:952976/148‑1 (MQ=255)
                                                                                                             cTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                  >  2:3536735/1‑148 (MQ=255)
                                                                                                             cTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                  <  1:1330968/148‑1 (MQ=255)
                                                                                                                 tCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTCCCACGGTGCCTCGCgttgtt                              >  2:2852721/1‑148 (MQ=255)
                                                                                                                 tCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                              >  1:610009/1‑148 (MQ=255)
                                                                                                                 tCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                              >  2:728131/1‑148 (MQ=255)
                                                                                                                 tCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                              >  1:609286/1‑148 (MQ=255)
                                                                                                                      ttCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCGTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAg                         <  2:1185460/148‑1 (MQ=255)
                                                                                                                        cccACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCg                       >  1:2873335/1‑148 (MQ=255)
                                                                                                                         ccACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                      >  1:2989071/1‑148 (MQ=255)
                                                                                                                          cacaCTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAc                     <  1:1466519/148‑1 (MQ=255)
                                                                                                                          cacaCTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAc                     <  2:710660/148‑1 (MQ=255)
                                                                                                                            cacTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCg                   >  1:867689/1‑148 (MQ=255)
                                                                                                                              cTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGcc                 >  1:685420/1‑148 (MQ=255)
                                                                                                                                   tCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTTTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCa            <  1:2002464/148‑1 (MQ=255)
                                                                                                                                      cgcgcggggCTGGTGTCATCGCGCTCTGCGCTGTTTATTTGTTTTGTGCCACGGTCACTTGAACGGGCGGTTCCCGCCAGTTGCACCAGAGCATCAAACACCGCGTACCACGGTGCCTCGCGTTTTTTTACGCGACCGCCAAGCAggg         <  1:2685700/141‑1 (MQ=255)
                                                                                                                                      cgcgcCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGAAACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAggg         <  1:1895918/147‑1 (MQ=255)
                                                                                                                                      cgcgcCGGGCGGGTGTCATCGCGCTCTTCGCTTTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGGGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAggg         <  1:1894836/147‑1 (MQ=255)
                                                                                                                                       gcgcCGGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGc        <  2:2444462/148‑1 (MQ=255)
                                                                                                                                          ccTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTTAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCCCGGTGCCTCGCTTTTTTTTAAGCGACCGCCAAGCAGGGCGAg     <  1:534366/148‑1 (MQ=255)
                                                                                                                                          ccTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAg     >  2:798978/1‑148 (MQ=255)
                                                                                                                                            tGGCTGGTGTCATCGCGCTCTTCGCTTTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa   <  2:1622673/148‑1 (MQ=255)
                                                                                                                                            tGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa   <  2:1323916/148‑1 (MQ=255)
                                                                                                                                            tGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGGACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa   <  2:2142863/148‑1 (MQ=255)
                                                                                                                                            gggCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTGGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACACCGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa   <  1:760962/147‑1 (MQ=255)
                                                                                                                                             ggCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGAt  >  2:152956/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGAT  >  CP009273/1530885‑1531173

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