Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,339	T→A	7.1%	Q15L (CAG→CTG)	trmE ←	GTPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,339	0	T	A	7.1%	184.4 / 4.8	61	Q15L (CAG→CTG)	trmE	GTPase
Reads supporting (aligned to +/- strand):  ref base T (24/32);  new base A (2/2);  total (26/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

AATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTT  >  NC_002947/6191‑6488                                                                                                                                                       |                                                                                                                                                       aaTCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGTAATGGACCGGAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGCGATCCACA‑GCCGGCCTGGTCCAGCTAGACGTTCTGATTACCTGACAATTCCCACACCGCCCCGGCCC‑Tg                                                                                                                                                       >  2:1592469/1‑150 (MQ=255)   taccTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGACGGCCTTG‑GCAGCCAGAGGGCCGGATAACCTGACCATTCCCACACCGCCCCGGCCC‑Tggg                                                                                                                                                     >  2:204140/3‑150 (MQ=255)   tCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑Tggg                                                                                                                                                     <  2:918245/150‑1 (MQ=255)    cccTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGc                                                                                                                                                    <  2:271800/150‑1 (MQ=255)     ccTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGcc                                                                                                                                                   >  2:1216973/1‑150 (MQ=255)        tCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTg                                                                                                                                                <  2:1571008/150‑1 (MQ=255)          gTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGc                                                                                                                                              <  2:672023/150‑1 (MQ=255)               gCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATgg                                                                                                                                         >  1:1110368/1‑150 (MQ=255)                cACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCAGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCACGGCCC‑AGGGCCGGGGCGGTgtg                                                                                                                                        <  2:85188/150‑3 (MQ=255)                 aCCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACATCGCCACGGCCCAGGGGCCGGGGCggtgtg                                                                                                                                        <  2:660909/150‑6 (MQ=255)                           cATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGttt                                                                                                                             <  2:1110368/150‑1 (MQ=255)                                tcgCGGAATGGATCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGc                                                                                                                        >  2:828428/1‑150 (MQ=255)                                    ggAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGTCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGc                                                                                                                                                    <  1:1408493/118‑1 (MQ=255)                                    ggAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGTCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGc                                                                                                                                                    >  2:1408493/1‑118 (MQ=255)                                              gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                >  1:397872/1‑144 (MQ=255)                                              gTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGtt                                                                                                                <  2:397872/144‑1 (MQ=255)                                               tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                         <  1:240670/150‑1 (MQ=255)                                               tAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTc                                                                                                         <  1:1216973/150‑1 (MQ=255)                                                     gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                               >  1:974361/1‑138 (MQ=255)                                                     gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                               <  2:974361/138‑1 (MQ=255)                                                     gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                                   <  2:1554029/150‑1 (MQ=255)                                                     gCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                                   >  1:649861/1‑150 (MQ=255)                                                        tGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACACCCTGACAAATACCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGTCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCctc                                                                                                <  1:524011/150‑1 (MQ=255)                                                           tgcgGGGTGAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                             <  1:1096010/149‑1 (MQ=255)                                                           cgcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                             <  2:1539124/150‑1 (MQ=255)                                                           cgcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTAc                                                                                             <  2:1205701/150‑1 (MQ=255)                                                            gcgGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACg                                                                                            >  2:1060617/1‑150 (MQ=255)                                                                       gTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGAc                                                                                 <  1:662671/150‑1 (MQ=255)                                                                                      aaCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                                  >  2:1273055/1‑150 (MQ=255)                                                                                      aaCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                                  >  1:1173728/1‑150 (MQ=255)                                                                                      aaCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGc                                                                  >  2:1153201/1‑150 (MQ=255)                                                                                       aCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                                 >  1:335093/1‑150 (MQ=255)                                                                                       aCAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGcc                                                                 >  1:996866/1‑150 (MQ=255)                                                                                           gccggccTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                               >  1:591483/1‑100 (MQ=255)                                                                                           gccggccTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTc                                                                                                               <  2:591483/100‑1 (MQ=255)                                                                                                  ttG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑gggg                                                     >  2:592273/1‑150 (MQ=255)                                                                                                      gCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCC‑CTAGTGGGGCGt                                                  <  1:381709/150‑1 (MQ=255)                                                                                                         gCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGtttt                                               <  2:649861/150‑1 (MQ=255)                                                                                                            agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GTGGCGTTTTTTa                                            <  1:1060617/150‑1 (MQ=255)                                                                                                            agagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                            <  1:322473/150‑1 (MQ=255)                                                                                                             gagGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTAt                                           <  1:1061592/150‑1 (MQ=255)                                                                                                                   cGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAg                                     <  1:98430/150‑1 (MQ=255)                                                                                                                    ggACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑g                                                        <  1:49992/130‑1 (MQ=255)                                                                                                                    ggACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAg‑g                                                        >  2:49992/1‑130 (MQ=255)                                                                                                                         aCCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGc                               <  2:996866/150‑1 (MQ=255)                                                                                                                            tGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTa                            <  2:1173728/150‑1 (MQ=255)                                                                                                                             gACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                                     >  1:723282/1‑61 (MQ=255)                                                                                                                             gACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGcaca                                                                                                                     <  2:723282/61‑1 (MQ=255)                                                                                                                             gACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTAt                           >  2:243027/1‑150 (MQ=255)                                                                                                                               cAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTa                         <  1:828428/150‑1 (MQ=255)                                                                                                                               cAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTa                         <  1:1592469/150‑1 (MQ=255)                                                                                                                                 aTTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGc                       <  1:592273/150‑1 (MQ=255)                                                                                                                                      cacaCCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTc                  <  1:243027/150‑1 (MQ=255)                                                                                                                                       acacCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCa                 >  1:123075/1‑150 (MQ=255)                                                                                                                                         atcGCCACGGCCC‑AGGGCCGGGGCGGTGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAag               >  2:575248/3‑150 (MQ=255)                                                                                                                                         atcGCCACGGCCC‑AGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTa                                            <  2:84631/119‑1 (MQ=255)                                                                                                                                         atcGCCACGGCCC‑AGGGCCGGGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGACGGTTTTTa                                            >  1:84631/3‑121 (MQ=255)                                                                                                                                          ccGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                                   <  2:1187961/66‑1 (MQ=255)                                                                                                                                          ccGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGc                                                                                                   >  1:1187961/1‑66 (MQ=255)                                                                                                                                               ccGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAagcagcag         >  1:428037/1‑150 (MQ=255)                                                                                                                                                    cc‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttt    >  1:841760/1‑150 (MQ=255)                                                                                                                                                       tGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttttt  <  1:204140/150‑1 (MQ=255)                                                                                                                                                       tGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCtttttt  <  2:335093/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       AATCCCTTCGTCCAGCACCAGCCCTTCATCGTCGCGGAATGGACCGTAATGGGCATGCCGCGGGGTCAGGGTGCGGCCGGTGATCAACAGGCCGGCCTTG‑GCAGCCAGAGGGCCGGACAACCTGACAATTCCCACACCGCCCCGGCCC‑TGGGCCGTGGCGATGGCGGCGATGGTTTCACGCACAGTGTTCATGCTCGAAAGCCTCTACGACAAAAACGACAGATAGCAAAAACGCCCCACTAG‑GGGGCGTTTTTTATTCACAGGCTAGCGTATTAGCCCGTCAAGCAGCAGCTTTTTT  >  NC_002947/6191‑6488

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: ‑