Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,433:1 +T 6.8% 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 6.8% 285.9 / 14.2 74intergenic (‑51/+38)trmE/oxaAGTPase/membrane protein insertase
Reads supporting (aligned to +/- strand):  ref base . (33/36);  new base T (2/3);  total (35/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

CCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGCAGGGAA  >  NC_002947/6286‑6581
                                                                                                                                                    |                                                                                                                                                    
ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggg                                                                                                                                                    >  2:1008967/1‑150 (MQ=255)
ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggg                                                                                                                                                    <  2:820496/150‑1 (MQ=255)
ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgg                                                                                                                                                     <  2:608917/149‑1 (MQ=255)
ccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAgg                                                                                                                                                     >  1:608917/1‑149 (MQ=255)
    ggCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCg                                                                                                                                                >  1:1124580/1‑150 (MQ=255)
     gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAggggcgtt                                                                                                                                               <  2:111823/150‑5 (MQ=255)
        gCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCTAGTGGGGCGtttt                                                                                                                                            <  2:132762/150‑1 (MQ=255)
           agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                                                                                                         <  1:1179569/150‑1 (MQ=255)
           agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTa                                                                                                                                         <  2:532696/150‑1 (MQ=255)
           agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAAATAGGAAAAACTCCCCACTAGGGGGCGTTTTTTa                                                                                                                                         >  2:1301901/1‑150 (MQ=255)
                  cGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAg                                                                                                                                  <  1:883371/150‑1 (MQ=255)
                   ggACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAgg                                                                                                                                 >  2:223712/1‑150 (MQ=255)
                         ccTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCg                                                                                                                           >  2:769030/1‑150 (MQ=255)
                           tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTa                                                                                                                         >  1:28319/1‑150 (MQ=255)
                            gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTAt                                                                                                                        >  2:546312/1‑150 (MQ=255)
                              cagttCCCACATCGCCCCGTGCCTGGGCCGTGGCGATGGCGGCGATGCTGTCACGCACAGTGTTAATGCTCTACAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                                                                                                                      <  2:253635/147‑1 (MQ=255)
                                aTTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc                                                                                                                    >  1:68235/1‑150 (MQ=255)
                                 ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcaca                                                                                                                                   >  1:867266/1‑134 (MQ=255)
                                 ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcaca                                                                                                                                   <  2:867266/134‑1 (MQ=255)
                                     cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                                                                                                               >  2:941389/1‑150 (MQ=255)
                                        acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                                                                                                            <  1:370018/150‑1 (MQ=255)
                                             cccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCAACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagca                                                                                                       <  1:462735/150‑1 (MQ=255)
                                                 gCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtt                                                                                                   <  1:1008967/150‑1 (MQ=255)
                                                    cTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                                                <  2:790146/150‑1 (MQ=255)
                                                       ggCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTgg                                                                                             <  2:28319/150‑1 (MQ=255)
                                                       ggCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTgg                                                                                             <  2:1423628/150‑1 (MQ=255)
                                                        gCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGt                                                                                            <  2:1139432/150‑1 (MQ=255)
                                                        gCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGt                                                                                            >  2:309113/1‑150 (MQ=255)
                                                        gCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGt                                                                                            >  1:1334090/1‑150 (MQ=255)
                                                        gCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGt                                                                                            <  1:1329997/150‑1 (MQ=255)
                                                           gTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAgc                                                                                         >  2:389063/1‑150 (MQ=255)
                                                                gATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCAAAAGCCTCTACGACAAACACTTCAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCtt                                                                                    <  1:945596/150‑1 (MQ=255)
                                                                        cgATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACgg                                                                            <  2:68235/150‑1 (MQ=255)
                                                                        cgATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACgg                                                                            <  2:1073438/150‑1 (MQ=255)
                                                                             gTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                                                           <  2:1323426/146‑1 (MQ=255)
                                                                             gTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                                                           >  1:1323426/1‑146 (MQ=255)
                                                                                tCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGAt                                                                    <  1:309113/150‑1 (MQ=255)
                                                                                   cGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTgctcta                                                                 <  2:1124580/150‑6 (MQ=255)
                                                                                      acaGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCa                                                              <  2:515439/150‑1 (MQ=255)
                                                                                          tgtTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCActgc                                                          <  2:455411/150‑1 (MQ=255)
                                                                                          tgtTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCActgc                                                          <  1:769030/150‑1 (MQ=255)
                                                                                                tGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGaa                                                    >  1:22255/1‑150 (MQ=255)
                                                                                                tGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGaa                                                    >  2:799391/1‑150 (MQ=255)
                                                                                                 gCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGaaa                                                   >  2:135819/1‑150 (MQ=255)
                                                                                                  cTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAAt                                                  <  1:389063/150‑1 (MQ=255)
                                                                                                   tCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATc                                                 >  2:970906/1‑150 (MQ=255)
                                                                                                     gAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttttt                                                                                               <  2:1145763/102‑1 (MQ=255)
                                                                                                     gAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttttt                                                                                               >  1:1145763/1‑102 (MQ=255)
                                                                                                      aaaGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGAc                                              >  1:270555/1‑150 (MQ=255)
                                                                                                          ccTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGc                                          <  2:270555/150‑1 (MQ=255)
                                                                                                           ctctACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTa                                                                                           >  1:363527/1‑100 (MQ=255)
                                                                                                           ctctACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTa                                                                                           <  2:363527/100‑1 (MQ=255)
                                                                                                                    aaaaaCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAACTCCAGTGGTAGCTGCTTCGATACGGCGAGTGATGTACTACTGATCCGAACTCGACACGCAgtggttgag                                >  1:449396/1‑147 (MQ=255)
                                                                                                                      aaaCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACaa                              <  1:1301901/150‑1 (MQ=255)
                                                                                                                        aCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                                                           <  2:655175/103‑1 (MQ=255)
                                                                                                                        aCGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                                                           >  1:655175/1‑103 (MQ=255)
                                                                                                                          gACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCa                          >  1:572489/1‑150 (MQ=255)
                                                                                                                           aCAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCCCTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAg                         <  1:1152728/150‑1 (MQ=255)
                                                                                                                           aCAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACACAg                         >  2:1201/1‑150 (MQ=255)
                                                                                                                             aGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACtgctgc                                                       <  2:263688/118‑1 (MQ=255)
                                                                                                                             aGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACtgctgc                                                       >  1:263688/1‑118 (MQ=255)
                                                                                                                             aGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTa                       >  1:987803/1‑150 (MQ=255)
                                                                                                                             aGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTa                       >  2:387319/1‑150 (MQ=255)
                                                                                                                             aGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTa                       <  1:546312/150‑1 (MQ=255)
                                                                                                                               aTAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTc                                  <  1:49996/137‑1 (MQ=255)
                                                                                                                               aTAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTc                                  >  2:49996/1‑137 (MQ=255)
                                                                                                                                 ataaaaaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGc                   >  1:802261/4‑150 (MQ=255)
                                                                                                                                  taaaaaaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATcagc                                              <  2:991044/120‑4 (MQ=255)
                                                                                                                                  taaaaaaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCggc                                              >  1:991044/3‑120 (MQ=255)
                                                                                                                                   cAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCCTCAAGCAGCAGCTTTTTTGGTGCATGCTTCCATAAAACCACTCACTGACCACGACGGTGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCAc                 <  1:387319/150‑1 (MQ=255)
                                                                                                                                     aaaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCa               >  2:146011/1‑150 (MQ=255)
                                                                                                                                      aaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAg              >  2:883414/1‑150 (MQ=255)
                                                                                                                                       aaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCAc                 <  1:489168/146‑1 (MQ=255)
                                                                                                                                       aaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTATTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCAc                 >  2:489168/1‑146 (MQ=255)
                                                                                                                                            cccACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACaa                              <  1:1289800/128‑1 (MQ=255)
                                                                                                                                            cccACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACaa                              >  2:1289800/1‑128 (MQ=255)
                                                                                                                                            cccACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGc        <  2:22255/150‑1 (MQ=255)
                                                                                                                                               aCTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGCAgg     >  1:757326/1‑150 (MQ=255)
                                                                                                                                                  aGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGCAGGGaa  <  2:987803/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGCAGGGAA  >  NC_002947/6286‑6581

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