Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,531,992 T→A 12.0% intergenic (+32/‑136) PP_4018 → / → topB acetyltransferase/DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,531,9920TA12.0% 367.8 / 30.7 140intergenic (+32/‑136)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (77/45);  new base A (9/7);  total (86/54)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 140 total aligned reads displayed.

AGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4531844‑4532137
                                                                                                                                                    |                                                                                                                                                  
aGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTa                                                                                                                                                   >  1:1099018/1‑150 (MQ=255)
 gcgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAc                                                                                                                                                  >  1:3859652/1‑150 (MQ=255)
 gcgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAc                                                                                                                                                  <  2:3602378/150‑1 (MQ=255)
   gcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAccc                                                                                                                                                <  1:5061509/150‑1 (MQ=255)
    cccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcg                                                                                                                                               >  1:3112795/1‑150 (MQ=255)
     ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                              >  1:4416399/1‑150 (MQ=255)
     ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                              >  2:5029601/1‑150 (MQ=255)
     ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTACCcgc                                                                                                                                              >  1:2321213/1‑150 (MQ=255)
       aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                            <  2:4130364/150‑1 (MQ=255)
        ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                           <  1:75755/150‑1 (MQ=255)
        ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                           <  1:2147538/150‑1 (MQ=255)
        ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAATCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                           <  2:3055997/150‑1 (MQ=255)
        ggTGCACGCCACGCCGTTCTAGGAGCGGCTGGGCTTTCGCGGAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGGGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                           <  2:3757838/150‑1 (MQ=255)
         gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTTCCCGCGAag                                                                                                                                          >  1:1850389/1‑150 (MQ=255)
          tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                         <  1:2826418/150‑1 (MQ=255)
          tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                         <  1:3556566/150‑1 (MQ=255)
          ggCACGCCACGCCGTTCTCTGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                         <  2:4993553/149‑1 (MQ=255)
             acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaacgccc                                                                                                                                      >  1:35709/1‑145 (MQ=255)
             acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGc                                                                                                                                      >  1:541138/1‑150 (MQ=255)
               gccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                                    >  2:4685326/1‑150 (MQ=255)
                    gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGACGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGACGAGGCCGGCTCa                                                                                                                               >  2:1694339/1‑150 (MQ=255)
                           taggcgcggcgGGGCTTGCGCTTAGTCTGCGAGGATTTCCGCTAAGCTGGAACCCCGCTTGTCGACTTGGTGCGCGACTCGCGCGCCTTACCCATCCCCGGCCCTGGCCTATCCGCGGTTTTACCGGCTAATAGGCCGGCTCAGGCGata                                                                                                                        <  2:3674718/145‑1 (MQ=255)
                           tATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCCCGAAGAGGCCGGCTCAGGCGata                                                                                                                        >  1:1669150/1‑150 (MQ=255)
                              gAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATAc                                                                                                                     >  1:3240720/1‑150 (MQ=255)
                                     tGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATccccc                                                                                                              >  1:4337442/1‑149 (MQ=255)
                                        gCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTg                                                                                                           >  2:552566/1‑150 (MQ=255)
                                         ttttCGCGTTGTCAGCGAGTATTTCCTCGAAGCTGGCAGCCCGCATGTCGACTTGGTGCGCGACTCGCGCGCCTGATCCATCACGTGCACTGGCCTCTTGGCGGGTTTACCGGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGa                                                                                                          <  2:3627333/149‑1 (MQ=255)
                                          tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                         >  2:4095015/1‑150 (MQ=255)
                                          tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCACCGCTGAc                                                                                                         >  2:899878/1‑150 (MQ=255)
                                          tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATACCCGCTGAc                                                                                                         >  2:1009949/1‑150 (MQ=255)
                                          tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATACCCGCAGAc                                                                                                         >  1:3916562/1‑150 (MQ=255)
                                          tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAACCGCGAAGAGGCCCGCTAACGCGAAATACATCCCCCCAGAc                                                                                                         >  1:1405918/1‑150 (MQ=255)
                                                  aGTCAGCGAGTATTGCCTCGACGCTGGCAGCCCTCTTGTCTACTGGGTGCGCGACGCGCGCGCCTTAACCACCCCCTGCACTGGCCTCTGCGGGGTTTTCCCCGCTATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAcc                                                                                                 <  2:4121945/150‑1 (MQ=255)
                                                      aGCGAGGATTTCCTCGACGCTGGCACCCCGAATGCCGACTTGGTGCGCGACTCGCGCGCGTGATCAACCACGTGAACTGGCCTCTTCGGGGTTTTACCGGCTACGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                             <  2:2029104/150‑1 (MQ=255)
                                                         gtgcatttCCTCGATGGGGGCCCCCCGCTTGTCGACCTGGGGCGCCCCGCGCGCGCCTTATCCATCCCGGGCTGTGGCCTGTTCGGGGTTTAACCGGCTACGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcat                                                                                          <  2:541138/144‑1 (MQ=255)
                                                          aGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatc                                                                                         >  1:5381596/1‑150 (MQ=255)
                                                             attcccTCGACGCGGGAACCCCTCTTGCCTACAGGGGGCGCCACTCGCGCGCCTGAACAATCCCCTGCCCGGCCCTCTTCGCGGGTTTACCCGCTATGAGCCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccat                                                                                      <  2:4595537/146‑1 (MQ=255)
                                                             aaTCCCGCGAAGGGGCCATCCCGCATTTCGACATGGTGCGGGACTCGGGCGCCTGATCCAGCCCCTGCACTGCCCTCTCCGCGGGTTCACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccat                                                                                      <  2:4595531/150‑1 (MQ=255)
                                                              aTTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACAAACCCCCTGACAAAAAACCAGTAcaaccaac                                                                                     >  1:1421661/1‑148 (MQ=255)
                                                                  cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGa                                                                                 <  2:4416399/150‑1 (MQ=255)
                                                                  cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGCTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGa                                                                                 >  2:1208685/1‑150 (MQ=255)
                                                                      aaGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTTATCAATAACCGGAACTGGCCTCTTCGCGTGTTTACCCGCGAATAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaa                                                                             <  2:3548453/150‑1 (MQ=255)
                                                                          tGGAATCCCGATTTTCGAATTGTTGCGCTACTCGCGCGCCTGTTCAATCACCTGAACTGGCCTCTTCGCGGTTTTACCCGCGATGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACtt                                                                         <  2:1643655/150‑1 (MQ=255)
                                                                            gCACCCCGCATGGCGACTGGGGGCGCGACTCGCGCGCCTGATCCAGCACCTGCACTGGCCTCTCCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGc                                                                       <  2:4404494/150‑1 (MQ=255)
                                                                                cccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCCCGAAGAGGCCGGCTCAGGCGATAAACAACCCCGCTGGCAAAAAAACAGTAAATACATACAGATAAAAATTGCcact                                                                   >  1:1118291/1‑150 (MQ=255)
                                                                                 ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTg                                                                  >  2:4098090/1‑150 (MQ=255)
                                                                                 ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTg                                                                  >  2:3720838/1‑150 (MQ=255)
                                                                                          gACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCGATATAAATCCCCGCTGTCAAAAAAAAAGTACATACATACAGATAAAAATTTCCTCTGCGGCCCCGc                                                        >  1:3277797/1‑150 (MQ=38)
                                                                                          gACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTTCCCGCGAAGAGGCCCGCGCAGGCGCTATACATCCCCGCGGACAAAAAAACAGTAAATACATCCAGATAAAACTTGCCGCTGCGccccccc                                                        >  1:522853/1‑148 (MQ=255)
                                                                                          gACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCCGGTTTTCCCCCGCAGAGGCCGGCTCAGGCGATATAAAACCCCCCTGACAAAAAAACAGTACAACCAAACAGGAAAAACGTTGCTCTGCGCCCCCGc                                                        >  1:1000521/1‑150 (MQ=255)
                                                                                           aCATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCt                                                       >  2:731705/1‑150 (MQ=255)
                                                                                           aCATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAACCCCCAACAGGCCGGCTCAGGCGATAAAAATCCCCCCGGACAAAAAAACTGTACAAACATCCAGATAAAACATGCCACAGCGCCCCCCCt                                                       >  1:1198635/1‑150 (MQ=255)
                                                                                            cagggggcgcGACTCGCGCGCCGTAACCAGCACCGGCACGGGCCTCGTCGCGGGTTTCCCCGCGAAGAGGCCGGCTCAGGCGATATACAGCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCtt                                                      <  2:3052258/144‑1 (MQ=255)
                                                                                              gggTGCGCTACCCGCGCGCCTGATCCATCACCGGCAGGGGCCTCTTCGGGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    <  2:3764323/149‑1 (MQ=255)
                                                                                                 ggcgcgcCTCGCGCCCGTTAACAAACACCTGCACGGGCCTCTGGGGGGTTTTCCCCGCGACGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCa                                                 <  2:1000521/149‑1 (MQ=255)
                                                                                                      gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGAGAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCtt                                            >  2:4085728/1‑149 (MQ=255)
                                                                                                       tcTCGCGCGCCTGATCCATCACCTGGACTGGCCTCTGCGCGGGTTTTCCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATTAAACGTGCCTCTGCGCCCCCTCTTGCGCATCCTAt                                           <  2:4960739/149‑1 (MQ=255)
                                                                                                         tCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATgg                                         <  2:5355960/150‑1 (MQ=255)
                                                                                                         gcgcgcgCCTTCCCCAGCACCGGCCCGGGCCTCTTCGGGGGTTTCCCCGCGACGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATgg                                         <  2:2630212/149‑1 (MQ=255)
                                                                                                             gcgcCTTACCCACCCCCTTCAGTGGCCTCTGGGGGGTTTTACCCGCTAAGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCcca                                     <  1:2125885/150‑1 (MQ=255)
                                                                                                                  tGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAACTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCGTGCGCATCCTATGGCCCaccacc                                >  1:1246702/1‑150 (MQ=255)
                                                                                                                      ccATCACCGGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGa                            <  2:2711586/150‑1 (MQ=255)
                                                                                                                       cATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGTAGAGGCCGGCTCAGGCGATATACATACCCGCTGACAAAAAAACTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                           >  2:1320455/1‑150 (MQ=255)
                                                                                                                       cATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCCGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCCCTTGCGCATCCTATGGCCCACCACCCCGAt                           >  2:4697409/1‑150 (MQ=255)
                                                                                                                       cATCACCTGCACTGGCCTCTGCGCGGGTTTACCCGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                           <  2:55913/150‑1 (MQ=255)
                                                                                                                        aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCGTGCGCATCCTATGGCCCACCACCCCGATg                          >  2:3175688/1‑150 (MQ=255)
                                                                                                                        aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCCGCTCAGGCGAGATACATCCCCCCTGACAAAAAACCAGTACATCCATCCAGATAAAACATTCCTCAGCGCCCCCGCTTGCGCATCCAATGGCCCACCACCCCCATg                          >  1:4323114/1‑150 (MQ=255)
                                                                                                                         tCACCTGCAGTGGCCTCTTCGGGGTTTTCCCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGa                         <  1:4459837/150‑1 (MQ=255)
                                                                                                                             cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCg                     >  2:4490121/1‑150 (MQ=255)
                                                                                                                                   tGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCg               >  2:2372854/1‑150 (MQ=255)
                                                                                                                                      ccTCTTCGCGGTTTTACCCGCGACGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGcc            <  2:2380710/150‑1 (MQ=255)
                                                                                                                                       ctctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa           >  2:5097902/1‑150 (MQ=255)
                                                                                                                                       ctctTCGCGGGTAAACCCGCGAAGAGGCCAGTGCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa           >  1:640187/1‑150 (MQ=39)
                                                                                                                                        tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  2:3438040/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  1:2550996/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  2:1903915/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  1:4249342/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCcca                                     <  2:1448150/123‑1 (MQ=255)
                                                                                                                                        tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCcca                                     >  1:1448150/1‑123 (MQ=255)
                                                                                                                                        tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  1:3969511/1‑150 (MQ=255)
                                                                                                                                        tctTCGCGGGTAAACCCGCGAAGAGGCCAGTGCAGGTGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  2:4410561/1‑150 (MQ=39)
                                                                                                                                          ttcactggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc        >  1:1956168/7‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCCGATACAACTTGCCTCTGCGCCCCCGCTTGCGCAGCCTATGGCCCCCCACCCCGATGAAGCGTCCCCGGCCATgc        >  2:4950384/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc        >  2:1189883/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc        >  2:290456/1‑150 (MQ=255)
                                                                                                                                           tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    >  1:3706798/1‑105 (MQ=255)
                                                                                                                                           tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcg       >  1:589792/1‑150 (MQ=255)
                                                                                                                                            cgcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc      >  2:2438419/1‑150 (MQ=255)
                                                                                                                                            cgcgGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc      <  2:589792/150‑1 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc      >  1:157783/1‑150 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAcccc                              >  2:4493525/1‑125 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  1:37934/1‑150 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  1:1514305/1‑150 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  2:3588180/1‑150 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  1:3066280/1‑150 (MQ=255)
                                                                                                                                              cgGTTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt    <  2:4912902/150‑1 (MQ=255)
                                                                                                                                               gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAAc                                                                           >  2:1144632/1‑79 (MQ=255)
                                                                                                                                               gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTc                   >  2:2291509/1‑134 (MQ=255)
                                                                                                                                               gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTc                   <  1:2291509/134‑1 (MQ=255)
                                                                                                                                                ggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt  <  2:2066594/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4531844‑4532137

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