Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,922,164 C→A 92.5% P311P (CCG→CCT)  ptrB ← protease II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,922,1640CA92.5% 858.7 / 55.1 306P311P (CCG→CCT) ptrBprotease II
Reads supporting (aligned to +/- strand):  ref base C (13/10);  new base A (166/117);  total (179/127)
Fisher's exact test for biased strand distribution p-value = 8.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

Only 100 of 309 total aligned reads displayed.

ATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAA  >  CP009273/1922017‑1922313
                                                                                                                                                   |                                                                                                                                                     
aTAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAgg                                                                                                                                                     >  2:3814634/1‑150 (MQ=255)
 tAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTggagga                                                                                                                                                    <  1:1367358/150‑1 (MQ=255)
  aGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGaa                                                                                                                                                   >  1:705293/1‑150 (MQ=255)
   ggCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAAt                                                                                                                                                  >  2:1460726/1‑150 (MQ=255)
   ggCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAAt                                                                                                                                                  >  2:217968/1‑150 (MQ=255)
    gCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAAtt                                                                                                                                                 >  2:3361520/1‑150 (MQ=255)
       ggATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAAc                                                                                                                                              >  1:2637305/1‑150 (MQ=255)
        gATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAAct                                                                                                                                             <  1:333343/150‑1 (MQ=255)
        gATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAAct                                                                                                                                             <  2:2749634/150‑1 (MQ=255)
         atcatcAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAActc                                                                                                                                            >  1:3415244/1‑150 (MQ=255)
            atcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTc                                                                                                                                         >  2:2406024/1‑150 (MQ=255)
            atcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTc                                                                                                                                         >  1:2147474/1‑150 (MQ=255)
             tcaAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTcc                                                                                                                                        >  1:315645/1‑150 (MQ=255)
                 aGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACt                                                                                                                                    >  1:1649992/1‑150 (MQ=255)
                    cAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAAATGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGtt                                                                                                                                 >  1:1580777/1‑150 (MQ=255)
                     aaTACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTg                                                                                                                                <  2:854820/150‑1 (MQ=255)
                       tACCAATGACGTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCt                                                                                                                              <  2:132718/150‑1 (MQ=255)
                           aaTGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGACTTAACTCTTCCCACTGTTGCtcatc                                                                                                                          >  2:5127988/1‑150 (MQ=255)
                               aCTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGc                                                                                                                      >  2:4759095/1‑150 (MQ=255)
                                  tCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATa                                                                                                                   <  2:622513/150‑1 (MQ=255)
                                   cccGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAc                                                                                                                  >  1:4491602/1‑150 (MQ=255)
                                    ccGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACg                                                                                                                 <  1:4759095/150‑1 (MQ=255)
                                     cGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACgg                                                                                                                >  2:649730/1‑150 (MQ=255)
                                     cGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACgg                                                                                                                >  1:919748/1‑150 (MQ=255)
                                       ggTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGt                                                                                                              >  2:1167699/1‑150 (MQ=255)
                                         tCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTAc                                                                                                            >  1:1006904/1‑150 (MQ=255)
                                           ttGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACgg                                                                                                          <  1:2193128/150‑1 (MQ=255)
                                           ttGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACgg                                                                                                          <  2:3181144/150‑1 (MQ=255)
                                            tGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGt                                                                                                         >  1:252399/1‑150 (MQ=255)
                                              cGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGtat                                                                                                       <  2:3086023/150‑1 (MQ=255)
                                              cGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGtat                                                                                                       >  2:3594410/1‑150 (MQ=255)
                                               ggTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGtata                                                                                                      >  2:3293081/1‑150 (MQ=255)
                                               ggTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGtata                                                                                                      <  1:3051141/150‑1 (MQ=255)
                                                 ttAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAg                                                                                                    >  1:4407670/1‑150 (MQ=255)
                                                   aaTTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGcc                                                                                                  <  2:3710442/150‑1 (MQ=255)
                                                    aTTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCa                                                                                                 >  1:4958753/1‑150 (MQ=255)
                                                      ttGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCaaa                                                                                               >  2:1554043/1‑150 (MQ=255)
                                                        ggCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTATGGTATAAGCCAAAGt                                                                                             <  2:4497640/150‑1 (MQ=255)
                                                           gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGtttt                                                                                          <  2:4871538/150‑1 (MQ=255)
                                                           gcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGtttt                                                                                          <  1:4167710/150‑1 (MQ=255)
                                                              aaCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGc                                                                                       <  2:3436680/150‑1 (MQ=255)
                                                              aaCTGGTTAACCCGCGCTGACGATCTAAAACCACCAGCCAGTCGGTAAACAGCTTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGc                                                                                       >  1:4049962/1‑150 (MQ=255)
                                                               aCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGcc                                                                                      <  1:1296609/150‑1 (MQ=255)
                                                                  tgTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTg                                                                                   >  1:2259724/2‑150 (MQ=255)
                                                                  ggTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTg                                                                                   <  2:2584046/150‑1 (MQ=255)
                                                                    ttAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGc                                                                                 >  2:4168839/1‑150 (MQ=255)
                                                                    ttAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGc                                                                                 <  2:3099499/150‑1 (MQ=255)
                                                                        cccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGtt                                                                             <  2:1827100/150‑1 (MQ=255)
                                                                         ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTg                                                                            <  2:2515025/150‑1 (MQ=255)
                                                                          cgcgCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTgg                                                                           <  2:2359231/150‑1 (MQ=255)
                                                                            cgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                         <  1:1980546/150‑1 (MQ=255)
                                                                            cgcTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGaa                                                                         >  1:1210141/1‑150 (MQ=255)
                                                                               tGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGc                                                                      >  1:2950938/1‑150 (MQ=255)
                                                                                   gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                                  <  2:3071745/150‑1 (MQ=255)
                                                                                   gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                                  >  1:1537860/1‑150 (MQ=255)
                                                                                   gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGAt                                                                  <  1:283876/150‑1 (MQ=255)
                                                                                     tctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaa                                                                <  2:1137318/150‑1 (MQ=255)
                                                                                      ctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaa                                                               >  1:4206524/1‑150 (MQ=255)
                                                                                      ctTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaa                                                               >  1:1920126/1‑150 (MQ=255)
                                                                                         cAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                            >  1:550489/1‑150 (MQ=255)
                                                                                          aaCCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACg                                                           >  1:3619004/1‑150 (MQ=255)
                                                                                           accaccAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGa                                                          >  1:4269209/1‑150 (MQ=255)
                                                                                             caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATg                                                        <  1:5139238/150‑1 (MQ=255)
                                                                                               ccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCt                                                      <  2:91784/150‑1 (MQ=255)
                                                                                               ccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCt                                                      >  2:3842454/1‑150 (MQ=255)
                                                                                               ccagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCt                                                      >  2:197020/1‑150 (MQ=255)
                                                                                                cagctagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTg                                                     >  2:3890395/1‑150 (MQ=255)
                                                                                                cagccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTg                                                     >  1:1936997/1‑150 (MQ=255)
                                                                                                 tgccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTgg                                                    >  2:1087430/2‑150 (MQ=255)
                                                                                                 agccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTgg                                                    >  2:547933/1‑150 (MQ=255)
                                                                                                 agccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTgg                                                    <  2:295435/150‑1 (MQ=255)
                                                                                                  gccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACTGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGt                                                   <  1:3871989/150‑1 (MQ=255)
                                                                                                   ccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTa                                                  >  1:3440897/1‑150 (MQ=255)
                                                                                                       tCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGa                                              >  2:1011269/1‑150 (MQ=255)
                                                                                                          gTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCa                                           >  1:667872/1‑150 (MQ=255)
                                                                                                          gTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCa                                           <  1:3814634/150‑1 (MQ=255)
                                                                                                           tAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCaa                                          >  2:3551323/1‑150 (MQ=255)
                                                                                                            aaaCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAg                                         >  2:2216198/1‑150 (MQ=255)
                                                                                                              aCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGc                                       >  1:4749845/1‑150 (MQ=255)
                                                                                                                aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTg                                     <  2:214739/150‑1 (MQ=255)
                                                                                                                aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTg                                     >  1:1892496/1‑150 (MQ=255)
                                                                                                                aGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTg                                     <  1:2293318/150‑1 (MQ=255)
                                                                                                                     aaaCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                                >  2:3955451/1‑150 (MQ=255)
                                                                                                                     aaaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTc                                >  1:3939788/1‑150 (MQ=255)
                                                                                                                      aaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCg                               >  2:2337700/1‑150 (MQ=255)
                                                                                                                      aaCCCTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCg                               >  2:1256127/1‑150 (MQ=255)
                                                                                                                         ccTTCCAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGa                            >  1:3257707/1‑150 (MQ=255)
                                                                                                                             cccgCATGATGTGTTCGCGTGGAGGACTTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCtt                        <  2:4629106/147‑1 (MQ=255)
                                                                                                                             ccAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCtt                        <  1:2923280/150‑1 (MQ=255)
                                                                                                                              cAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                       <  2:1409750/150‑1 (MQ=255)
                                                                                                                              cAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                       <  2:2474478/150‑1 (MQ=255)
                                                                                                                              cAGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCttt                       <  1:2398464/150‑1 (MQ=255)
                                                                                                                               aGCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTg                      >  1:256764/1‑150 (MQ=255)
                                                                                                                                gCATGATGTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgc                     <  1:1460726/150‑1 (MQ=255)
                                                                                                                                      tgTTTTCGCGTGGCGGACTTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcg               >  2:324277/1‑150 (MQ=255)
                                                                                                                                      tgTTTTCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcg               >  2:4886443/1‑150 (MQ=255)
                                                                                                                                        ttttCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcggc             >  2:737505/1‑150 (MQ=255)
                                                                                                                                         tttCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTTCGGCGCGGCa            >  2:1884760/1‑150 (MQ=255)
                                                                                                                                         aatCGCGTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGAGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGAGGCGCGGCa            >  2:4011900/3‑150 (MQ=255)
                                                                                                                                            cgcgTGGAGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGACCGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaa         >  2:5003083/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
ATAGGCCGGATCATCAAAGGCAATACCAATGACTTCCCGGGTCTTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAA  >  CP009273/1922017‑1922313

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: ‑