Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 519,674 T→C 5.7% intergenic (+28/‑72) argC → / → erpA N‑acetyl‑gamma‑glutamyl‑phosphate reductase/iron‑sulfur cluster insertion protein ErpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947519,6740TC 5.7% 258.3 / 4.0 87intergenic (+28/‑72)argC/erpAN‑acetyl‑gamma‑glutamyl‑phosphate reductase/iron‑sulfur cluster insertion protein ErpA
Reads supporting (aligned to +/- strand):  ref base T (39/43);  new base C (3/2);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

GTGCTTTCGGTGATCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCACCCC  >  NC_002947/519527‑519816
                                                                                                                                                   |                                                                                                                                               
gtgCTTTCGGTGATCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGcgtc                                                                                                                                                >  2:108599/1‑149 (MQ=255)
        ggTGATCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGc                                                                                                                                       >  2:241363/1‑149 (MQ=255)
            aTCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCttt                                                                                                                                   >  7:151004/1‑149 (MQ=255)
              atccgatcTGTTCAAGGGTGCCTCCGTCCATGCGGTGCAGAAACTCAACATTCTGTTCGGTGTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGTCCGCCAGGATTGGAAAAGGCCCTCGTCGTGCGGGCCtt                                                                                                                                    <  5:213648/139‑1 (MQ=255)
               gACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTg                                                                                                                                <  1:241363/149‑1 (MQ=255)
               gACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTATGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGACAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTg                                                                                                                                >  8:167842/1‑149 (MQ=255)
                 cAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCACGACGCGGGCCttttccaa                                                                                                                              <  3:222509/149‑5 (MQ=255)
                  aaCCTGGTCAAGGGTGCCTCCGGCCAGTCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGCTGAGCGCATGGGCCTGTCCCATGGCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGTCCCGCGTCGTGCGGGCCTTTTTGCAg                                                                                                                             <  7:145855/149‑1 (MQ=255)
                  aaCCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAg                                                                                                                             >  7:34230/1‑149 (MQ=255)
                    ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTGGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCACGACGCGGGCCttttc                                                                                                                                 <  2:155654/143‑2 (MQ=255)
                    ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAgg                                                                                                                            >  2:221011/1‑148 (MQ=255)
                    ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGt                                                                                                                           <  2:249167/149‑1 (MQ=255)
                    ccTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCACGACGCGGGCCttttc                                                                                                                                 >  1:155654/1‑142 (MQ=255)
                     cTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTc                                                                                                                          >  2:292803/1‑149 (MQ=255)
                      tGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTcg                                                                                                                         >  1:93410/1‑149 (MQ=255)
                        gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCttt                                                                                                                                   >  1:10078/1‑137 (MQ=255)
                        gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCttt                                                                                                                                   <  2:10078/137‑1 (MQ=255)
                        gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTcgcg                                                                                                                       <  8:199651/149‑1 (MQ=255)
                        gTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTcgc                                                                                                                        <  1:162922/148‑1 (MQ=255)
                        gTCAAGGGTGCCTCCGGCCAGGAGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGACGCCAGAATTGGAAAAGTCCCGCGTCGTGCGGGCCtt                                                                                                                                    >  6:213648/1‑136 (MQ=255)
                              ggTGCCTCCGGCCAGGCGGTGCAGAACCTCAAGATTCTTTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTaaa                                                                                                                 <  4:171851/149‑1 (MQ=255)
                               gTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAg                                                                                                                <  8:34230/149‑1 (MQ=255)
                               gTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAg                                                                                                                >  2:336422/1‑149 (MQ=255)
                                   cTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGa                                                                                                                      <  3:5960/139‑1 (MQ=255)
                                   cTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGa                                                                                                                      >  4:5960/1‑139 (MQ=255)
                                    tCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGca                                                                                                           >  7:56190/1‑149 (MQ=255)
                                    tCCGGCCAGGCGGTGCAGAACCTCAACATTCTCTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGca                                                                                                           >  7:63808/1‑149 (MQ=255)
                                           aGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTc                                                                                                     <  1:221011/149‑1 (MQ=255)
                                             gCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCtt                                                                                                   >  2:244095/1‑148 (MQ=255)
                                                   cAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGAt                                                                                            >  6:90397/1‑149 (MQ=255)
                                                    aGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAgg                                                                                                                            <  7:118229/116‑1 (MQ=255)
                                                    aGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAgg                                                                                                                            >  8:118229/1‑116 (MQ=255)
                                                     gAACCTCAACATCCTGTTCGGGCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCGGCCGTAAGGCCGCCAGCATTCGAAACGGCCCGCGTCGTGCGGGCCTTTTTGCAGGGCGCGACTAAAGCCGCACAATTCTTGACCGAttt                                                                                          <  5:340762/149‑1 (MQ=255)
                                                         cTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTtctc                                                                                      >  7:69324/1‑149 (MQ=255)
                                                         cTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTtctc                                                                                      <  8:27784/149‑1 (MQ=255)
                                                            taCATTCTGTTCGGTCCGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                     <  4:111994/146‑1 (MQ=255)
                                                            aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGTCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTc                                                                                                     >  3:288695/1‑131 (MQ=255)
                                                            aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGTCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTc                                                                                                     <  4:288695/131‑1 (MQ=255)
                                                            aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCgg                                                                                    <  8:151004/148‑1 (MQ=255)
                                                            aaCATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                     <  1:292803/147‑1 (MQ=255)
                                                               aTTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCgg                                                                                    <  6:218215/145‑1 (MQ=255)
                                                               aTTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCgg                                                                                    >  5:218215/1‑145 (MQ=255)
                                                                  cTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                              <  5:90397/148‑1 (MQ=255)
                                                                  cTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                              <  8:56190/148‑1 (MQ=255)
                                                                  cTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                              <  1:108599/148‑1 (MQ=255)
                                                                     ttCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGa                                                                           >  3:46392/1‑148 (MQ=255)
                                                                         gTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                     <  4:14160/134‑1 (MQ=255)
                                                                         gTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCg                                                                                     >  3:14160/1‑134 (MQ=255)
                                                                         gTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATg                                                                      <  4:46392/149‑1 (MQ=255)
                                                                         gTCTGGATGAGCGCATGGGCCTGTCCCATGCCCGCTTGCTGCCATAAGGCCGCCAGGATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCTACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATg                                                                      >  6:166511/1‑149 (MQ=255)
                                                                          tCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATgc                                                                     <  6:58392/149‑1 (MQ=255)
                                                                              gATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                              >  6:290491/1‑136 (MQ=255)
                                                                              gATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGc                                                                              <  5:290491/136‑1 (MQ=255)
                                                                                    cgcATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGa                                                           <  4:186966/149‑1 (MQ=255)
                                                                                     gcATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATgcgc                                                                   >  1:96126/1‑140 (MQ=255)
                                                                                     gcATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATgcgc                                                                   <  2:96126/140‑1 (MQ=255)
                                                                                      cATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGCGt                                                         <  4:157513/149‑1 (MQ=255)
                                                                                       aTGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGtt                                                        <  2:93410/149‑1 (MQ=255)
                                                                                        tGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCATGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGttt                                                       <  1:244095/149‑1 (MQ=255)
                                                                                            gcTGTCCCATGCCGGCTTGCTGCCGTAAGGCATCCAGCATTGGAAAAGGCCGGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCTCAATTCTTGAACGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                    <  5:313965/147‑1 (MQ=255)
                                                                                              tGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATggc                                                 <  8:63808/149‑1 (MQ=255)
                                                                                                  ccATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATggcggc                                              >  3:255361/1‑148 (MQ=255)
                                                                                                   cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCtt                                                                                                   <  8:53608/94‑1 (MQ=255)
                                                                                                   cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCtt                                                                                                   >  7:53608/1‑94 (MQ=255)
                                                                                                   cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                    <  2:185380/141‑1 (MQ=255)
                                                                                                   cATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                    >  1:185380/1‑141 (MQ=255)
                                                                                                     tGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTAcc                                          >  3:170666/1‑149 (MQ=255)
                                                                                                       ccGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACcgc                                        >  1:83157/1‑149 (MQ=255)
                                                                                                        cGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                    >  6:247179/1‑136 (MQ=255)
                                                                                                        cGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTAt                                                    <  5:247179/136‑1 (MQ=255)
                                                                                                                   cGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAAt                             <  7:167842/148‑1 (MQ=255)
                                                                                                                     tAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCaa                          <  2:83157/149‑1 (MQ=255)
                                                                                                                     tAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCaa                          >  6:127938/1‑149 (MQ=255)
                                                                                                                         gccgccAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACAt                       >  8:262001/1‑148 (MQ=255)
                                                                                                                                  aTTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGaa             >  5:218282/1‑149 (MQ=255)
                                                                                                                                  aTTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGa              >  5:213392/1‑148 (MQ=255)
                                                                                                                                   tgcaaaaaGGCCCGCACGACGCGGGCCTTTTCGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAGTGCGCGTCATCGAGTTTTATGGCGGCTCCCGCGCCGGGAGAATCCGCATGAGCGTCgcgc            >  6:167741/5‑146 (MQ=255)
                                                                                                                                    tGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAAc           >  1:325168/1‑149 (MQ=255)
                                                                                                                                     ggAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAAcc          >  5:195904/1‑149 (MQ=255)
                                                                                                                                      gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCt                                             >  3:282827/1‑113 (MQ=255)
                                                                                                                                      gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCt                                             <  4:282827/113‑1 (MQ=255)
                                                                                                                                      gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAAcc          >  4:298844/1‑148 (MQ=255)
                                                                                                                                      gAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCt         <  7:78324/149‑1 (MQ=255)
                                                                                                                                      aaaaaGGCCCGCACGACGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCt         >  7:94558/2‑149 (MQ=255)
                                                                                                                                       aaaaGGGCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCtt        <  3:298844/149‑1 (MQ=255)
                                                                                                                                        aaaGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTc       <  6:195904/149‑1 (MQ=255)
                                                                                                                                            gCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCAccc   <  1:336422/149‑1 (MQ=255)
                                                                                                                                              ccGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCAcccc  >  5:198406/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GTGCTTTCGGTGATCGACAACCTGGTCAAGGGTGCCTCCGGCCAGGCGGTGCAGAACCTCAACATTCTGTTCGGTCTGGATGAGCGCATGGGCCTGTCCCATGCCGGCTTGCTGCCGTAAGGCCGCCAGCATTGGAAAAGGCCCGCGTCGTGCGGGCCTTTTTGCAGGTCGCGACTAAAGCCGCACAATTCTTGACCGATTTTCTCGGGAAAGCGGATAATGCGCGTCATCGAGTTTTATGGCGGCTACCGCGCCGGGAGAATCAACATGAGCGTCGAAACCTTCACCCC  >  NC_002947/519527‑519816

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