Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,433:1 +T 7.7% intergenic (‑51/+38) trmE ← / ← oxaA GTPase/membrane protein insertase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,4331.T 7.7% 348.6 / 21.0 91intergenic (‑51/+38)trmE/oxaAGTPase/membrane protein insertase
Reads supporting (aligned to +/- strand):  ref base . (35/49);  new base T (5/2);  total (40/51)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

CTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGG  >  NC_002947/6287‑6574
                                                                                                                                                      |                                                                                                                                              
cTTGGCAGCCAGAGGGCCGGACAACCTGACATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAgggg                                                                                                                                             >  8:266890/1‑149 (MQ=255)
   ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGc                                                                                                                                           >  5:178713/1‑149 (MQ=255)
   ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAAG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACGAgggggg                                                                                                                                           >  8:348659/1‑148 (MQ=255)
      aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtt                                                                                                                                       >  8:226157/1‑149 (MQ=255)
        ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                                                                                                     <  7:209605/149‑1 (MQ=255)
        ccAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCCTAGTGGGGCGtttt                                                                                                                                     <  6:316892/149‑1 (MQ=255)
         cAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGttttt                                                                                                                                    <  5:125251/149‑1 (MQ=255)
           gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                                                                                                  <  6:84587/149‑1 (MQ=255)
           gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                                                                                                  >  7:109553/1‑149 (MQ=255)
           gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                                                                                                  <  6:46067/149‑1 (MQ=255)
           gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                                                                                                  <  7:46721/149‑1 (MQ=255)
                 cGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcaca                                                                                                                            <  7:280350/149‑1 (MQ=255)
                   gACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATAGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAgg                                                                                                                          <  6:6681/149‑1 (MQ=255)
                    aCAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCAACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGCGGCGTTTTTTATTCACAGGc                                                                                                                         >  3:119584/1‑149 (MQ=255)
                     cAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGttttt                                                                                                                                    <  1:325452/137‑1 (MQ=255)
                        ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAACACGACAGATAGCACAAACGCCCCACTAGGGGGCGTGTTTTATTACCAGGCTAGc                                                                                                                     >  2:290469/1‑149 (MQ=255)
                        ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGc                                                                                                                     >  3:47450/1‑149 (MQ=255)
                        ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGc                                                                                                                     <  3:15233/149‑1 (MQ=255)
                           gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTa                                                                                                                  >  3:248352/1‑149 (MQ=255)
                             cAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTTATTCACAGGCTAGCGTa                                                                                                                  <  8:129526/148‑1 (MQ=255)
                             cAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTAt                                                                                                                 <  2:204801/148‑1 (MQ=255)
                              aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTAt                                                                                                                 >  3:274904/1‑147 (MQ=255)
                              aaTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                                                                                                               <  6:209929/149‑1 (MQ=255)
                               aTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAg                                                                                                              <  7:266890/149‑1 (MQ=255)
                                ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc                                                                                                             >  3:147566/1‑149 (MQ=255)
                                  cccACACCGCCCCGGCCCTTGGCCGTGGCGATTGCGGCGCTGGTTTCACGCCCAGTTTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAACACTACAGATAGCAAAACCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                                                                                                           <  2:171336/149‑1 (MQ=255)
                                   ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                                                                                                           <  8:109553/148‑1 (MQ=255)
                                    cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                                                                                                         >  8:82255/1‑149 (MQ=255)
                                     ccatcGCCACGGCCCAGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                                                                                                         >  3:268713/5‑148 (MQ=255)
                                     acacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                                                                                                        <  2:310167/149‑1 (MQ=255)
                                      cacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGCTTACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCg                                                                                                          <  1:115959/149‑1 (MQ=255)
                                        ccGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                                                                                                     >  8:333938/1‑149 (MQ=255)
                                          gCCCCGGCCCTGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                                                                                                     >  5:74243/1‑147 (MQ=255)
                                          gCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                                                                                                    >  2:260522/1‑148 (MQ=255)
                                          gCCACGGCCCAGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAAG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                                                                                                    <  1:260522/149‑1 (MQ=255)
                                              cGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag                                                                                               >  6:16589/1‑149 (MQ=255)
                                                   cTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttt                                                                                          <  4:147566/149‑1 (MQ=255)
                                                    cgggCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                                         <  8:164425/148‑1 (MQ=255)
                                                      gggggTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGCCAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                       >  7:46073/5‑149 (MQ=255)
                                                      ggCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                       <  8:145904/149‑1 (MQ=255)
                                                        ccGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTgg                                                                                      >  4:203128/1‑148 (MQ=255)
                                                              cGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAgctgc                                                                               <  4:211974/149‑1 (MQ=255)
                                                                    cggcgATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGAt                                                                         <  7:348659/149‑1 (MQ=255)
                                                                       cgATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACg                                                                      <  4:268713/149‑1 (MQ=255)
                                                                        gATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACgg                                                                     <  7:333938/149‑1 (MQ=255)
                                                                             tttCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGt                                                                >  7:11053/1‑149 (MQ=255)
                                                                               tCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCaa                                                                                                      >  7:117200/1‑109 (MQ=255)
                                                                               tCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCaa                                                                                                      <  8:117200/109‑1 (MQ=255)
                                                                                 aCGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGAt                                                             <  7:226157/148‑1 (MQ=255)
                                                                                 aCGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATg                                                            >  7:326634/1‑149 (MQ=255)
                                                                                  cGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGt                                                           <  4:274904/149‑1 (MQ=255)
                                                                                      caGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTAATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTAcc                                                        >  1:292956/1‑148 (MQ=255)
                                                                                          gtTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCAccgc                                                   >  1:134697/1‑149 (MQ=255)
                                                                                            tCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc                                                                                                             >  1:173937/1‑89 (MQ=255)
                                                                                            tCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc                                                                                                             <  2:173937/89‑1 (MQ=255)
                                                                                            tCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCActgctg                                                 >  2:14603/1‑149 (MQ=255)
                                                                                              aTGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCg                                               <  5:16589/149‑1 (MQ=255)
                                                                                               tGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGa                                              <  2:320588/149‑1 (MQ=255)
                                                                                                 cTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGaaa                                            <  4:248352/149‑1 (MQ=255)
                                                                                                  tCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAAt                                           >  8:341091/1‑149 (MQ=255)
                                                                                                   cGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGa                                                              <  5:192422/129‑1 (MQ=255)
                                                                                                   cGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGa                                                              >  6:192422/1‑129 (MQ=255)
                                                                                                      aaGCCTCTACG‑‑‑ACAAAAACGAAAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGAc                                       <  8:156320/149‑1 (MQ=255)
                                                                                                       aGCCTCTACG‑‑‑TCAAATACGAAAGATATCAAAATCGGCCCTCTAGTGGTCGTTTTTTATTCACAGGCTAGCTTATTAGCCCGTCATGCTGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAATTCGACa                                      <  5:187129/149‑1 (MQ=255)
                                                                                                        gCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAg                                     <  7:82255/149‑1 (MQ=255)
                                                                                                          ctctACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCAAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTAGGTAGCTTCTTAGATAGGGCGAGTGATGTACAACTGCTGCGAAATCGACAgg                                    >  4:354441/1‑148 (MQ=255)
                                                                                                          ctctACG‑‑‑ACAAAAACGACACATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGc                                   <  6:178713/149‑1 (MQ=255)
                                                                                                             tACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACGGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgt                                >  3:73723/1‑149 (MQ=255)
                                                                                                               cG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTNCCACTGCTGCGAAATCGACAGGCAgttg                              <  1:213798/149‑1 (MQ=255)
                                                                                                               cG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttg                              <  3:203128/149‑1 (MQ=255)
                                                                                                                g‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGCCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTATTTAGGTAGCTGCTTCGATACGGCGAGTGATGTACCAATGCTGCGAAAAGGACAGGAAgtttt                             >  5:319551/1‑147 (MQ=255)
                                                                                                                     cAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTc                           <  5:193751/149‑1 (MQ=255)
                                                                                                                      aaaaaCGACAGATAGCAAACACGCCCCAGTAGTGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTca                          <  7:141021/149‑1 (MQ=255)
                                                                                                                      aaaaaCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTca                          <  4:47450/149‑1 (MQ=255)
                                                                                                                        aaaCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgt                                >  7:227974/1‑141 (MQ=255)
                                                                                                                        aaaCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgt                                <  8:227974/141‑1 (MQ=255)
                                                                                                                          aCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAAc                      >  3:149474/1‑149 (MQ=255)
                                                                                                                              cAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAg                  <  1:206008/149‑1 (MQ=255)
                                                                                                                              cAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAg                  >  3:304515/1‑149 (MQ=255)
                                                                                                                               aGATAGCATAAACGCCCCTCTAGGTGGCGTTTTTTATTCACAGGCTAGCGTTTTAGCCCGTCAAGCAGCAGCTTTTTTGGTTTCTGCTTCGATACGGCGAGTGATGTTCCACTGGTGCGAAATCTAGAGGCTGTTGTTCACAACCCAg                  <  4:119584/148‑1 (MQ=255)
                                                                                                                                 aTAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTAc               <  8:183430/149‑1 (MQ=255)
                                                                                                                                   ataaaaaaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCGTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCGGCGAAATCGACAGGCAGTTGTTCCCTACCCAGTACAg             >  7:165400/4‑149 (MQ=255)
                                                                                                                                   ataaaaaaCGCCCCCTAGTGGGGCGGTTTTTATTGACAGGGTAGCGTATTAGCCGGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCGGCGAATTCGACAGGCAGTTGTTCACAACCCAGTCCAg             >  7:132869/4‑149 (MQ=255)
                                                                                                                                     aaaaaaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCa           >  6:189352/2‑149 (MQ=255)
                                                                                                                                      aaaaaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCAc          >  4:33657/1‑149 (MQ=255)
                                                                                                                                       aaaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgtt                            <  1:207047/130‑1 (MQ=255)
                                                                                                                                       aaaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAgttgtt                            >  2:207047/1‑130 (MQ=255)
                                                                                                                                        aaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                       <  8:265695/70‑1 (MQ=255)
                                                                                                                                        aaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                       >  7:265695/1‑70 (MQ=255)
                                                                                                                                        aaaCGCCCCCTAGTGGGGCGTTTTTTTTTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCAc          >  5:249353/1‑147 (MQ=255)
                                                                                                                                            gCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGAcc    >  3:332223/1‑149 (MQ=255)
                                                                                                                                             ccccACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCg   <  1:14603/149‑1 (MQ=255)
                                                                                                                                              cccACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCgg  <  8:326634/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                              
CTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACG‑‑‑ACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGG  >  NC_002947/6287‑6574

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