Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,429 Δ1 bp 6.9% intergenic (‑47/+42) trmE ← / ← oxaA GTPase/membrane protein insertase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,4290A. 6.9% 268.9 / 14.0 87intergenic (‑47/+42)trmE/oxaAGTPase/membrane protein insertase
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base . (2/4);  total (44/43)
Fisher's exact test for biased strand distribution p-value = 4.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.68e-01

GCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGCAGG  >  NC_002947/6281‑6578
                                                                                                                                                    |                                                                                                                                                       
gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCAc                                                                                                                                                        >  1:1880246/1‑150 (MQ=255)
gccggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCAc                                                                                                                                                        >  1:1196828/1‑150 (MQ=255)
   ggccTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg                                                                                                                                                    <  1:126363/150‑1 (MQ=255)
            aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGttt                                                                                                                                          >  1:1376419/1‑150 (MQ=255)
            aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCTAGTGGGGCGttt                                                                                                                                          <  2:1648701/150‑1 (MQ=255)
            aGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCTAGTGGGGCGttt                                                                                                                                          <  1:422727/150‑1 (MQ=255)
             gCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                                                                                                         >  2:491065/1‑150 (MQ=255)
                 gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAGACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTAt                                                                                                                                     <  2:597605/150‑1 (MQ=255)
                             aCCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACGAGGGGGCGTTTTTTATTCACAGGCTAGc                                                                                                                         <  1:567135/150‑1 (MQ=255)
                                tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcaca                                                                                                                                >  2:146047/1‑140 (MQ=255)
                                tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcaca                                                                                                                                <  1:146047/140‑1 (MQ=255)
                                tGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTCTTATTCACAGGCTAGCGTa                                                                                                                      <  2:1196828/150‑1 (MQ=255)
                                 gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcac                                                                                                                                 >  2:713176/1‑138 (MQ=255)
                                 gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTcac                                                                                                                                 <  1:713176/138‑1 (MQ=255)
                                 gACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTAt                                                                                                                     <  1:711775/150‑1 (MQ=255)
                                      ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGcc                                                                                                                <  1:1496720/150‑1 (MQ=255)
                                      ttCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGcc                                                                                                                >  1:947675/1‑150 (MQ=255)
                                       tCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                                                                                                               <  2:487784/150‑1 (MQ=255)
                                       tCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGccc                                                                                                               <  2:947675/150‑1 (MQ=255)
                                        ccgccatcGCCACGGCCCAGGGCCGGGGCGGTGTGGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCg                                                                                                              >  2:1694098/8‑150 (MQ=38)
                                         cgccatcGCCACGGCCCAGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCCTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAg                                                                                                                  <  1:1570722/139‑1 (MQ=255)
                                         cgccatcGCCACGGCCCAGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAg                                                                                                                  >  2:1570722/7‑145 (MQ=255)
                                         ccACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGt                                                                                                             >  1:1283040/1‑150 (MQ=255)
                                          cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGttt                                                                                                                                          <  1:1103994/120‑1 (MQ=255)
                                          cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGttt                                                                                                                                          >  2:1103994/1‑120 (MQ=255)
                                          cacaCCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                                                                                                            <  1:557778/150‑1 (MQ=255)
                                           acacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCa                                                                                                           >  2:1369426/1‑150 (MQ=255)
                                           acacCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCa                                                                                                           >  2:1691334/1‑150 (MQ=255)
                                             acCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag                                                                                                         >  1:1154729/1‑150 (MQ=255)
                                              ccGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                                                                                                        <  1:1369426/150‑1 (MQ=255)
                                               cGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagca                                                                                                       >  1:522369/1‑150 (MQ=255)
                                                 ccccGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagc                                                                                                     <  2:1283040/150‑1 (MQ=255)
                                                    cGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGt                                                                                                                                            <  2:292362/108‑1 (MQ=255)
                                                    cGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGt                                                                                                                                            >  1:292362/1‑108 (MQ=255)
                                                        ccTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttt                                                                                              >  1:1594550/1‑150 (MQ=255)
                                                         cTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                                             >  1:1259962/1‑150 (MQ=255)
                                                             gCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGt                                                                                         <  2:522369/150‑1 (MQ=255)
                                                                 tGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                           >  1:844041/1‑144 (MQ=255)
                                                                 tGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                           <  2:844041/144‑1 (MQ=255)
                                                                      aTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                                             >  1:1772805/1‑137 (MQ=255)
                                                                      aTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                                             <  2:1772805/137‑1 (MQ=255)
                                                                       tGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCg                                                                               >  1:1263496/1‑150 (MQ=255)
                                                                            gcgATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAg                                                                                       <  2:1495589/137‑1 (MQ=255)
                                                                            gcgATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAg                                                                                       >  1:1495589/1‑137 (MQ=255)
                                                                                   tttCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTg                                                                   >  1:1260965/1‑150 (MQ=255)
                                                                                   tttCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCTAGTGGGGCGttttt                                                                                                                                        >  1:1931702/1‑81 (MQ=38)
                                                                                   tttCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCCTAGTGGGGCGttttt                                                                                                                                        <  2:1931702/81‑1 (MQ=38)
                                                                                     tCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGAt                                                                 <  2:1594550/150‑1 (MQ=255)
                                                                                        cGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGt                                                               >  1:1594148/1‑150 (MQ=255)
                                                                                         gCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTAc                                                             >  1:1452536/1‑150 (MQ=255)
                                                                                           acaGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                                                                                                         <  1:1074166/72‑1 (MQ=255)
                                                                                           acaGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGtttt                                                                                                                                         >  2:1074166/1‑72 (MQ=255)
                                                                                            caGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCAc                                                          <  2:1376419/150‑1 (MQ=255)
                                                                                                  tCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACtgctgc                                                    <  2:1880246/150‑1 (MQ=255)
                                                                                                    aTGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                                             >  2:1772784/1‑107 (MQ=255)
                                                                                                    aTGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCttttt                                                                                             <  1:1772784/107‑1 (MQ=255)
                                                                                                    aTGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGa                                                  >  2:1750290/1‑150 (MQ=255)
                                                                                                        tCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                                                                                                        >  2:1764858/1‑92 (MQ=255)
                                                                                                        tCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagc                                                                                                        <  1:1764858/92‑1 (MQ=255)
                                                                                                        tCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGa                                                                              <  2:1796752/118‑1 (MQ=255)
                                                                                                        tCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGa                                                                              >  1:1796752/1‑118 (MQ=255)
                                                                                                        tCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTAc                                                             <  1:1831254/135‑1 (MQ=255)
                                                                                                        tCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTAc                                                             >  2:1831254/1‑135 (MQ=255)
                                                                                                            aaGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACa                                          <  1:491065/150‑1 (MQ=255)
                                                                                                             aGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACa                                          <  2:1594148/150‑1 (MQ=255)
                                                                                                               ccTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGc                                       >  1:587629/1‑150 (MQ=255)
                                                                                                                 tctACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                           <  1:1470584/96‑1 (MQ=255)
                                                                                                                 tctACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTg                                                                                           >  2:1470584/1‑96 (MQ=255)
                                                                                                                      gACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                                                                                                            <  2:1304287/74‑1 (MQ=255)
                                                                                                                      gACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                                                                                                            >  1:1304287/1‑74 (MQ=255)
                                                                                                                       aCAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACtgctgc                                                    >  1:1123701/1‑129 (MQ=255)
                                                                                                                       aCAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTc                               <  2:587629/150‑1 (MQ=255)
                                                                                                                       aCAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGCTAGCTGCTTCGATACGGCGAGTGATGTACCACtgctgc                                                    <  2:1123701/129‑1 (MQ=255)
                                                                                                                                aCAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAg                      >  2:28697/1‑150 (MQ=255)
                                                                                                                                    aTAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACa                  >  2:618894/1‑150 (MQ=255)
                                                                                                                                       gCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCTGCAGCTTTTCTCGCCTGCGCTTCGATACCGCGAGTGTTGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCa               <  1:1691334/150‑1 (MQ=255)
                                                                                                                                       gCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAg                                                                                       <  1:1213328/78‑1 (MQ=255)
                                                                                                                                       gCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAg                                                                                       >  2:1213328/1‑78 (MQ=255)
                                                                                                                                        cAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCAc              <  2:1146864/150‑1 (MQ=255)
                                                                                                                                         aaaaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                                                                                                                   <  1:1935964/48‑1 (MQ=255)
                                                                                                                                         aaaaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTa                                                                                                                   >  2:1935964/1‑48 (MQ=255)
                                                                                                                                         aaaaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCAcc             >  2:1765407/1‑150 (MQ=255)
                                                                                                                                          aaaaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAAccc                        >  1:1526544/1‑138 (MQ=255)
                                                                                                                                          aaaaCGCCCCCTAGTGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAAccc                        <  2:1526544/138‑1 (MQ=255)
                                                                                                                                           aaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                                                        <  1:1185641/89‑1 (MQ=255)
                                                                                                                                           aaaCGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGc                                                                        >  2:1185641/1‑89 (MQ=255)
                                                                                                                                                 cccACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGc     >  1:923220/1‑150 (MQ=255)
                                                                                                                                                    aCTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGCAgg  >  2:1722775/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GCCGGCCTTGGCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCCTGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAGGGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTTGGTAGCTGCTTCGATACGGCGAGTGATGTACCACTGCTGCGAAATCGACAGGCAGTTGTTCACAACCCAGTACAGCACCAGACCGGCAGG  >  NC_002947/6281‑6578

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: