Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,531,991 T→A 10.7% intergenic (+31/‑137) 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,9910TA10.7% 395.6 / 31.0 140intergenic (+31/‑137)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (77/48);  new base A (9/6);  total (86/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.92e-01

Only 100 of 144 total aligned reads displayed.

TCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4531842‑4532137
                                                                                                                                                     |                                                                                                                                                   
tCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGtt                                                                                                                                                     >  1:2630212/1‑150 (MQ=255)
tCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGtt                                                                                                                                                     >  2:5272433/1‑150 (MQ=255)
tCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGtt                                                                                                                                                     >  2:3039291/1‑150 (MQ=255)
  aGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTa                                                                                                                                                   >  1:1099018/1‑150 (MQ=255)
   gcgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAc                                                                                                                                                  <  2:3602378/150‑1 (MQ=255)
   gcgcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAc                                                                                                                                                  >  1:3859652/1‑150 (MQ=255)
     gcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAccc                                                                                                                                                <  1:5061509/150‑1 (MQ=255)
       ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                              >  1:4416399/1‑150 (MQ=255)
         aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                            <  2:4130364/150‑1 (MQ=255)
          ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                           <  1:75755/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)
            tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                         <  1:3556566/150‑1 (MQ=255)
            ggCACGCCACGCCGTTCTCTGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                         <  2:4993553/149‑1 (MQ=255)
               acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaacgccc                                                                                                                                      >  1:35709/1‑145 (MQ=255)
               acgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGGTTAACCGCGAAGAGGc                                                                                                                                      >  1:541001/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)
                                            tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                         >  2:4095015/1‑150 (MQ=255)
                                            tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCACCGCTGAc                                                                                                         >  2:899878/1‑150 (MQ=255)
                                            tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATACCCGCAGAc                                                                                                         >  1:3916562/1‑150 (MQ=255)
                                            tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAACCGCGAAGAGGCCCGCTAACGCGAAATACATCCCCCCAGAc                                                                                                         >  1:1405918/1‑150 (MQ=255)
                                                  gTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACaaaaaa                                                                                                   >  2:3160897/1‑150 (MQ=255)
                                                    aGTCAGCGAGTATTGCCTCGACGCTGGCAGCCCTCTTGTCTACTGGGTGCGCGACGCGCGCGCCTTAACCACCCCCTGCACTGGCCTCTGCGGGGTTTTCCCCGCTATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAcc                                                                                                 <  2:4121945/150‑1 (MQ=255)
                                                     gTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGTTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCt                                                                                                <  2:4986441/150‑1 (MQ=255)
                                                        aGCGAGGATTTCCTCGACGCTGGCACCCCGAATGCCGACTTGGTGCGCGACTCGCGCGCGTGATCAACCACGTGAACTGGCCTCTTCGGGGTTTTACCGGCTACGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                             <  2:2029104/150‑1 (MQ=255)
                                                         gCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCCCGGGTTTACCCGCGAAGAGGCCCGCTCAGGCGATATACATCCCCCCAGACCAAAAACCAGTAc                                                                                            >  1:2021153/1‑150 (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)
                                                                         tgCTGGCTTCCCGCATGTCTACTTGTTGCGCTACTCGCGCGCCTGATCAATCACCTGCACTGGCCTCTTCGCGGTTTTACCCGCTTATAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaaa                                                                            <  2:4514810/149‑1 (MQ=255)
                                                                            tGGAATCCCGATTTTCGAATTGTTGCGCTACTCGCGCGCCTGTTCAATCACCTGAACTGGCCTCTTCGCGGTTTTACCCGCGATGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACtt                                                                         <  2:1643655/150‑1 (MQ=255)
                                                                              gaaTCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCAGTGGCCTCTACGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGc                                                                       <  2:5222848/148‑1 (MQ=255)
                                                                              gCACCCCGCATGGCGACTGGGGGCGCGACTCGCGCGCCTGATCCAGCACCTGCACTGGCCTCTCCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGc                                                                       <  2:4404494/150‑1 (MQ=255)
                                                                                  cccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCCCGAAGAGGCCGGCTCAGGCGATAAACAACCCCGCTGGCAAAAAAACAGTAAATACATACAGATAAAAATTGCcact                                                                   >  1:1118291/1‑150 (MQ=255)
                                                                                   ccGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTg                                                                  >  2:4098090/1‑150 (MQ=255)
                                                                                    cGCTTGTCGACAGGGGGCGCTACGCGCGCGCCTTAGCCATCACGTGCCCGGGCCTCTTGGGGGGTTTCCCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTgc                                                                 <  2:2526494/150‑1 (MQ=255)
                                                                                            gACTTGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGc                                                        <  2:3542125/150‑1 (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)
                                                                                                tGGGGCGCGACTCGGGCGCCTTTTCCAACACCTGCAGTGGCCTCTGCGGGGTTTTCCCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    <  1:2191737/150‑1 (MQ=255)
                                                                                                gggTGCGCTACCCGCGCGCCTGATCCATCACCGGCAGGGGCCTCTTCGGGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    <  2:3764323/149‑1 (MQ=255)
                                                                                                   ggcgcgcCTCGCGCCCGTTAACAAACACCTGCACGGGCCTCTGGGGGGTTTTCCCCGCGACGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCa                                                 <  2:1000521/149‑1 (MQ=255)
                                                                                                     cgcgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCGTCGCGGGTTTTCCCCCGAAGAGGCCGGCGCAGGCGATATACATCCCCCCGGACAAAAAAACAGTAAAACCAACCAGATAAAAATTTCCTCAGCCCCCCCGCTTGCGCaac                                               >  1:1499681/1‑148 (MQ=255)
                                                                                                        gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGAGAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCtt                                            >  2:4085728/1‑149 (MQ=255)
                                                                                                           gcgcgcgCCTTCCCCAGCACCGGCCCGGGCCTCTTCGGGGGTTTCCCCGCGACGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATgg                                         <  2:2630212/149‑1 (MQ=255)
                                                                                                               gcgcCTTACCCACCCCCTTCAGTGGCCTCTGGGGGGTTTTACCCGCTAAGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCcca                                     <  1:2125885/150‑1 (MQ=255)
                                                                                                                    tGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAAACTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCGTGCGCATCCTATGGCCCaccacc                                >  1:1246702/1‑150 (MQ=255)
                                                                                                                        ccATCCCGTGCACGGGCCTCTTGGGGGGTTTACTGGCGATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGa                            <  1:3109762/150‑1 (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)
                                                                                                                          aTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCCGCTCAGGCGAGATACATCCCCCCTGACAAAAAACCAGTACATCCATCCAGATAAAACATTCCTCAGCGCCCCCGCTTGCGCATCCAATGGCCCACCACCCCCATg                          >  1:4323114/1‑150 (MQ=255)
                                                                                                                           tCACCTGCAGTGGCCTCTTCGGGGTTTTCCCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGa                         <  1:4459837/150‑1 (MQ=255)
                                                                                                                               cTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCg                     >  2:4490121/1‑150 (MQ=255)
                                                                                                                               cTGCACTGGCCTCTTAGCGGGTTTAACCGCGAAGAGGCCGGCTCAGGCGATATAAATACCCCCAGACAAAAAACCTGTAAAAACATCCAGATAAAAATTGCCTCTGCGCCCCCCCGTGCGCATCCAATGGCCCACCACCCCGATGAAGCg                     >  1:4404500/1‑150 (MQ=255)
                                                                                                                                   cgtGGCCTCTCCGGGGTTTTACCGGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTccc                 <  1:2801769/148‑1 (MQ=255)
                                                                                                                                     tGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCg               >  2:2372854/1‑150 (MQ=255)
                                                                                                                                        ccTCTTCGCGGTTTTACCCGCGACGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGcc            <  2:2380710/150‑1 (MQ=255)
                                                                                                                                         ctctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa           >  2:3603991/1‑150 (MQ=255)
                                                                                                                                         ctctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa           >  2:5097902/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  2:1903915/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  1:2550996/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  1:4249342/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  2:3438040/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTTTACCCGCGAAGAGGCCAGTGCAGGTGATGTACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  1:616717/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCcca                                     <  2:1448150/123‑1 (MQ=255)
                                                                                                                                          tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  1:3969511/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTAAACCCGCGAAGAGGCCAGTGCAGGTGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt          >  2:4410561/1‑150 (MQ=39)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCCGATACAACTTGCCTCTGCGCCCCCGCTTGCGCAGCCTATGGCCCCCCACCCCGATGAAGCGTCCCCGGCCATgc        >  2:4950384/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCCTgc        >  2:4950017/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc        >  2:1189741/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc        >  2:290456/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcg       >  2:5430727/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    >  1:3706798/1‑105 (MQ=255)
                                                                                                                                             tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                    <  2:3706798/105‑1 (MQ=255)
                                                                                                                                             tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcg       >  1:589792/1‑150 (MQ=255)
                                                                                                                                              cgcgGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc      <  2:589792/150‑1 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc      >  1:157783/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  1:37934/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  2:1006301/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  2:3588180/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  1:1514305/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc     >  1:3066280/1‑150 (MQ=255)
                                                                                                                                                 gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAAc                                                                           >  2:1144632/1‑79 (MQ=255)
                                                                                                                                                 gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTc                   <  1:2291509/134‑1 (MQ=255)
                                                                                                                                                 gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTc                   >  2:2291509/1‑134 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt  <  2:2066594/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGT  >  NC_002947/4531842‑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: