Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,531,991 T→A 11.2% 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,9910TA11.2% 412.1 / 35.7 148intergenic (+31/‑137)PP_4018/topBacetyltransferase/DNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (75/54);  new base A (9/8);  total (86/62)
Fisher's exact test for biased strand distribution p-value = 7.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.73e-01

Only 100 of 152 total aligned reads displayed.

TCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGTTT  >  NC_002947/4531842‑4532139
                                                                                                                                                     |                                                                                                                                                    
tCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGtt                                                                                                                                                      >  1:1414312/1‑150 (MQ=255)
acAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGtt                                                                                                                                                      <  2:3813632/149‑1 (MQ=255)
     gcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAccc                                                                                                                                                 <  1:2752534/150‑1 (MQ=255)
     gcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAccc                                                                                                                                                 <  2:4033568/150‑1 (MQ=255)
     gcCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAAGTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTAccc                                                                                                                                                 <  2:2346758/150‑1 (MQ=255)
      cccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcg                                                                                                                                                <  2:3259042/150‑1 (MQ=255)
       ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTTCCcgc                                                                                                                                               >  1:3542196/1‑150 (MQ=255)
       ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                               <  2:2405735/150‑1 (MQ=255)
       ccAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCcgc                                                                                                                                               <  2:3181976/150‑1 (MQ=255)
        cAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCTTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCcgcg                                                                                                                                              <  2:572145/150‑1 (MQ=39)
        cAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCGCCTGAGCCGGCCTCTTCGCGGGTAAACCcgcg                                                                                                                                              <  2:2896655/150‑1 (MQ=39)
         aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                             >  2:2907510/1‑150 (MQ=255)
         aGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGa                                                                                                                                             <  2:1436286/150‑1 (MQ=255)
          ggTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTAAACCCGCGaa                                                                                                                                            <  1:3354844/150‑1 (MQ=255)
          ggTGCAAGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  2:852102/150‑1 (MQ=255)
           gTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAag                                                                                                                                           >  2:1763491/1‑150 (MQ=255)
            tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                          <  2:1261433/150‑1 (MQ=255)
            tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAaga                                                                                                                                          <  2:1356203/150‑1 (MQ=255)
            tGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACCGGCCTCTTCGCGGGTAAACCCGCGAaga                                                                                                                                          <  1:3564885/150‑1 (MQ=255)
             gCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAagag                                                                                                                                         <  1:3465906/150‑1 (MQ=255)
              cacgccacgcCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTGTACCCGCGAAGAgg                                                                                                                                        >  1:3922779/1‑150 (MQ=255)
               acgccacgccGTTCTATGGGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGc                                                                                                                                       >  1:268285/1‑150 (MQ=255)
                cgccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGcc                                                                                                                                      <  2:1588516/150‑1 (MQ=255)
                  ccacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTAACCGCGAAGAGGCCgg                                                                                                                                    >  1:4193189/1‑150 (MQ=255)
                   cacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGc                                                                                                                                   >  1:2934213/1‑150 (MQ=255)
                   cacgccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCCAAGAGGCCGGc                                                                                                                                   >  1:1809789/1‑150 (MQ=255)
                      gccGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCa                                                                                                                                >  2:2841924/1‑150 (MQ=255)
                      gccGTTCTAAGAGCGGCTGGGCTTTCGCGTAGTCAGAGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCa                                                                                                                                >  2:3168281/1‑150 (MQ=255)
                             tATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCCCGCAGAGGCCCGCGCAGGCGcta                                                                                                                         >  1:224751/1‑150 (MQ=255)
                             tATGAGCGGATGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAATATGGCATCCCGCATGTCGACATGGTGAGAGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGCGGGTAAACCCCAAAAGCCCGGCTCACGCGcta                                                                                                                         >  1:224710/1‑150 (MQ=255)
                               tGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAtata                                                                                                                       >  2:608684/1‑150 (MQ=255)
                               tGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAtata                                                                                                                       >  1:1404696/1‑150 (MQ=255)
                                      cTGGGCTTTCGCGTAGTCAGCGAGGAACTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATcccc                                                                                                                >  2:2788662/1‑150 (MQ=255)
                                        gggCTTTCGCGTTGTCAGCGAGGAATTCCTCGACGCTGGCATCCCGCTTGTCGCCCTGGGGCGCGCCTCGGGCGCGTGATCCATCACCTGCACTGGCCTCTTCGGGGGTTTCCCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGc                                                                                                              <  2:2511217/150‑1 (MQ=255)
                                            tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                          >  1:3255131/1‑150 (MQ=255)
                                            tttCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGAc                                                                                                          >  2:885330/1‑150 (MQ=255)
                                              tCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCCGGCGATATACATCCCCGCTGACaa                                                                                                        >  2:1692813/1‑150 (MQ=255)
                                                gcgTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGAGATACATCCCCGCTGACaaaa                                                                                                      >  2:2885122/1‑150 (MQ=255)
                                                       cAGCGTGTAATTCCTCGAAGGTGGCATCCCGAATGTCGACTTGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGCCCTCTTCGCGGTTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGt                                                                                               <  2:3542196/150‑1 (MQ=255)
                                                        aGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTa                                                                                              >  2:4119681/1‑150 (MQ=255)
                                                           gAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCATCGCGGGGTGACCCGCGAAGAGGCCGGCTCAGGCGATATACCTCCCCGCTGACAAAAAACCAGTAcat                                                                                           >  1:2253824/1‑150 (MQ=255)
                                                                agtCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTAcatccatc                                                                                      <  2:4121482/148‑1 (MQ=255)
                                                                 ttCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTTCCCGCGAAGAGGCCGGCGCAGGCGATATACATCCCCCCTGACAAAAAACCTGTAcagacaacc                                                                                     >  1:3834538/1‑143 (MQ=255)
                                                                 ttCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAATAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACaaccatcc                                                                                     >  1:1475395/1‑150 (MQ=255)
                                                                  tCCTCGACGGTGGAAACCCTCTTGTCGTCTGGGGGCGCGACTCGGGCGCTTGATCCATCACGTGCACTGGCCTCTTCGGGGTTTTCCCGGCTATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACAtccatcca                                                                                    <  1:2953672/150‑1 (MQ=255)
                                                                    cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGATGCCGGCGCAGGCGATATACAACCCCCCGGACAAAAAAACAGTAAAtacaaccaca                                                                                  >  1:2067081/1‑148 (MQ=255)
                                                                    cTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCGTACCGGGTTTTACCGCGAAATGGCCGGCGCAGGGGATATAAATACCCGCGGTCAAAAAACCTGTACAgccagccaca                                                                                  >  1:2263220/1‑148 (MQ=255)
                                                                     tCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTAGCGGGTTTACCCGCGAACAGGCCGGCTCAGGCGATATACATCCCCGCTGACCAAAAAACTGTACATCCAACCAGAt                                                                                 >  1:110692/1‑150 (MQ=255)
                                                                       gAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTGGCGGGTTTAACCGCGAAGAGGGCGGCTCAGGCGGTTTTAATACCCGCTTACAAAAAAACTGTAAATACATACACATaa                                                                               >  1:3038543/1‑150 (MQ=255)
                                                                         aggtGGAATCCCGCATGTCGACTTGGTGGGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATaaaa                                                                             <  2:3340419/147‑1 (MQ=255)
                                                                            tGGCACCCCGATTGTCGAATTGGGGCGCTACTCGGGCGCGTTATCAATAACCTGCACTGGCCTATTCGGGGTTTTACCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACtt                                                                          <  2:443916/150‑1 (MQ=255)
                                                                                ctCCCTCTTGTCTACTGGGTGCGGGAGGCGCGCGCGTTATCCTTCACCTGCAGGGCCCTCTTCGGGGTTTTACCGGCTAAGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCct                                                                      <  2:2003774/149‑1 (MQ=255)
                                                                                                gggggcgcGACTCGCGCGCCGGACCCATCCCCGGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                      <  2:3920686/146‑1 (MQ=255)
                                                                                                gggTGCGCGCCTCGCGCGCCTGATCAATCACCTGCAGTGGCCTGTTGGGGGGTTTACCGGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTgc                                                      <  1:2865122/149‑1 (MQ=255)
                                                                                                       cgACTCGCGCGCCTGATCCTTCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCt                                               <  2:320131/150‑1 (MQ=255)
                                                                                                       cgACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCt                                               <  2:2578432/150‑1 (MQ=255)
                                                                                                        gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACATGCCTCTGCGCCCCCGCTTGCGCATCCTa                                              >  2:1009742/1‑150 (MQ=255)
                                                                                                        gACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCCGGTTTACCCCCGAAGAGGCCGGCTCAGGCGATATAAAAACCCGCTGACAAAAAACCAGTAAATCCATACAGATAAAAAATTCCACTGCGCCCCCCCTTGCGCATccaa                                              >  1:2039560/1‑148 (MQ=255)
                                                                                                                cgcCTGATCCATCACCTGCACTGGCCTCTGCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCccac                                      <  2:4463944/150‑1 (MQ=255)
                                                                                                                    tGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCaccacc                                  >  1:1465497/1‑150 (MQ=255)
                                                                                                                     gctCCATCCCCTGCACGGGCCTCTTCGCGGGTTTACCGGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                                 <  2:3267471/148‑1 (MQ=255)
                                                                                                                     gATCCATCACCTGCACGGGCCTCTTGGGGGTTTTACCCGCTATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCAccc                                 <  2:176532/150‑1 (MQ=255)
                                                                                                                       tCCATCACCTGCACGGGCCTCTTGGGGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCg                               <  2:657578/150‑1 (MQ=255)
                                                                                                                        ccATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCAGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGa                              >  1:2164822/1‑150 (MQ=255)
                                                                                                                         cATCACCTGCACTGGCCTCTTGGCGGTTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                             <  2:2876599/150‑1 (MQ=255)
                                                                                                                         cATCACCTGCACTGGCCTCTTCGGGGGTTTACCGGCGAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGAt                             <  2:2447546/150‑1 (MQ=255)
                                                                                                                            cACCTGCAGTGGCCTGTTGGGGTTTTTACCCGCTAAGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  2:1993407/150‑1 (MQ=255)
                                                                                                                            cACCTGCAGTGGCCTCTTGGGGGTTTTCCCGGCTATGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  2:1901618/150‑1 (MQ=255)
                                                                                                                            cACCTGCACTGGCCTCTTCGGGGTTTTACCCGCTAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  2:4045453/150‑1 (MQ=255)
                                                                                                                            cACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGACGTGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  2:819955/150‑1 (MQ=255)
                                                                                                                            cACCTGAACTGGCCTCTGCGGGGTTTTACCCGCGAATAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  2:522966/150‑1 (MQ=255)
                                                                                                                                        ccTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGACGCGTCCCCGGcc              >  2:1860675/1‑150 (MQ=255)
                                                                                                                                        ccTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGcc              >  1:3182922/1‑150 (MQ=255)
                                                                                                                                         ctctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCa             >  2:3970111/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  1:3479544/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  1:3174268/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  2:1039496/1‑150 (MQ=255)
                                                                                                                                          tctTCGCGGGTAAACCCGCGAAGAGGCCAGTGCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCAt            >  2:3816380/1‑150 (MQ=39)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          <  2:4435873/150‑1 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:1118527/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  2:3438460/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  2:2912882/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:215217/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:787946/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:787995/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGcc                                         <  1:573012/119‑1 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  1:1738174/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgc          >  2:4479340/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGaa                          <  1:2815640/133‑1 (MQ=255)
                                                                                                                                              cgcgGGTAAACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATgcgc        >  2:2236108/1‑150 (MQ=255)
                                                                                                                                               gcgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGcc       >  1:2431176/1‑150 (MQ=255)
                                                                                                                                                ggggTTTCCCCGCGACGCGGCCGGCTCAGGCGATATACCTCCCCGCTGACAAAAAACCTGTACAGCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt      <  2:3495401/149‑1 (MQ=255)
                                                                                                                                                cgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt      >  2:1676957/1‑150 (MQ=255)
                                                                                                                                                cgGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCt      <  1:4119681/150‑1 (MQ=255)
                                                                                                                                                 gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTg     >  1:3069954/1‑150 (MQ=255)
                                                                                                                                                 gggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTg     >  1:3030189/1‑150 (MQ=255)
                                                                                                                                                 gggTTCCCCCGCGACGAGCCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTg     <  1:2418976/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAGCGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt    <  2:2651487/150‑1 (MQ=255)
                                                                                                                                                  ggTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGt    >  1:4481222/1‑150 (MQ=255)
                                                                                                                                                    tttCCCCGCCACGCGCCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGttt  <  1:960064/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGCGCCCAGGTGCACGCCACGCCGTTCTATGAGCGGCTGGGCTTTCGCGTAGTCAGCGAGGAATTCCTCGAAGCTGGCATCCCGCATGTCGACATGGTGCGCGACTCGCGCGCCTGATCCATCACCTGCACTGGCCTCTTCGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCGATATACATCCCCGCTGACAAAAAACCTGTACATCCATCCAGATAAAACTTGCCTCTGCGCCCCCGCTTGCGCATCCTATGGCCCACCACCCCGATGAAGCGTCCCCGGCCATGCGCCTGTTT  >  NC_002947/4531842‑4532139

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