Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,077,788 C→G 6.1% R74G (CGT→GGT)  cbtA → CP4‑44 prophage; cytoskeleton‑binding toxin CbtA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,077,7880CG 6.1% 380.4 / 10.2 113R74G (CGT→GGT) cbtACP4‑44 prophage; cytoskeleton‑binding toxin CbtA
Reads supporting (aligned to +/- strand):  ref base C (49/56);  new base G (2/5);  total (52/61)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

Only 100 of 114 total aligned reads displayed.

TACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGC  >  NC_000913/2077639‑2077936
                                                                                                                                                     |                                                                                                                                                    
tACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGc                                                                                                                                                      <  2:362411/150‑1 (MQ=255)
 aCTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCg                                                                                                                                                     <  1:260281/150‑1 (MQ=255)
 aCTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCg                                                                                                                                                     >  1:1393088/1‑150 (MQ=255)
  ctgctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGt                                                                                                                                                    >  2:534468/1‑150 (MQ=255)
  ctgctgTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGt                                                                                                                                                    >  1:501057/1‑150 (MQ=255)
       gTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATCGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGa                                                                                                                                               <  1:1555806/150‑1 (MQ=255)
        tCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAc                                                                                                                                              <  2:824375/150‑1 (MQ=255)
         cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGTAGCATATTGAGGCAGTCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                             >  2:1230260/1‑150 (MQ=255)
         cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                             >  1:1047245/1‑150 (MQ=255)
         cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                             >  2:443395/1‑150 (MQ=255)
         cccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                             <  2:1497309/150‑1 (MQ=255)
          ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGAcc                                                                                                                                             >  1:470135/1‑149 (MQ=255)
          ccGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCa                                                                                                                                            >  1:965482/1‑150 (MQ=255)
             aCTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTCGGTACGCACCAGcg                                                                                                                                         <  2:1318007/150‑2 (MQ=255)
              ctgctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTCGGTACGCACCAGcgc                                                                                                                                        <  1:70111/150‑3 (MQ=255)
               tgctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTCGGTACGCACCAGcgcg                                                                                                                                       <  1:1409469/150‑4 (MQ=255)
                gctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTCGGTACGCACCAGcgcgt                                                                                                                                      <  2:1287301/150‑5 (MQ=255)
                gctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTCGGTACGCACCAGcgcgt                                                                                                                                      <  1:465841/150‑5 (MQ=255)
                 ctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGa                                                                                                                                     <  2:1288198/150‑1 (MQ=255)
                 ctgGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGa                                                                                                                                     >  2:1171848/1‑150 (MQ=255)
                   ggACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAtt                                                                                                                                   <  2:391902/150‑1 (MQ=255)
                   ggACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGAtt                                                                                                                                   <  2:386111/150‑1 (MQ=255)
                   ggAACAGCACTAAGGCCTCACACTGAATGACACCCATTTTGCCGATTAACGTGTGATTGAGCAGCATATTGAGGCAGGAATTTAACTGTGTGATGCGGTGAACTTTCTCGTGGATAAATACGCGCTGGTGGGTACCGACCAGCCGGGAtt                                                                                                                                   >  2:895466/1‑150 (MQ=255)
                    gACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTc                                                                                                                                  <  2:74304/150‑1 (MQ=255)
                        aGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAgcg                                                                                                                              >  1:1568201/1‑150 (MQ=255)
                        aGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAgcg                                                                                                                              <  2:1047245/150‑1 (MQ=255)
                           aCTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCt                                                                                                                           <  2:528708/150‑1 (MQ=255)
                            cTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTg                                                                                                                          <  2:208964/150‑1 (MQ=255)
                            cTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTACACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTg                                                                                                                          >  2:1175447/1‑150 (MQ=255)
                             tATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                                                         <  2:1262195/150‑1 (MQ=255)
                             tATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGt                                                                                                                         <  2:56667/150‑1 (MQ=255)
                              aTGGCCTCACACTGAATGACATACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGTCAGGCATTTCACTGTGTGATTCGGTGAACTTTGTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGGCTGTa                                                                                                                        <  2:1227189/150‑1 (MQ=255)
                                ggCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTAcc                                                                                                                      <  1:1248123/150‑1 (MQ=255)
                                 gCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTAccc                                                                                                                     >  2:59331/1‑150 (MQ=255)
                                  ccTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCg                                                                                                                    >  2:372703/1‑150 (MQ=255)
                                     cacaCTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGctc                                                                                                                 >  1:284669/1‑150 (MQ=255)
                                       cacTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCtctc                                                                                                               <  1:1171848/150‑1 (MQ=255)
                                       cacTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCtctc                                                                                                               >  1:1543019/1‑150 (MQ=255)
                                        acTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCa                                                                                                              <  1:1212287/150‑1 (MQ=255)
                                        acTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCa                                                                                                              <  2:210804/150‑1 (MQ=255)
                                         cTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAg                                                                                                             <  2:824830/150‑1 (MQ=255)
                                           gAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGtt                                                                                                           <  1:1290703/150‑1 (MQ=255)
                                            aaTGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTta                                                                                                          <  2:1518479/150‑1 (MQ=255)
                                             aTGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAAGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaa                                                                                                         <  2:1656186/150‑1 (MQ=255)
                                              tGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaat                                                                                                        <  1:1316516/150‑1 (MQ=255)
                                              tGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaat                                                                                                        <  1:243883/150‑1 (MQ=255)
                                               gACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTtaata                                                                                                       >  2:710795/1‑150 (MQ=255)
                                                 cacacCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATaaa                                                                                                     >  1:1190187/1‑150 (MQ=255)
                                                 cacacCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATaaa                                                                                                     >  1:595948/1‑150 (MQ=255)
                                                   cacCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACa                                                                                                   >  1:1305727/1‑150 (MQ=255)
                                                     ccTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGc                                                                                                 >  2:779405/1‑150 (MQ=255)
                                                        tttGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATc                                                                                              >  1:650499/1‑150 (MQ=255)
                                                         ttGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCg                                                                                             <  2:965482/150‑1 (MQ=255)
                                                              gATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATAtc                                                                                        >  1:1147927/1‑150 (MQ=255)
                                                                    cGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCgg                                                                                  >  2:636141/1‑150 (MQ=255)
                                                                      tgtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGc                                                                                >  1:691942/1‑150 (MQ=255)
                                                                       gtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCt                                                                               >  2:446973/1‑150 (MQ=255)
                                                                       gtgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCt                                                                               <  1:766720/150‑1 (MQ=255)
                                                                        tgATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTc                                                                              <  2:1524330/150‑1 (MQ=255)
                                                                           ttGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCa                                                                           <  2:1404657/150‑1 (MQ=255)
                                                                              agcagcATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGctct                                                                                                                >  2:1155401/1‑110 (MQ=255)
                                                                              agcagcATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGctct                                                                                                                <  1:1155401/110‑1 (MQ=255)
                                                                              agcagcATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAggg                                                                        <  1:534468/150‑1 (MQ=255)
                                                                               gcagcaTATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGc                                                                       <  1:625277/150‑1 (MQ=255)
                                                                                  gcaTATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGAc                                                                    >  1:1550718/1‑150 (MQ=255)
                                                                                   caTATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGAcc                                                                   <  2:1568201/150‑1 (MQ=255)
                                                                                       ttGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGcc                                                               >  1:196094/1‑150 (MQ=255)
                                                                                        tGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCt                                                              >  1:1432675/1‑150 (MQ=255)
                                                                                         gAGGCAGGCATTCCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtg                                                             <  1:1175447/150‑1 (MQ=255)
                                                                                           ggcaggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtgat                                                           <  1:1317285/150‑1 (MQ=255)
                                                                                            gcaggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCtgatg                                                          <  1:988407/150‑1 (MQ=255)
                                                                                             caggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGc                                                                       >  2:688333/1‑136 (MQ=255)
                                                                                             caggcaTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGc                                                                       <  1:688333/136‑1 (MQ=255)
                                                                                                 caTTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCcg                                                     >  1:1126570/1‑150 (MQ=255)
                                                                                                     tCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGAc                                                 <  2:329170/150‑1 (MQ=255)
                                                                                                     tCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGAc                                                 <  2:501057/150‑1 (MQ=255)
                                                                                                       aCTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACaa                                               >  2:101841/1‑150 (MQ=255)
                                                                                                        cTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAAt                                              >  1:1579200/1‑150 (MQ=255)
                                                                                                         tgtgtgATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAAtt                                             <  2:288813/150‑1 (MQ=255)
                                                                                                           tgtgATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAc                                           <  2:1550718/150‑1 (MQ=255)
                                                                                                           tgtgATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTAc                                           <  2:1305727/150‑1 (MQ=255)
                                                                                                            gtgATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACa                                          >  2:1193343/1‑150 (MQ=255)
                                                                                                                tGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAAc                                      >  2:1377951/1‑150 (MQ=255)
                                                                                                                   ggTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGt                                   <  2:1190187/150‑1 (MQ=255)
                                                                                                                           ttCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACa                           <  2:1543019/150‑1 (MQ=255)
                                                                                                                             ctcGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGCCAATTACAGAACGGGAAATAACAtt                         >  1:301059/1‑150 (MQ=255)
                                                                                                                               cGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTAc                       <  2:1131129/150‑1 (MQ=255)
                                                                                                                               cGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCCCAACAATAACAGACCGTTAATAACAATAcc                      >  2:1083294/1‑150 (MQ=255)
                                                                                                                                 tgtaaaaaTACGAGATGGTGAATACAGAACAGACGTGATTCAGAGCCTGTAACCGCTCTCATTTAATAAACAGAATCGATATCCTCCTGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTAccc                     >  2:1273504/4‑150 (MQ=255)
                                                                                                                                 tGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTAcc                      <  2:1569094/149‑1 (MQ=255)
                                                                                                                                          cccggcTGGTCGGTACGCACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGt            >  1:330084/5‑150 (MQ=255)
                                                                                                                                          aCGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGt            >  1:680147/1‑150 (MQ=255)
                                                                                                                                           cgcgCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTa           <  1:12586/150‑1 (MQ=255)
                                                                                                                                              gcTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATcc        <  1:929143/150‑1 (MQ=255)
                                                                                                                                               cTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCg       <  1:446973/150‑1 (MQ=255)
                                                                                                                                                 ggTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGa     >  1:333206/1‑150 (MQ=255)
                                                                                                                                                 ggTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGa     >  2:162436/1‑150 (MQ=255)
                                                                                                                                                   tGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAgg   >  1:1474448/1‑150 (MQ=255)
                                                                                                                                                   ggCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAgg   <  2:1151171/149‑1 (MQ=255)
                                                                                                                                                    gCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGc  >  2:1408161/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TACTGCTGTCCCGACTGCTGGACCAGCACTATGGCCTCACACTGAATGACACACCTTTTGCCGATGAACGTGTGATTGAGCAGCATATTGAGGCAGGCATTTCACTGTGTGATGCGGTGAACTTTCTCGTGGAAAAATACGCGCTGGTGCGTACCGACCAGCCGGGATTCAGCGCCTGTACCCGCTCTCAGTTAATAAACAGCATCGATATCCTCCGGGCTCGCAGGGCGACCGGCCTGATGACCCGCGACAATTACAGAACGGTAAATAACATTACCCTGGGTAAGTATCCGGAGGC  >  NC_000913/2077639‑2077936

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: