Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,311,745 T→C 6.7% intergenic (‑207/+71) PP_1144 ← / ← rapA membrane protein/RNA polymerase‑associated protein RapA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,311,7450TC 6.7% 224.5 / 4.4 75intergenic (‑207/+71)PP_1144/rapAmembrane protein/RNA polymerase‑associated protein RapA
Reads supporting (aligned to +/- strand):  ref base T (33/37);  new base C (3/2);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-01

GCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGC  >  NC_002947/1311601‑1311894
                                                                                                                                                |                                                                                                                                                     
gCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGt                                                                                                                                                  <  2:87149/150‑1 (MQ=255)
       aggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTgg                                                                                                                                                >  2:320157/1‑145 (MQ=255)
       aggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTgg                                                                                                                                                <  1:320157/145‑1 (MQ=255)
        ggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCa                                                                                                                                          >  2:374763/1‑150 (MQ=255)
        ggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCa                                                                                                                                          >  1:783610/1‑150 (MQ=255)
        ggaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGCTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGccc                                                                                                                                          >  2:5509/1‑149 (MQ=255)
         gaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCAc                                                                                                                                         >  2:971411/1‑150 (MQ=255)
         gaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCAc                                                                                                                                         <  1:1324310/150‑1 (MQ=255)
         gaaggaGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCAc                                                                                                                                         >  1:782918/1‑150 (MQ=255)
               gTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCtt                                                                                                                                   <  1:971411/150‑1 (MQ=255)
                   ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACCGCAGTGGATGc                                                                                                                                            <  1:1035599/137‑1 (MQ=255)
                   ttAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACCGCAGTGGATGc                                                                                                                                            >  2:1035599/1‑137 (MQ=255)
                        ccTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgtgt                                                                                                                          >  1:1466485/1‑150 (MQ=255)
                           tCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTcg                                                                                                                       <  2:1226082/150‑1 (MQ=255)
                              tGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGATCAGGGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                    <  1:1277472/150‑1 (MQ=39)
                                    aTGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGc                                                                                                              <  2:645148/150‑1 (MQ=255)
                                      ggCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGccc                                                                                                            <  1:953170/150‑1 (MQ=255)
                                      ggCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGccc                                                                                                            <  1:1551697/150‑1 (MQ=255)
                                         acaATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCt                                                                                                         >  1:305874/1‑150 (MQ=255)
                                          caATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCa                                                                                                                                          >  1:75558/1‑116 (MQ=255)
                                          caATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCa                                                                                                                                          <  2:75558/116‑1 (MQ=255)
                                            aTGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGt                                                                                                                                                  <  1:1421758/106‑1 (MQ=255)
                                            aTGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGt                                                                                                                                                  >  2:1421758/1‑106 (MQ=255)
                                                 ttCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACaat                                                                                                 <  1:975008/150‑2 (MQ=255)
                                                  tCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGACGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAgg                                                                                                <  2:782918/150‑1 (MQ=255)
                                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                    <  1:1386032/128‑1 (MQ=255)
                                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                    >  2:1386032/1‑128 (MQ=255)
                                                    tACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAgggg                                                                                              >  2:30438/1‑150 (MQ=255)
                                                       gATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCTCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCCTCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGgcgc                                                                                           <  2:848853/150‑1 (MQ=255)
                                                              ttATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGcc                                                                                                             <  1:1210736/125‑1 (MQ=255)
                                                              ttATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGcc                                                                                                             >  2:1210736/1‑125 (MQ=255)
                                                               tataACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGc                                                                                                              >  1:648218/1‑123 (MQ=255)
                                                               tataACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGc                                                                                                              <  2:648218/123‑1 (MQ=255)
                                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                    >  2:1211138/1‑111 (MQ=255)
                                                                     tCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCgg                                                                                                                    <  1:1211138/111‑1 (MQ=255)
                                                                      ccAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagc                                                                            <  2:897625/150‑1 (MQ=255)
                                                                        aaTCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCggg                                                                                        >  1:1363392/1‑136 (MQ=255)
                                                                        aaTCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCggg                                                                                        <  2:1363392/136‑1 (MQ=255)
                                                                            gCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAgccagcca                                                                      >  2:982986/1‑150 (MQ=255)
                                                                                  aGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCa                                                                >  2:769184/1‑150 (MQ=255)
                                                                                    gAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCAAGCCACCAGCAcg                                                              >  1:242636/1‑150 (MQ=255)
                                                                                     aaTAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCAcgc                                                             <  2:181004/150‑1 (MQ=255)
                                                                                      aTAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgt                                                                                                                            <  2:1118066/86‑1 (MQ=255)
                                                                                      aTAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGgt                                                                                                                            >  1:1118066/1‑86 (MQ=255)
                                                                                       tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTa                                                                              <  2:677142/131‑1 (MQ=255)
                                                                                       tAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTa                                                                              >  1:677142/1‑131 (MQ=255)
                                                                                         aaTTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAAt                                                         <  2:783610/150‑1 (MQ=255)
                                                                                              ggggCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCtt                                                    >  2:1597839/1‑150 (MQ=255)
                                                                                               gggCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTc                                                   <  1:374763/150‑1 (MQ=255)
                                                                                                  cccTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAg                                                <  1:5509/150‑1 (MQ=255)
                                                                                                   ccTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagccagcc                                                                       >  1:361882/1‑126 (MQ=255)
                                                                                                   ccTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTagccagcc                                                                       <  2:361882/126‑1 (MQ=255)
                                                                                                          ccTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAggg                                        >  1:686859/1‑150 (MQ=255)
                                                                                                            ttCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcg                                      >  1:1617758/1‑150 (MQ=255)
                                                                                                              ccTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcggc                                    <  1:769184/150‑1 (MQ=255)
                                                                                                                 gtgGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCtt                                 <  1:538901/150‑1 (MQ=255)
                                                                                                                  tgGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGc                                                                 >  2:1428924/1‑117 (MQ=255)
                                                                                                                  tgGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGc                                                                 <  1:1428924/117‑1 (MQ=255)
                                                                                                                  tgGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTtc                                >  2:670247/1‑150 (MQ=255)
                                                                                                                         gggCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCa                         <  2:1617758/150‑1 (MQ=255)
                                                                                                                              aGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGcc                    >  1:395961/1‑150 (MQ=255)
                                                                                                                                cccGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAg                  >  2:992683/1‑150 (MQ=255)
                                                                                                                                  cgcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGc                     <  2:1034545/145‑1 (MQ=255)
                                                                                                                                  cgcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGc                     >  1:1034545/1‑145 (MQ=255)
                                                                                                                                   gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcgg                                     <  1:1212417/128‑1 (MQ=255)
                                                                                                                                   gcgAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGggcgg                                     >  2:1212417/1‑128 (MQ=255)
                                                                                                                                        gCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCtt          <  2:305874/150‑1 (MQ=255)
                                                                                                                                        gCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTc                                                   >  2:137781/1‑109 (MQ=255)
                                                                                                                                        gCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTc                                                   <  1:137781/109‑1 (MQ=255)
                                                                                                                                        gCATCCACCGCAGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCtt          >  2:1204454/1‑150 (MQ=255)
                                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCggg                                                                                        >  2:557907/1‑71 (MQ=255)
                                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCggg                                                                                        <  1:557907/71‑1 (MQ=255)
                                                                                                                                         cATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTc         <  1:1597839/150‑1 (MQ=255)
                                                                                                                                           tCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTcgc       <  2:1428849/150‑1 (MQ=255)
                                                                                                                                              aCTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGcgcc    >  1:527159/1‑150 (MQ=255)
                                                                                                                                                tGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGcgccgc  >  2:1629081/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GCAGTATAGGAAGGAGTACTTAAGCCTTCCTGACTAATGGCACAATGATTTCTACGATTTGGTTATAACTCCAATCGCAGTGAGGAATAAATTTGGGGCCCTAATGCCTTCCTGTGGGAGCGGGCGAGCCCGCGAAGCATCCACTGCGGTGGATGCCACCGGCTTTGCCGGTGTTCGCGGGCGTGCCCGCTCCCACAAGGGGCGCGGGGCCTGCGCTAGCCAGCCACCAGCACGCGAATGGCTTCCAGGCGCAGGGCGGCCTTCTCCAGCATTGCCAGGCCATCTTCGCGCCGC  >  NC_002947/1311601‑1311894

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