Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,969,020 A→C 6.4% I368I (ATA→ATC PP_3498 → transposase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,969,0200AC 6.4% 309.1 / 14.6 125I368I (ATA→ATCPP_3498transposase
Reads supporting (aligned to +/- strand):  ref base A (59/58);  new base C (5/3);  total (64/61)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.94e-01

Only 100 of 126 total aligned reads displayed.

AGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAGG  >  NC_002947/3968871‑3969160
                                                                                                                                                     |                                                                                                                                            
aGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATtata                                                                                                                                              <  1:552015/150‑1 (MQ=11)
   cTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGCCGACAGCTACCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATgg                                                                                                                                           <  2:594285/150‑1 (MQ=11)
    tgatgaTCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATggg                                                                                                                                          >  1:1797749/1‑150 (MQ=32)
    tgatgaTCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATggg                                                                                                                                          >  2:2820988/1‑150 (MQ=32)
     gatgatCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGc                                                                                                                                         >  2:2489785/1‑150 (MQ=32)
       tgatCCAGGGCGCGCATTCCGTCTTGANCCAGCTAAGGCCNGATTCCNATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGccc                                                                                                                                       >  1:2920076/1‑150 (MQ=17)
             aGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGNGATCAGNCTGGCTAATCGANATCTGCGGATTATATGGGCCCTGTTAc                                                                                                                                 <  1:2349792/150‑1 (MQ=17)
              gggCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACa                                                                                                                                >  1:1402432/1‑150 (MQ=32)
                    gcATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGa                                                                                                                          >  2:1744391/1‑150 (MQ=32)
                      aTTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATc                                                                                                                        >  1:1180225/1‑150 (MQ=32)
                        tCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAg                                                                                                                      <  1:293590/150‑1 (MQ=32)
                            tCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGttt                                                                                                                  >  2:1972939/1‑150 (MQ=32)
                                  gCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc                                                                                                            <  2:2800440/150‑1 (MQ=32)
                                  gCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAAc                                                                                                            <  1:1972939/150‑1 (MQ=32)
                                   ccAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg                                                                                                           >  2:1787834/1‑150 (MQ=32)
                                    cAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGc                                                                                                          >  1:2150303/1‑150 (MQ=32)
                                    cAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGATGGCTATACCGCCTTGGTCGCAGTGAGGCGGCGTTCAGTCTTGCTAGTCGATATCTGGGGATTATATAGGTCCTGGTACAGAGCGATCAGGTTTACCCACGc                                                                                                          >  2:2942922/1‑150 (MQ=11)
                                      gCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCaa                                                                                                        >  2:2359072/1‑150 (MQ=32)
                                        tAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAc                                                                                                      <  1:2489785/150‑1 (MQ=32)
                                         aaGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAAcc                                                                                                     >  2:1280424/1‑150 (MQ=32)
                                            gCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaa                                                                                                  <  2:3045501/150‑1 (MQ=32)
                                             ccTGATTCCGATCAGCCAGATGATCACCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                                                 >  2:1666850/1‑150 (MQ=17)
                                               tGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaacga                                                                                               >  1:2613184/1‑150 (MQ=32)
                                               tGATTCCGATCAGCCAGATGATCGCCACCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaacga                                                                                               >  1:2690687/1‑150 (MQ=17)
                                                 atttcgattagccagcctgatcgccgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcgAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGaaa                                                                                             >  2:2441401‑M2/86‑150 (MQ=255)
                                                      gATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACagag                                                                                                                             >  2:1936279/1‑113 (MQ=33)
                                                      gATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACagag                                                                                                                             <  1:1936279/113‑1 (MQ=32)
                                                      gATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                                                        >  1:1778294/1‑150 (MQ=32)
                                                        tCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                      <  1:172174/150‑1 (MQ=32)
                                                        tCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGAAAAAAAAAAGAAAAACCGGa                                                                                      <  2:2920076/150‑1 (MQ=11)
                                                          aGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCt                                                                                                                                      >  2:1088781/1‑100 (MQ=34)
                                                          aGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCt                                                                                                                                      <  1:1088781/100‑1 (MQ=33)
                                                             cAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAggcgg                                                                                 <  2:899576/150‑1 (MQ=32)
                                                                  gatCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAAAGAAAAACCGGAGGCGGCTATg                                                                            >  1:2941175/1‑150 (MQ=17)
                                                                      gcccccTCTTGCACTGGCTGTTCCGCCTTGGGCGCAAGTAGGCGCCGATCAGGCTGGCTAATAGAAATCTGCGGGTCATATCGTCCCTGTTACAGAGCGATCAGGATTACGAACGCAAACCGAACGAAACAACGGAGGCGGCTAtaatcc                                                                        >  2:2218199/1‑145 (MQ=11)
                                                                       ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGtt                                                                                                                                   >  1:2042248/1‑90 (MQ=35)
                                                                       ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGtt                                                                                                                                   <  2:2042248/90‑1 (MQ=35)
                                                                       ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg                                                                                                           >  2:894487/1‑114 (MQ=32)
                                                                       ccgccTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACg                                                                                                           <  1:894487/114‑1 (MQ=33)
                                                                       ccgccTCTTGCGATGGCTATCCCGCCATGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGACATGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAAAGAAAAACCGGAGGCGGCTATGAGCct                                                                       >  2:1403015/1‑150 (MQ=11)
                                                                         gcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaggctggctaatcgAAATCTGCGGATTATATGGGCCCTg                                                                                                                                     >  1:500986‑M2/62‑86 (MQ=255)
                                                                         cgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcgAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTa                                                                                                                 <  1:770337‑M2/45‑1 (MQ=255)
                                                                         cgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctaatcgAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCaa                                                                                                             <  1:677978‑M2/49‑1 (MQ=255)
                                                                         cgcctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctaatcgAAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCaa                                                                                                             >  2:677978‑M2/62‑110 (MQ=255)
                                                                          cctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaggctggctaatcgAAATCTGCGGATTATCTGGGCCCTGTCACAGAGCGATCAGGtt                                                                                                                   >  1:1731921‑M2/61‑103 (MQ=255)
                                                                          cctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcaggctggctaatcgAAATCTGCGGATTATCTGGGCCCTGTCACAGAGCGATCAGGtt                                                                                                                   <  2:1731921‑M2/43‑1 (MQ=255)
                                                                           ctctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                                                 <  1:393648/120‑1 (MQ=32)
                                                                           ctctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACcgaac                                                                                                 >  2:393648/1‑120 (MQ=32)
                                                                           ctctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAAGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                                   <  1:2112097/150‑1 (MQ=17)
                                                                           cctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcgGAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                                   >  1:2613510‑M2/60‑150 (MQ=255)
                                                                           cctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcgGAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                                   <  1:2465046‑M2/91‑1 (MQ=255)
                                                                            tctTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGAc                                                                  >  2:2775609/1‑150 (MQ=32)
                                                                            ctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcgGAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                      >  2:282920‑M2/59‑130 (MQ=255)
                                                                            ctccttgcgcccaaggcgggatagccagcgcaagaggcggcgatcatctggctgatcgGAATCTGCGGATTATCTGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                      <  1:282920‑M2/72‑1 (MQ=255)
                                                                              ttGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCa                                                                >  2:1255327/1‑150 (MQ=32)
                                                                                    tGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGtt                                                                                                                                   >  2:400967/1‑77 (MQ=35)
                                                                                      gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                                                        <  1:1225973/118‑1 (MQ=33)
                                                                                      gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCg                                                                                        >  2:1225973/1‑118 (MQ=32)
                                                                                      gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAg                                                                                     <  1:1527478/121‑1 (MQ=33)
                                                                                      gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAg                                                                                     >  2:1527478/1‑121 (MQ=32)
                                                                                      gCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAAAGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTa                                                        <  2:2941175/150‑1 (MQ=17)
                                                                                        tATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                                                      <  1:2775609/150‑1 (MQ=32)
                                                                                        tATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                                                      <  2:2397001/150‑1 (MQ=32)
                                                                                        tATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGAAGGCGGCTATCAGCCTCTGGCCAccctgttact                                                      >  2:613250/1‑143 (MQ=11)
                                                                                        tATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGATGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAt                                                      <  1:2624460/150‑1 (MQ=17)
                                                                                           cccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGcc                                                   >  1:1335327/1‑150 (MQ=32)
                                                                                            ccGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCa                                                  >  2:2934822/1‑150 (MQ=32)
                                                                                             cGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCAc                                                 >  2:2208307/1‑150 (MQ=32)
                                                                                               ccTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGa                                                                                      <  1:2025224/111‑1 (MQ=33)
                                                                                                cTTGGGCGCAAGGAGGCGGCGATCAGGATGGCTAGTCGAAATGTGCGGGTTATATGAGCCCTGTTACAGAGCGCTAAGGGTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCGAACgg                                              >  1:1651836/1‑150 (MQ=11)
                                                                                                  tGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAACCCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCtggg                                            <  2:1651836/150‑1 (MQ=17)
                                                                                                  tGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCgggg                                            >  2:2518901/1‑150 (MQ=32)
                                                                                                   gggCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCgg                                              <  2:1317683/147‑1 (MQ=32)
                                                                                                   gggCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGgc                                           <  2:2150303/150‑1 (MQ=32)
                                                                                                   gccagcgcaagaggcggcgatcatctggctgatcgAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGgc                                           <  1:2441401‑M2/115‑1 (MQ=255)
                                                                                                      cgcAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCt                                        <  1:2820988/150‑1 (MQ=32)
                                                                                                         aaGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAgg                                     <  2:1797749/150‑1 (MQ=32)
                                                                                                               gcggcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCTGCTATGAGCct                                                                       <  1:743301/110‑1 (MQ=17)
                                                                                                               gcggcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCct                                                                       >  2:743301/1‑110 (MQ=33)
                                                                                                               gcggcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                                   <  1:798769/114‑1 (MQ=32)
                                                                                                               gcggcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGa                                                                   >  2:798769/1‑114 (MQ=33)
                                                                                                               gcggcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCCACCGGGTCATGCCAGCGGGGCGATAGGCGCTcc                               >  2:1478109/1‑150 (MQ=255)
                                                                                                                  gcgATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGa                                                                           >  2:1459774/1‑103 (MQ=33)
                                                                                                                    gATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAAccc                          >  1:466530/1‑150 (MQ=32)
                                                                                                                        aGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGc                      >  2:1421819/1‑150 (MQ=32)
                                                                                                                         ggctggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTc                                <  2:2597170/139‑1 (MQ=33)
                                                                                                                         ggctggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTc                                >  1:2597170/1‑139 (MQ=33)
                                                                                                                            tggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGt                  >  1:790545/1‑150 (MQ=32)
                                                                                                                             ggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTc                                <  2:1418233/135‑1 (MQ=33)
                                                                                                                             ggctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTc                                >  1:1418233/1‑135 (MQ=33)
                                                                                                                              gctAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGTTAGGCGCGCCAACCCCTGCTAGtta                <  2:1402432/150‑3 (MQ=12)
                                                                                                                                tAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAg                                      >  2:1562782/1‑126 (MQ=33)
                                                                                                                                tAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAg                                      <  1:1562782/126‑1 (MQ=33)
                                                                                                                                 aaTCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACa             <  1:1421819/150‑1 (MQ=32)
                                                                                                                                     gAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTa                                                        >  1:988523‑M2/2‑103 (MQ=255)
                                                                                                                                      aaaTCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAc        <  1:1280424/150‑1 (MQ=32)
                                                                                                                                      aaaTCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGAc        <  1:2208307/150‑1 (MQ=32)
                                                                                                                                           tGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAg   <  2:2613184/150‑1 (MQ=32)
                                                                                                                                            gCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAgg  <  2:466530/150‑1 (MQ=32)
                                                                                                                                                aTTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAAGGCCACCGGGGCGCTAg                                      >  1:1573902/1‑110 (MQ=18)
                                                                                                                                                     |                                                                                                                                            
AGCCTGATGATCCAGGGCGCGCATTCCGTCTTGAGCCAGCTAAGGCCTGATTCCGATCAGCCAGATGATCGCCGCCTCTTGCGCTGGCTATCCCGCCTTGGGCGCAAGGAGGCGGCGATCAGGCTGGCTAATCGAAATCTGCGGATTATATGGGCCCTGTTACAGAGCGATCAGGTTTACCAACGCAAACCGAACGAAAAACCGGAGGCGGCTATGAGCCTCTGACCAACCGGGTAATGCCACCGGGGCGCTAGGCGCTCCAACCCCTGCTAGTGAACATTGACCCTAGG  >  NC_002947/3968871‑3969160

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: