Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,825,702 C→A 70.6% K732N (AAG→AAT yicI ← putative alpha‑glucosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,825,7020CA70.6% 479.8 / inf 333K732N (AAG→AATyicIputative alpha‑glucosidase
Reads supporting (aligned to +/- strand):  ref base C (45/53);  new base A (132/103);  total (177/156)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

GACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTCG  >  CP009273/3825553‑3825849
                                                                                                                                                     |                                                                                                                                                   
  cggaCCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGttatt                                                                                                                                                   <  2:3429656/147‑1 (MQ=255)
    tACCCGGCATTCAATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTtcg                                                                                                                                                 <  1:2631539/150‑1 (MQ=255)
     aCCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgc                                                                                                                                                <  2:4405412/150‑1 (MQ=255)
       ccGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcct                                                                                                                                              <  2:3043902/150‑1 (MQ=255)
         ggCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcg                                                                                                                                            >  2:4765671/1‑150 (MQ=255)
         ggCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcg                                                                                                                                            >  2:1237570/1‑150 (MQ=255)
         ggCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcg                                                                                                                                            >  1:2034123/1‑150 (MQ=255)
           cATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATtcgcctcgcc                                                                                                                                          >  2:1019766/1‑150 (MQ=255)
            aTTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGccc                                                                                                                                         <  1:3796660/150‑1 (MQ=255)
             ttCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCg                                                                                                                                        <  2:679344/150‑1 (MQ=255)
                cATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCAc                                                                                                                                     >  1:424880/1‑150 (MQ=255)
                 aTCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCAcc                                                                                                                                    <  1:522352/150‑1 (MQ=255)
                  tCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCa                                                                                                                                   >  2:1768982/1‑150 (MQ=255)
                       aGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTcac                                                                                                                              >  1:1061784/1‑150 (MQ=255)
                       aGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTcac                                                                                                                              >  2:2260188/1‑150 (MQ=255)
                           tACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTa                                                                                                                          <  1:210141/150‑1 (MQ=255)
                             cAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAAt                                                                                                                        <  1:4414626/150‑1 (MQ=255)
                              aaCGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATc                                                                                                                       <  1:4794726/150‑1 (MQ=255)
                                cGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCgt                                                                                                                     >  1:44816/1‑150 (MQ=255)
                                cGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCgt                                                                                                                     >  1:3380397/1‑150 (MQ=255)
                                 gTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCgtg                                                                                                                    >  1:1352246/1‑150 (MQ=255)
                                  tAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCgtgt                                                                                                                   <  2:1255378/150‑1 (MQ=255)
                                    aTTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTg                                                                                                                 >  2:753752/1‑150 (MQ=255)
                                     ttGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGc                                                                                                                <  1:896767/150‑1 (MQ=255)
                                      tGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGcc                                                                                                               >  1:1569228/1‑150 (MQ=255)
                                         cAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGt                                                                                                            <  2:3946451/150‑1 (MQ=255)
                                             gcATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGt                                                                                                        >  1:2073905/1‑150 (MQ=255)
                                               aTTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgc                                                                                                      >  1:3932919/1‑150 (MQ=255)
                                                ttCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgct                                                                                                     <  1:3217623/150‑1 (MQ=255)
                                                 tCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgctg                                                                                                    <  2:1243225/150‑1 (MQ=255)
                                                  cccTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGtgctgc                                                                                                   >  2:19714/1‑150 (MQ=255)
                                                   ccTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTgctgct                                                                                                  >  2:3657412/1‑150 (MQ=255)
                                                     ttGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCttt                                                                                                >  2:1436819/1‑150 (MQ=255)
                                                     ttGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCttt                                                                                                >  1:1726788/1‑150 (MQ=255)
                                                       gAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTta                                                                                              <  2:4123684/150‑1 (MQ=255)
                                                       gAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTGTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTta                                                                                              <  1:4581302/150‑1 (MQ=255)
                                                        aGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaa                                                                                             >  2:4489975/1‑150 (MQ=255)
                                                         ggCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTtaaa                                                                                            >  1:672986/1‑150 (MQ=255)
                                                                 caccaGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGAt                                                                                    >  1:775364/1‑150 (MQ=255)
                                                                 caccaGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGAt                                                                                    >  1:1657683/1‑150 (MQ=255)
                                                                     aGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCAcc                                                                                >  2:3152595/1‑150 (MQ=255)
                                                                       ccccTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGa                                                                              <  1:3922077/150‑1 (MQ=255)
                                                                         ccTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATc                                                                            >  2:2535183/1‑150 (MQ=255)
                                                                           tGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCg                                                                          <  2:4779010/150‑1 (MQ=255)
                                                                           tGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCg                                                                          >  1:2206344/1‑150 (MQ=255)
                                                                             cTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTc                                                                        <  1:3053852/150‑1 (MQ=255)
                                                                              tCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTca                                                                       <  1:3643096/150‑1 (MQ=255)
                                                                              tCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTca                                                                       <  1:3343551/150‑1 (MQ=255)
                                                                               cACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTcag                                                                      >  2:573963/1‑150 (MQ=255)
                                                                                 cTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTcagca                                                                    <  1:817466/150‑1 (MQ=255)
                                                                                  tttCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTcagcag                                                                   >  2:4086466/1‑150 (MQ=255)
                                                                                  tttCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTcagcag                                                                   >  2:1627227/1‑150 (MQ=255)
                                                                                    tCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCg                                                                 <  1:1037865/150‑1 (MQ=255)
                                                                                     cAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCgg                                                                >  1:1156574/1‑150 (MQ=255)
                                                                                        ccTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCa                                                             >  2:1271842/1‑150 (MQ=255)
                                                                                        ccTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCa                                                             >  2:3734142/1‑150 (MQ=255)
                                                                                           gCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACtt                                                          <  2:2674305/150‑1 (MQ=255)
                                                                                              aaCCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTcac                                                       >  2:1871055/1‑150 (MQ=255)
                                                                                               aCCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTcaca                                                      <  1:1445611/150‑1 (MQ=255)
                                                                                               aCCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTcaca                                                      >  1:3873399/1‑150 (MQ=255)
                                                                                                 cGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGa                                                    >  2:4709828/1‑150 (MQ=255)
                                                                                                 cGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGa                                                    >  2:1698758/1‑150 (MQ=255)
                                                                                                  gTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGAc                                                   <  2:1041469/150‑1 (MQ=255)
                                                                                                  gTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGAc                                                   >  1:1168536/1‑150 (MQ=255)
                                                                                                   tCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACg                                                  <  1:1266858/150‑1 (MQ=255)
                                                                                                   tCTTGCAGACCATTGACTGTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACg                                                  <  1:4737145/150‑1 (MQ=255)
                                                                                                    cTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACgg                                                 <  2:2652727/150‑1 (MQ=255)
                                                                                                      tGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCt                                               <  2:4205436/150‑1 (MQ=255)
                                                                                                      tGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCt                                               >  1:921526/1‑150 (MQ=255)
                                                                                                          gACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCAt                                           <  1:2342319/150‑1 (MQ=255)
                                                                                                           aCCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATg                                          >  2:4145364/1‑150 (MQ=255)
                                                                                                                ttACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGt                                     >  1:2753141/1‑150 (MQ=255)
                                                                                                                   cTTTCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCtt                                  >  1:3963837/1‑150 (MQ=255)
                                                                                                                    tttCACAACATTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTg                                 <  1:1887324/150‑1 (MQ=255)
                                                                                                                      tCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCa                               >  2:2345282/1‑150 (MQ=255)
                                                                                                                        acaacaTTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGa                             <  1:2490126/150‑1 (MQ=255)
                                                                                                                         caacaTTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCGGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAt                            >  2:3573241/1‑150 (MQ=255)
                                                                                                                         caacaTTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAt                            <  2:2837755/150‑1 (MQ=255)
                                                                                                                         caacaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGAt                            <  2:3871042/150‑1 (MQ=255)
                                                                                                                           acaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTg                          >  1:1598297/1‑150 (MQ=255)
                                                                                                                           acaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTg                          <  1:4086466/150‑1 (MQ=255)
                                                                                                                           acaTTGCGCAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTg                          >  1:582437/1‑150 (MQ=255)
                                                                                                                              gtGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAag                       <  1:4651881/149‑1 (MQ=255)
                                                                                                                                  gcAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGt                   <  2:854936/150‑1 (MQ=255)
                                                                                                                                  gcAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGt                   >  2:525658/1‑150 (MQ=255)
                                                                                                                                   cAGGCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTg                  <  2:2131853/150‑1 (MQ=255)
                                                                                                                                     ggCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGa                >  2:3367723/1‑150 (MQ=255)
                                                                                                                                     ggCACAGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGa                <  2:2634224/150‑1 (MQ=255)
                                                                                                                                       cacaGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACATTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAAt              <  2:493694/150‑1 (MQ=255)
                                                                                                                                         caGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGc            >  1:2743406/1‑150 (MQ=255)
                                                                                                                                         caGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGc            <  1:288524/150‑1 (MQ=255)
                                                                                                                                         caGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGc            <  1:2716443/150‑1 (MQ=255)
                                                                                                                                          aGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGcg           <  1:995780/150‑2 (MQ=255)
                                                                                                                                          aGTGTCCAGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGcc           >  2:1592348/1‑150 (MQ=255)
                                                                                                                                              tCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGc       <  2:1339863/150‑1 (MQ=255)
                                                                                                                                                cAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCt     >  2:1462896/1‑150 (MQ=255)
                                                                                                                                                 aGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCtt    >  1:4108648/1‑150 (MQ=255)
                                                                                                                                                 aGTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCtt    <  2:1156631/150‑1 (MQ=255)
                                                                                                                                                  gTTATTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTc   >  1:4356059/1‑150 (MQ=255)
                                                                                                                                                   ttcttcGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTCg  >  2:939853/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GACGTACCCGGCATTCCATCAATAGATTACAACGTAATTGTCAGCGCATTCCCTTGAGGCTTCACCACCAGCCCCTGCTCACTTTCAGCCTGCGAACCGTCTTGCAGACCATTTACTTTCACAACATTGCGCAGGCACAGTGTCCAGTTCTTCGCCTCGCCCGCACCAGTCACAGTAATCGTGTTGCCAGTACGTGCTGCTTTTAAAGTAAAGATCACCGATCCGTCAGCAGCGGGCACTTCACAGACGGCTTCATGCCCGTCTTGCAGATTGAAGAGGTGGAATGCCGTGCCTTCG  >  CP009273/3825553‑3825849

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: