Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,096,445 A→T 6.9% Q25L (CAA→CTA)  mutY → adenine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,096,4450AT 6.9% 896.7 / 40.6 276Q25L (CAA→CTA) mutYadenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base A (125/132);  new base T (4/15);  total (129/147)
Fisher's exact test for biased strand distribution p-value = 2.98e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 280 total aligned reads displayed.

CGGAAAGTTCCGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCC  >  CP009273/3096299‑3096594
                                                                                                                                                  |                                                                                                                                                     
   aaaGTTCCGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGa                                                                                                                                                 >  1:1106101/1‑150 (MQ=255)
   aaaGTTCCGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGa                                                                                                                                                 >  2:2647270/1‑150 (MQ=255)
    aaGTTCCGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGAc                                                                                                                                                <  2:1759960/150‑1 (MQ=255)
       ttCCGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAg                                                                                                                                             <  1:3080836/150‑1 (MQ=255)
         ccGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGAc                                                                                                                                           >  1:1178320/1‑150 (MQ=255)
         ccGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTAAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGAc                                                                                                                                           <  2:798746/150‑1 (MQ=255)
          cGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACg                                                                                                                                          <  1:786600/150‑1 (MQ=255)
           ggTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGc                                                                                                                                         >  1:3409903/1‑150 (MQ=255)
           ggTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGc                                                                                                                                         >  2:2776623/1‑150 (MQ=255)
           ggTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGc                                                                                                                                         >  2:1602025/1‑150 (MQ=255)
           ggTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGc                                                                                                                                         >  2:3862421/1‑150 (MQ=255)
            gTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGcc                                                                                                                                        >  1:1816032/1‑150 (MQ=255)
              ttACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCt                                                                                                                                      <  2:2830355/150‑1 (MQ=255)
                 cacCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACa                                                                                                                                   <  1:664007/150‑1 (MQ=255)
                    ccTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAg                                                                                                                                <  2:3496963/150‑1 (MQ=255)
                    ccTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAg                                                                                                                                <  1:3549000/150‑1 (MQ=255)
                     cTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGt                                                                                                                               <  1:4506993/150‑1 (MQ=255)
                       gCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCTAATTGACAAGACGCCCTACAAAGTAt                                                                                                                             <  1:2861227/150‑1 (MQ=255)
                       gCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTAt                                                                                                                             <  2:9554/150‑1 (MQ=255)
                        ccGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATg                                                                                                                            >  1:844472/1‑150 (MQ=255)
                        ccGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATg                                                                                                                            <  1:2694441/150‑1 (MQ=255)
                          gTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGc                                                                                                                          <  1:3716806/150‑1 (MQ=255)
                          gTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGc                                                                                                                          <  1:846695/150‑1 (MQ=255)
                          gTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGc                                                                                                                          <  2:4287843/150‑1 (MQ=255)
                              cTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCtctc                                                                                                                      >  2:2036077/1‑150 (MQ=255)
                              cTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCtctc                                                                                                                      <  2:253888/150‑1 (MQ=255)
                               tgtgCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCTAATTGACAAGACGCCCTACAAAGTATGGCTCTCa                                                                                                                     <  1:2176428/150‑1 (MQ=255)
                               tctgctgcAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCTAATTGACAAGACGCCCTACAAAGTATGGCTCTCa                                                                                                                     <  1:4415319/148‑1 (MQ=255)
                                       aaTCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGAt                                                                                                             >  1:551270/1‑150 (MQ=255)
                                        aTCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATg                                                                                                            <  1:1294698/150‑1 (MQ=255)
                                         tCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCTAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGt                                                                                                           <  1:913355/150‑1 (MQ=255)
                                         tCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGt                                                                                                           >  1:1004307/1‑150 (MQ=255)
                                         tCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGt                                                                                                           <  2:241155/150‑1 (MQ=255)
                                         tCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGt                                                                                                           >  1:1889779/1‑150 (MQ=255)
                                           ttGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTg                                                                                                         <  2:454957/150‑1 (MQ=255)
                                            tGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGc                                                                                                        <  2:687898/150‑1 (MQ=255)
                                            tGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGc                                                                                                        <  2:4221826/150‑1 (MQ=255)
                                                    acaacaGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACt                                                                                                <  2:2686093/150‑1 (MQ=255)
                                                     caacaGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACTAACTc                                                                                               <  1:3489415/150‑1 (MQ=255)
                                                       acaGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAg                                                                                             >  1:3457292/1‑150 (MQ=255)
                                                          gTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGtt                                                                                          >  1:2668944/1‑150 (MQ=255)
                                                           tGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTg                                                                                         >  1:3797162/1‑150 (MQ=255)
                                                           tGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTg                                                                                         <  1:3895336/150‑1 (MQ=255)
                                                             aaTTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCg                                                                                       >  2:2384810/1‑150 (MQ=255)
                                                              aTTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGa                                                                                      <  2:13652/150‑1 (MQ=255)
                                                              aTTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGa                                                                                      <  2:3723156/150‑1 (MQ=255)
                                                                tCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGAcc                                                                                    <  2:1710280/150‑1 (MQ=255)
                                                                  ggTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGt                                                                                  >  1:2518458/1‑150 (MQ=255)
                                                                   gTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGtt                                                                                 <  1:1673407/150‑1 (MQ=255)
                                                                     gACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTAt                                                                               <  2:511428/150‑1 (MQ=255)
                                                                      aCCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATc                                                                              <  2:1363857/150‑1 (MQ=255)
                                                                       ccATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCTAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATcc                                                                             <  1:2611613/150‑1 (MQ=255)
                                                                       ccATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATcc                                                                             <  1:1576847/150‑1 (MQ=255)
                                                                            cAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTAt                                                                        <  2:1106101/150‑1 (MQ=255)
                                                                             aaGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTAtt                                                                       <  1:3619808/150‑1 (MQ=255)
                                                                              aGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTAttt                                                                      <  1:1605673/150‑1 (MQ=255)
                                                                                cGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCGGGTTGCGACCGTTATCCCCTATTTTg                                                                    <  2:2292392/150‑1 (MQ=255)
                                                                                cGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTg                                                                    >  1:4255251/1‑150 (MQ=255)
                                                                                 gTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACACACTCAGGTTGCGACCGTTATCCCCTATTTTGa                                                                   >  2:1788307/1‑150 (MQ=255)
                                                                                 gTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGa                                                                   >  2:44101/1‑150 (MQ=255)
                                                                                 ctCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGa                                                                   >  2:1786825/2‑150 (MQ=255)
                                                                                    gCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACg                                                                >  2:1861362/1‑150 (MQ=255)
                                                                                        ttttCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCttt                                                            >  2:4156140/1‑150 (MQ=255)
                                                                                          ttCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTAt                                                          >  1:2492559/1‑150 (MQ=255)
                                                                                          ttCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTAt                                                          <  2:4102519/150‑1 (MQ=255)
                                                                                             aGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGgc                                                       >  2:1423272/1‑150 (MQ=255)
                                                                                                    gTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTc                                                >  2:3650357/1‑150 (MQ=255)
                                                                                                    gTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTc                                                >  2:4353986/1‑150 (MQ=255)
                                                                                                     ttCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTcc                                               <  1:334042/150‑1 (MQ=255)
                                                                                                     ttCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTcc                                               <  1:1589558/150‑1 (MQ=255)
                                                                                                       cTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACACAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCg                                             <  2:1536684/150‑1 (MQ=255)
                                                                                                           aCTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCTAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATTTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACgg                                         <  1:14540/150‑1 (MQ=255)
                                                                                                            cTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGt                                        <  2:1466209/150‑1 (MQ=255)
                                                                                                            cTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGt                                        <  2:2977797/150‑1 (MQ=255)
                                                                                                             tGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTg                                       <  1:1566634/150‑1 (MQ=255)
                                                                                                              ggTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGa                                      >  2:2706382/1‑150 (MQ=255)
                                                                                                              ggTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGa                                      >  2:3191711/1‑150 (MQ=255)
                                                                                                               gTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGAc                                     <  1:1009196/150‑1 (MQ=255)
                                                                                                                   gATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAt                                 >  2:218796/1‑150 (MQ=255)
                                                                                                                   gATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAt                                 >  2:1969558/1‑150 (MQ=255)
                                                                                                                   gATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAt                                 >  2:168772/1‑150 (MQ=255)
                                                                                                                     tAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtct                               <  1:1989162/150‑1 (MQ=255)
                                                                                                                     tAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtct                               <  1:2230633/150‑1 (MQ=255)
                                                                                                                      aaaTACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGAtctc                              <  2:451264/150‑1 (MQ=255)
                                                                                                                             ggCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATg                       <  2:4293600/150‑1 (MQ=255)
                                                                                                                               cGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATgcg                     >  1:2489435/1‑150 (MQ=255)
                                                                                                                                  aaaaCTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGcgccg                  >  1:3052657/1‑150 (MQ=255)
                                                                                                                                   aaaCTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGcgccgc                 >  2:2086139/1‑150 (MQ=255)
                                                                                                                                   aaaCTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGcgccgc                 >  2:562355/1‑150 (MQ=255)
                                                                                                                                    aaCTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCt                >  2:4035065/1‑150 (MQ=255)
                                                                                                                                     aCTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTc               >  2:1036812/1‑150 (MQ=255)
                                                                                                                                       tctGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTcga             >  1:1052146/1‑150 (MQ=255)
                                                                                                                                         tGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTcgacg           <  1:2647270/150‑1 (MQ=255)
                                                                                                                                         tGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGATTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTcgacg           >  1:588969/1‑150 (MQ=255)
                                                                                                                                            ccTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAg        >  2:1275107/1‑150 (MQ=255)
                                                                                                                                              tGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGAGATGTTGCACCAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGtt      >  2:2902735/1‑150 (MQ=255)
                                                                                                                                               ggCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTtc     >  1:844553/1‑150 (MQ=255)
                                                                                                                                                gCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTtct    <  1:27929/150‑1 (MQ=255)
                                                                                                                                                 cAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTGTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCCATGCGCCGCTCGACGAAGTtctc   >  2:3462081/1‑150 (MQ=255)
                                                                                                                                                  aaaTTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTcc  >  2:2830228/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGAAAGTTCCGGTTTACACCCTGCCGTCGCTGTGCTGCAATCTTGCCCCCAACAACAGTGAATTCGGTGACCATGCAAGCGTCGCAATTTTCAGCCCAGGTTCTGGACTGGTACGATAAATACGGGCGAAAAACTCTGCCCTGGCAAATTGACAAGACGCCCTACAAAGTATGGCTCTCAGAAGTGATGTTGCAACAAACTCAGGTTGCGACCGTTATCCCCTATTTTGAACGCTTTATGGCGCGCTTCCCGACGGTGACCGATCTCGCCAATGCGCCGCTCGACGAAGTTCTCC  >  CP009273/3096299‑3096594

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