Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,498,013 C→A 5.9% R14S (CGT→AGT)  insQ → IS609 transposase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,498,0130CA 5.9% 1277.6 / 34.4 338R14S (CGT→AGT) insQIS609 transposase B
Reads supporting (aligned to +/- strand):  ref base C (153/165);  new base A (7/13);  total (160/178)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 344 total aligned reads displayed.

CACGCCGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGCAACAGTCACTG  >  CP009273/1497864‑1498162
                                                                                                                                                     |                                                                                                                                                     
cACGCCGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCAGTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCATTc                                                                                                                                                       >  2:2554752/1‑150 (MQ=255)
  cGCCGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGt                                                                                                                                                     >  2:2054534/1‑150 (MQ=255)
  cGCCGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGt                                                                                                                                                     <  1:1407115/150‑1 (MQ=255)
  cGCCGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGt                                                                                                                                                     <  1:1410860/150‑1 (MQ=255)
    ccGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGc                                                                                                                                                   >  2:1410756/1‑150 (MQ=255)
    ccGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGc                                                                                                                                                   >  2:1251806/1‑150 (MQ=255)
      ggATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGctc                                                                                                                                                 <  1:3236919/150‑1 (MQ=255)
       gATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGctct                                                                                                                                                <  2:852468/150‑1 (MQ=255)
        atatCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTg                                                                                                                                               >  1:142498/1‑150 (MQ=255)
          atCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGc                                                                                                                                             >  1:2785487/1‑150 (MQ=255)
              gTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGt                                                                                                                                         >  1:237213/1‑150 (MQ=255)
               tttCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTc                                                                                                                                        <  1:3471542/150‑1 (MQ=255)
                  cTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGa                                                                                                                                     >  2:3959520/1‑150 (MQ=255)
                    tttttCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAAt                                                                                                                                   <  1:442841/150‑1 (MQ=255)
                         cacaGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGaa                                                                                                                              >  2:802259/1‑150 (MQ=255)
                         cacaGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGaa                                                                                                                              >  2:3085163/1‑150 (MQ=255)
                          acaGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCAGGCACGTCAGAATGAGAAt                                                                                                                             <  1:4457677/150‑1 (MQ=255)
                           caGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATc                                                                                                                            >  2:1128272/1‑150 (MQ=255)
                              aCCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATg                                                                                                                         >  1:420374/1‑150 (MQ=255)
                                cAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAg                                                                                                                       >  2:2676427/1‑150 (MQ=255)
                                 aaaGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAgg                                                                                                                      >  1:2052875/1‑150 (MQ=255)
                                 aaaGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAgg                                                                                                                      >  1:1403405/1‑150 (MQ=255)
                                  aaGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGc                                                                                                                     >  2:3576032/1‑150 (MQ=255)
                                  aaGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGc                                                                                                                     >  2:1704423/1‑150 (MQ=255)
                                    gTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCg                                                                                                                   >  2:3341529/1‑150 (MQ=255)
                                    gTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCg                                                                                                                   >  2:210917/1‑150 (MQ=255)
                                     tATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCgg                                                                                                                  <  1:625129/150‑1 (MQ=255)
                                     tATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCgg                                                                                                                  <  1:720523/150‑1 (MQ=255)
                                      aTCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGt                                                                                                                 >  2:188442/1‑150 (MQ=255)
                                      aTCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGt                                                                                                                 <  2:4335513/150‑1 (MQ=255)
                                            tGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTaataaa                                                                                                           <  1:2251428/150‑1 (MQ=255)
                                            tGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTaataaa                                                                                                           <  1:2118391/150‑1 (MQ=255)
                                            tGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTaataaa                                                                                                           <  1:4114489/150‑1 (MQ=255)
                                             gACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTaataaat                                                                                                          >  2:699671/1‑150 (MQ=255)
                                             gACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTaataaat                                                                                                          <  1:1729973/150‑1 (MQ=255)
                                                cTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACa                                                                                                       <  2:2392507/150‑1 (MQ=255)
                                                 tgtgAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACAt                                                                                                      >  2:624092/1‑150 (MQ=255)
                                                  gtgAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATc                                                                                                     <  1:1813920/150‑1 (MQ=255)
                                                   tgAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATcc                                                                                                    >  1:3458862/1‑150 (MQ=255)
                                                    gAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATccc                                                                                                   >  2:3752333/1‑150 (MQ=255)
                                                    gAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATccc                                                                                                   >  1:1997387/1‑150 (MQ=255)
                                                     aaaCGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCt                                                                                                  >  2:558003/1‑150 (MQ=255)
                                                        cGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTAc                                                                                               >  2:538677/1‑150 (MQ=255)
                                                        cGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATAGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTAc                                                                                               <  2:3851028/150‑1 (MQ=255)
                                                        cGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATAGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTAc                                                                                               <  2:1590512/150‑1 (MQ=255)
                                                         gATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATAGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACg                                                                                              <  2:3870466/150‑1 (MQ=255)
                                                          aTTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACgg                                                                                             <  2:338983/150‑1 (MQ=255)
                                                           ttACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACggg                                                                                            <  1:1114276/150‑1 (MQ=255)
                                                                ggCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAAt                                                                                       <  2:2356740/150‑1 (MQ=255)
                                                                ggCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAAt                                                                                       <  2:1187388/150‑1 (MQ=255)
                                                                  cATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATgg                                                                                     <  1:4539950/150‑1 (MQ=255)
                                                                   aTTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGc                                                                                    >  1:2437669/1‑150 (MQ=255)
                                                                   aTTGAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGc                                                                                    >  1:4299749/1‑150 (MQ=255)
                                                                    tttAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCt                                                                                   <  1:183169/150‑1 (MQ=255)
                                                                      tAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTc                                                                                 <  1:920854/150‑1 (MQ=255)
                                                                            ccAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCctggct                                                                           <  2:4064646/150‑1 (MQ=255)
                                                                             cAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCctggctg                                                                          <  1:1280751/150‑1 (MQ=255)
                                                                                   aGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAg                                                                    >  2:2831610/1‑150 (MQ=255)
                                                                                   aGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAg                                                                    >  1:4248096/1‑150 (MQ=255)
                                                                                   aGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAg                                                                    >  1:3885396/1‑150 (MQ=255)
                                                                                     aCCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTg                                                                  <  1:2579372/150‑1 (MQ=255)
                                                                                      cccGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTgg                                                                 <  1:2898193/150‑1 (MQ=255)
                                                                                      cccGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTgg                                                                 >  2:3183298/1‑150 (MQ=255)
                                                                                      cccGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTgg                                                                 <  2:2624207/150‑1 (MQ=255)
                                                                                        cGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCTTGTCGTTTCGTTTTCAATCGTGCTCTGGCACGGCAGAATGAGAATCATGAGGCCGGTAATCAATACATCCCGTACGGGAACATGGCTTCCTGGCTGGTTGAGTGGaa                                                               >  2:148291/1‑150 (MQ=255)
                                                                                         ggtggtCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGaaa                                                              >  1:3582952/1‑150 (MQ=255)
                                                                                          gtggtCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGaaaa                                                             >  2:4500167/1‑150 (MQ=255)
                                                                                            ggtCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAAt                                                           >  1:219208/1‑150 (MQ=255)
                                                                                            ggtCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAAt                                                           >  1:511014/1‑150 (MQ=255)
                                                                                              tCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGc                                                         <  2:1501529/150‑1 (MQ=255)
                                                                                               cAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGcc                                                        >  2:1295483/1‑150 (MQ=255)
                                                                                                aaCAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCa                                                       <  1:1251780/150‑1 (MQ=255)
                                                                                                    ggAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATAGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGa                                                   <  2:1923965/150‑1 (MQ=255)
                                                                                                       gtgtGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAAc                                                <  1:1982362/150‑1 (MQ=255)
                                                                                                        tgtgAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACAGCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACg                                               <  1:3136480/150‑1 (MQ=255)
                                                                                                          tgAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCa                                             >  1:1538004/1‑150 (MQ=255)
                                                                                                            aaaTGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAAt                                           >  2:3082960/1‑150 (MQ=255)
                                                                                                              aTGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATAGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATgg                                         >  1:2950592/1‑150 (MQ=255)
                                                                                                                gAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCt                                       <  2:261441/150‑1 (MQ=255)
                                                                                                                  ggCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTa                                     <  2:3222410/150‑1 (MQ=255)
                                                                                                                  ggCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTa                                     <  1:1644109/150‑1 (MQ=255)
                                                                                                                  ggCGCTTCGCAGGAGCATGTCGTTTCTTTTTCAATCGTGCTCTGGCACTTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTa                                     <  2:3660154/150‑1 (MQ=255)
                                                                                                                       ttCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGAt                                <  1:1007897/150‑1 (MQ=255)
                                                                                                                        tCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGAtt                               <  1:802103/150‑1 (MQ=255)
                                                                                                                        tCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGAtt                               <  1:3085163/150‑1 (MQ=255)
                                                                                                                          gCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATtct                             >  2:3946099/1‑150 (MQ=255)
                                                                                                                              gcgcATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCt                         <  2:542370/150‑1 (MQ=255)
                                                                                                                                gcATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTca                       <  1:3909483/150‑1 (MQ=255)
                                                                                                                                   tGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAg                    >  1:3981721/1‑150 (MQ=255)
                                                                                                                                   tGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACAACCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAg                    >  1:2672337/1‑150 (MQ=255)
                                                                                                                                     tcgtttcgttTTCAATAGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGcc                  >  2:2794860/1‑150 (MQ=255)
                                                                                                                                      cgtttcgtttTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCa                 <  1:449226/150‑1 (MQ=255)
                                                                                                                                       gtttcgtttTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCAt                >  2:724360/1‑150 (MQ=255)
                                                                                                                                        tttcgtttTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGTAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCAtt               <  1:1028452/150‑1 (MQ=255)
                                                                                                                                          tcgtttTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGc             <  2:1016459/150‑1 (MQ=255)
                                                                                                                                           cgtttTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGCa            >  2:1820211/1‑150 (MQ=255)
                                                                                                                                                 cAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGCAACAGTc      >  2:2235225/1‑150 (MQ=255)
                                                                                                                                                   aTCGTGCTCTGGCAAGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGCAACAGTCAc    <  1:868070/150‑1 (MQ=255)
                                                                                                                                                    tCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGCAACAGTCACt   <  1:1208487/150‑1 (MQ=255)
                                                                                                                                                     cGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGCAACAGTCACTg  >  1:1729407/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACGCCGGATATCGGTTTCTTTTTTCACAGACCAAAGTATCATTTGACCTGTGAAACGATTACAGGCATTTAAATTCCAGTTAAGACCCGGTGGTCAACAGGAGTGTGAAATGAGGCGCTTCGCAGGCGCATGTCGTTTCGTTTTCAATCGTGCTCTGGCACGTCAGAATGAGAATCATGAGGCCGGTAATAAATACATCCCTTACGGGAAAATGGCTTCCTGGCTGGTTGAGTGGAAAAATGCCACTGAAACGCAATGGCTTAAAGATTCTCCCTCACAGCCATTGCAACAGTCACTG  >  CP009273/1497864‑1498162

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: