Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,531,029 G→T 84.3% 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.3% 580.9 / 123.6 277R180S (CGC→AGC) narWnitrate reductase 2 (NRZ), delta subunit (assembly subunit)
Reads supporting (aligned to +/- strand):  ref base G (23/21);  new base T (124/109);  total (147/130)
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.67e-01

Only 100 of 287 total aligned reads displayed.

CTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATC  >  CP009273/1530882‑1531174
                                                                                                                                                   |                                                                                                                                                 
cTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCgg                                                                                                                                                  >  1:301291/1‑149 (MQ=255)
 tCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCctggctgg                                                                                                                                                 <  1:549708/149‑1 (MQ=255)
 tCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCggg                                                                                                                                                 <  1:3121913/149‑1 (MQ=255)
  cccGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGgt                                                                                                                                                <  1:2537899/149‑1 (MQ=255)
     gCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTc                                                                                                                                             <  1:2102763/149‑1 (MQ=255)
             gTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGct                                                                                                                                     >  2:1219796/1‑149 (MQ=255)
              tCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGctc                                                                                                                                    >  1:4531415/1‑149 (MQ=255)
              tCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGctc                                                                                                                                    <  1:4374905/149‑1 (MQ=255)
                 aCATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTc                                                                                                                                 >  2:4661642/1‑149 (MQ=255)
                 aCATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTc                                                                                                                                 >  1:1987334/1‑149 (MQ=255)
                  cATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCg                                                                                                                                >  1:1060241/1‑149 (MQ=255)
                    tatTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCt                                                                                                                              <  2:1283505/149‑1 (MQ=255)
                    tatTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCt                                                                                                                              <  1:967286/149‑1 (MQ=255)
                     atTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTg                                                                                                                             >  1:3826719/1‑149 (MQ=255)
                       tGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGtt                                                                                                                           >  2:1962514/1‑149 (MQ=255)
                       tGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGtt                                                                                                                           <  2:2966399/149‑1 (MQ=255)
                       tGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGtt                                                                                                                           <  1:3372664/149‑1 (MQ=255)
                           gcgcGACATCCTGGCTACAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGAGCTCTTCGCTGTTTACt                                                                                                                       <  2:834812/149‑1 (MQ=255)
                           gcgcGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACt                                                                                                                       >  2:4329063/1‑149 (MQ=255)
                           gcgcGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACt                                                                                                                       >  2:606891/1‑149 (MQ=255)
                             gcgACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTg                                                                                                                     >  2:1224701/1‑149 (MQ=255)
                              cgACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGt                                                                                                                    <  2:2262059/149‑1 (MQ=255)
                               gACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGtt                                                                                                                   >  1:917360/1‑149 (MQ=255)
                                 cATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGtttt                                                                                                                 <  2:1017376/149‑1 (MQ=255)
                                  aTCCTGGCTAAAGCGTCGCTGATATTGATTAATCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCAACAGTTCCTCTTCCCACAATGCATCAAGCGCCTGGCGGGTGTCATAGAGCTCTTCGCTGTTTACATGTTTtg                                                                                                                >  1:3795060/1‑149 (MQ=255)
                                  aTCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtg                                                                                                                >  2:3146238/1‑149 (MQ=255)
                                    ccTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTtgtg                                                                                                              <  1:1019470/149‑1 (MQ=255)
                                     cTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGa                                                                                                             <  1:4315860/149‑1 (MQ=255)
                                      tGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGac                                                                                                            >  2:806901/1‑149 (MQ=255)
                                         cTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACt                                                                                                         <  2:499731/149‑1 (MQ=255)
                                          tAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTg                                                                                                        >  2:4252902/1‑149 (MQ=255)
                                          tAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTg                                                                                                        >  2:4407690/1‑149 (MQ=255)
                                          tAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTg                                                                                                        >  1:758342/1‑149 (MQ=255)
                                                gTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGAGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACtt                                                                                                  >  2:247221/1‑149 (MQ=255)
                                                 tCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTg                                                                                                 >  2:582702/1‑149 (MQ=255)
                                                 tCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTg                                                                                                 <  2:3070760/149‑1 (MQ=255)
                                                 tCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTg                                                                                                 <  2:3468974/149‑1 (MQ=255)
                                                   gCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGaa                                                                                               >  1:1755611/1‑149 (MQ=255)
                                                   gCTGATAGTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCCAGAAACTTCACCTGTTCCTCTTCCCCCACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTGACTTGTTTTGTGACACTGTCACTTGaa                                                                                               <  2:734007/149‑1 (MQ=255)
                                                    cTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGaaa                                                                                              <  2:3894422/149‑1 (MQ=255)
                                                      gATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAgg                                                                                            <  1:3521087/149‑1 (MQ=255)
                                                      gATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAgg                                                                                            <  1:3565868/149‑1 (MQ=255)
                                                       atatTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAggg                                                                                           >  1:1653671/1‑149 (MQ=255)
                                                        tatTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGa                                                                                          <  2:3181612/149‑1 (MQ=255)
                                                          ttGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGc                                                                                        <  2:603460/149‑1 (MQ=255)
                                                          ttGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGc                                                                                        >  1:4309625/1‑149 (MQ=255)
                                                             aTTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGGTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTc                                                                                     <  2:4169767/149‑1 (MQ=255)
                                                              ttAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                    <  1:3792339/149‑1 (MQ=255)
                                                              ttAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                    <  2:4531166/149‑1 (MQ=255)
                                                              ttAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTcc                                                                                    >  2:1674945/1‑149 (MQ=255)
                                                                  gcggcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGcc                                                                                <  2:1630095/149‑1 (MQ=255)
                                                                  gcggcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGcc                                                                                <  2:3390912/149‑1 (MQ=255)
                                                                    ggcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAg                                                                              >  2:3399804/1‑149 (MQ=255)
                                                                     gcgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGt                                                                             >  2:2159963/1‑149 (MQ=255)
                                                                      cgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGtt                                                                            <  1:4441612/149‑1 (MQ=255)
                                                                      cgAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGtt                                                                            <  2:1465831/149‑1 (MQ=255)
                                                                                 acGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCa                                                                 <  2:2833311/149‑1 (MQ=255)
                                                                                 acGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCa                                                                 >  2:1772497/1‑149 (MQ=255)
                                                                                 acGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCa                                                                 >  2:1809272/1‑149 (MQ=255)
                                                                                 acGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCa                                                                 >  1:2924436/1‑149 (MQ=255)
                                                                                 acGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCa                                                                 >  2:2663510/1‑149 (MQ=255)
                                                                                     cGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCaa                                                             <  1:121564/149‑1 (MQ=255)
                                                                                       tGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAaac                                                           >  1:2825799/1‑149 (MQ=255)
                                                                                            ttATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAAcgc                                                      <  2:2240902/149‑1 (MQ=255)
                                                                                              aTCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATTGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGt                                                    >  1:1009456/1‑149 (MQ=255)
                                                                                               tCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTa                                                   <  2:3797531/149‑1 (MQ=255)
                                                                                                     aTAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACgg                                             >  1:4427480/1‑149 (MQ=255)
                                                                                                      tAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGt                                            >  2:3515473/1‑149 (MQ=255)
                                                                                                       aaaCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTg                                           >  2:3837133/1‑149 (MQ=255)
                                                                                                        aaCTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGc                                          >  1:3319136/1‑149 (MQ=255)
                                                                                                          cTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCt                                        <  1:1402910/149‑1 (MQ=255)
                                                                                                          cTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCt                                        >  1:3678355/1‑149 (MQ=255)
                                                                                                           ttCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTc                                       <  2:2363837/149‑1 (MQ=255)
                                                                                                             cACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgc                                     <  1:623582/149‑1 (MQ=255)
                                                                                                              aCCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGAGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTcgcg                                    >  2:136046/1‑149 (MQ=255)
                                                                                                                     ccTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGtttt                             >  1:342734/1‑149 (MQ=255)
                                                                                                                     ccTCTTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGtttt                             >  1:59235/1‑149 (MQ=255)
                                                                                                                     ccTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGAGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGtttt                             >  1:2211792/1‑149 (MQ=255)
                                                                                                                      ctctTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTa                            >  1:2421197/1‑149 (MQ=255)
                                                                                                                      ctctTCCCACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTa                            >  1:4256262/1‑149 (MQ=255)
                                                                                                                       tctTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGAGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTaa                           >  1:3346016/1‑149 (MQ=255)
                                                                                                                           cccACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGa                       <  1:502452/149‑1 (MQ=255)
                                                                                                                            ccACACTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAc                      >  1:3183281/1‑149 (MQ=255)
                                                                                                                             cacaCTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGAcc                     <  1:1598552/149‑1 (MQ=255)
                                                                                                                               cacTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGc                   <  2:2896240/149‑1 (MQ=255)
                                                                                                                                acTGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGCGCTTCCCGCCAGTTGCAGCAGAGCATCAAACACCGCGTACCACGGTGCCTCGCGGTGTTGTAAGCGACCGcc                  >  1:745075/1‑149 (MQ=255)
                                                                                                                                acTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGcc                  <  2:3252780/149‑1 (MQ=255)
                                                                                                                                  tGCATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCaa                >  2:361359/1‑149 (MQ=255)
                                                                                                                                    cATCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCCGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGc              >  2:2874006/1‑149 (MQ=255)
                                                                                                                                     aTCAAGCGCCTGGGTGGTGTAATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTACCCGCCAGTTGCAGCAGAGCAACAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCa             <  1:1291545/149‑1 (MQ=255)
                                                                                                                                      tCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg            <  2:2623391/149‑1 (MQ=255)
                                                                                                                                      tCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg            <  2:4691110/149‑1 (MQ=255)
                                                                                                                                      tCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg            >  1:4168500/1‑149 (MQ=255)
                                                                                                                                      tCAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAg            <  1:2563193/149‑1 (MQ=255)
                                                                                                                                       cAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTTTTTACTTGTTTTGTGAAACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAgg           <  1:4521439/149‑1 (MQ=255)
                                                                                                                                       cAAGCGCCTGGCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAgg           >  1:291922/1‑149 (MQ=255)
                                                                                                                                         aGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGc         >  1:1115298/1‑149 (MQ=255)
                                                                                                                                              ctggctggTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa    >  2:4288267/1‑149 (MQ=255)
                                                                                                                                              ctggctggTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACGGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGa    <  1:2607406/149‑1 (MQ=255)
                                                                                                                                                ggCTGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAACGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATc  <  2:4188886/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CTCCCGCACTGATGTCGACATATTGCGGCGCGACATCCTGGCTAAAGCGTCGCTGATATTGATTAAGCGGCGAGCTGTCACACGCCGTGGCGTTATCTTCAATAAACTTCACCTGTTCCTCTTCCCACACTGCATCAAGCGCCTGGCGGGTGTCATCGCGCTCTTCGCTGTTTACTTGTTTTGTGACACTGTCACTTGAAAGGGAGCTTCCCGCCAGTTGCAGCAGAGCATCAAACAACGCGTACCACGGTGCCTCGCGTTGTTTTAAGCGACCGCCAAGCAGGGCGAGGATC  >  CP009273/1530882‑1531174

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: