Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,531,029 G→T 84.2% 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,0290GT84.2% 432.9 / 95.0 218R180S (CGC→AGC) narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base G (14/21);  new base T (76/107);  total (90/128)
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.06e-01

Only 100 of 221 total aligned reads displayed.

CCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGG  >  CP009273/1530884‑1531176
                                                                                                                                                 |                                                                                                                                                   
cccGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCctggctgg                                                                                                                                                   <  1:2603932/148‑1 (MQ=255)
  cGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGgtg                                                                                                                                                 <  2:195392/148‑1 (MQ=255)
    cACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTc                                                                                                                                               >  2:2995194/1‑148 (MQ=255)
     aCTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCa                                                                                                                                              <  1:2155353/148‑1 (MQ=255)
       tGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATc                                                                                                                                            >  1:314543/1‑148 (MQ=255)
         aTGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATcgc                                                                                                                                          >  1:1024616/1‑148 (MQ=255)
         aTGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATcgc                                                                                                                                          >  2:3038439/1‑148 (MQ=255)
          tGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATcgcg                                                                                                                                         <  1:2527301/148‑1 (MQ=255)
             cGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGctc                                                                                                                                      <  1:1359372/148‑1 (MQ=255)
             cGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGctc                                                                                                                                      >  1:1801166/1‑148 (MQ=255)
               aCATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCtt                                                                                                                                    <  2:1340576/148‑1 (MQ=255)
               aCATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCtt                                                                                                                                    <  1:612875/148‑1 (MQ=255)
                  tatTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGc                                                                                                                                 >  1:1490248/1‑148 (MQ=255)
                   atTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCt                                                                                                                                >  1:1724766/1‑148 (MQ=255)
                    ttGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTg                                                                                                                               <  2:1545448/148‑1 (MQ=255)
                    ttGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTg                                                                                                                               <  1:2648352/148‑1 (MQ=255)
                     tGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGGTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGt                                                                                                                              <  2:2463322/148‑1 (MQ=255)
                        ggcgCGACATCCTGGCTACAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTa                                                                                                                           <  2:921848/148‑1 (MQ=255)
                        ggcgCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTa                                                                                                                           <  1:2454897/148‑1 (MQ=255)
                         gcgcGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTAc                                                                                                                          <  2:686196/148‑1 (MQ=255)
                          cgcgACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACt                                                                                                                         <  2:1569860/148‑1 (MQ=255)
                          cgcgACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACt                                                                                                                         <  1:2761277/148‑1 (MQ=255)
                            cgACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACCCTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGGTTTCTTg                                                                                                                       >  1:570377/1‑148 (MQ=255)
                             gACCTCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGt                                                                                                                      <  2:2550006/148‑1 (MQ=255)
                              aCATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGtt                                                                                                                     <  1:2007709/148‑1 (MQ=255)
                               cATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGttt                                                                                                                    <  1:2961509/148‑1 (MQ=255)
                               cATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGttt                                                                                                                    <  2:3115116/148‑1 (MQ=255)
                                 tCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtg                                                                                                                  <  1:829450/148‑1 (MQ=255)
                                 tCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTAACTTGTTTtg                                                                                                                  >  2:1191937/1‑148 (MQ=255)
                                 tCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtg                                                                                                                  <  2:3012343/148‑1 (MQ=255)
                                  ccTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtgt                                                                                                                 <  2:950048/148‑1 (MQ=255)
                                   cTGGCTAAAGCGTCGCTTATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtgtg                                                                                                                <  1:79162/148‑1 (MQ=255)
                                   cTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtgtg                                                                                                                <  1:3038439/148‑1 (MQ=255)
                                       cTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGacac                                                                                                            >  2:3110738/1‑148 (MQ=255)
                                           aGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTc                                                                                                        >  2:1155693/1‑148 (MQ=255)
                                           aGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTc                                                                                                        >  2:3231730/1‑148 (MQ=255)
                                           aGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACAGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGta                                                                                                        >  1:879801/1‑147 (MQ=255)
                                             cGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCAc                                                                                                      <  2:1651241/148‑1 (MQ=255)
                                                cGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCACTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTg                                                                                                   >  2:66521/1‑148 (MQ=255)
                                                      tatTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAggg                                                                                             <  2:483901/148‑1 (MQ=255)
                                                        ttGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAg                                                                                           <  2:2479370/148‑1 (MQ=255)
                                                            ttAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTc                                                                                       >  1:2281978/1‑148 (MQ=255)
                                                             tAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                      <  2:624990/148‑1 (MQ=255)
                                                             tAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                      <  1:1373982/148‑1 (MQ=255)
                                                             tAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                      >  2:718060/1‑148 (MQ=255)
                                                             tAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                      >  2:1525019/1‑148 (MQ=255)
                                                             tAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                      >  2:1395987/1‑148 (MQ=255)
                                                             tAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                      >  2:1396303/1‑148 (MQ=255)
                                                                 cggcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGcc                                                                                  >  1:156774/1‑148 (MQ=255)
                                                                  ggcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCa                                                                                 <  2:1544436/148‑1 (MQ=255)
                                                                  ggcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCa                                                                                 <  2:1418068/148‑1 (MQ=255)
                                                                      aGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTg                                                                             >  1:726478/1‑148 (MQ=255)
                                                                      aGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGGTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTg                                                                             >  1:2385452/1‑148 (MQ=255)
                                                                         tGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTgcag                                                                          <  2:1943381/148‑1 (MQ=255)
                                                                         tGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTgcag                                                                          <  2:2644698/148‑1 (MQ=255)
                                                                          gTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTgcagc                                                                         >  1:2774678/1‑148 (MQ=255)
                                                                            cacaCGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGcagcag                                                                       <  2:1762979/148‑1 (MQ=255)
                                                                              cacGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCagag                                                                     >  2:2253274/1‑148 (MQ=255)
                                                                              cacGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCagag                                                                     >  2:2039639/1‑148 (MQ=255)
                                                                               acGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGc                                                                    <  1:1880069/148‑1 (MQ=255)
                                                                                   cGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCa                                                                >  2:2330103/1‑148 (MQ=255)
                                                                                     tGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCaaa                                                              >  2:1462441/1‑148 (MQ=255)
                                                                                       gCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAaaca                                                            >  1:3100054/1‑148 (MQ=255)
                                                                                        cGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAaacaa                                                           <  1:2421328/148‑1 (MQ=255)
                                                                                           tATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAcgc                                                        >  2:1703942/1‑148 (MQ=255)
                                                                                            aTCTTCAATAAACTTCACCTGTTCCTCGTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAcgcg                                                       <  2:151925/148‑1 (MQ=255)
                                                                                             tCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGt                                                      >  2:2776119/1‑148 (MQ=255)
                                                                                               ttCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTAACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTAc                                                    <  2:2826315/148‑1 (MQ=255)
                                                                                                 cAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCa                                                  <  2:1490248/148‑1 (MQ=255)
                                                                                                    tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                               >  1:1685030/1‑148 (MQ=255)
                                                                                                    tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                               <  2:200682/148‑1 (MQ=255)
                                                                                                    tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                               <  2:823649/148‑1 (MQ=255)
                                                                                                      aaCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTg                                             <  1:3240092/148‑1 (MQ=255)
                                                                                                       aCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGc                                            <  1:1529866/148‑1 (MQ=255)
                                                                                                         ttCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCt                                          <  2:775392/148‑1 (MQ=255)
                                                                                                           cACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcg                                        <  1:2330103/148‑1 (MQ=255)
                                                                                                             ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgcg                                      >  1:157068/1‑148 (MQ=255)
                                                                                                             ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgcg                                      >  1:2017249/1‑148 (MQ=255)
                                                                                                             ccTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgcg                                      <  2:830949/148‑1 (MQ=255)
                                                                                                              cTGTTCCTCATCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgt                                     >  2:659487/1‑148 (MQ=255)
                                                                                                               tGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgtt                                    >  1:1093155/1‑148 (MQ=255)
                                                                                                                ggtCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttg                                   <  2:421063/146‑1 (MQ=255)
                                                                                                                  tCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCgttgtt                                 <  2:1324235/148‑1 (MQ=255)
                                                                                                                   catctTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGttt                                <  1:362139/146‑1 (MQ=255)
                                                                                                                    ctctTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGtttt                               <  1:2308015/148‑1 (MQ=255)
                                                                                                                     tctTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTa                              <  1:1872487/148‑1 (MQ=255)
                                                                                                                       ttCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCACTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAg                            <  1:66521/148‑1 (MQ=255)
                                                                                                                        tCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCGTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGc                           <  2:2551161/148‑1 (MQ=255)
                                                                                                                          ccACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                         >  2:1884637/1‑148 (MQ=255)
                                                                                                                          ccACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                         >  1:2790182/1‑148 (MQ=255)
                                                                                                                            acacTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAcc                       >  2:2529368/1‑148 (MQ=255)
                                                                                                                                 gCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCaa                  >  2:206298/1‑148 (MQ=255)
                                                                                                                                  cATCCAGCGCCGGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTGCCCGCCAGTTGCAGAAGAGCATCAAACAACGCGTCCCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAg                 <  2:1130982/148‑1 (MQ=255)
                                                                                                                                  cATCAAGCGCCGGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAg                 <  2:2782762/148‑1 (MQ=255)
                                                                                                                                  cATCAAGCGCCGGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAg                 <  2:2782756/148‑1 (MQ=255)
                                                                                                                                    tCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCa               >  2:2595469/1‑148 (MQ=255)
                                                                                                                                    tCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCa               >  1:132004/1‑148 (MQ=255)
                                                                                                                                    tCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCa               <  2:3236483/148‑1 (MQ=255)
                                                                                                                                     cacgcgcCGGGGGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGAAGAGCATAAAACAACGCGTACCACGGTGCCTCGCTTTTTTTTAAGCGACCGCCAAGCAg              <  1:1468389/145‑1 (MQ=255)
                                                                                                                                             tGGCTGGTGTCATCGCGCTCTTCGCTGTTTAATTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa      >  2:2138253/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATCGG  >  CP009273/1530884‑1531176

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: