Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,554,138	G→A	6.8%	A56A (GCG→GCA)	sbcB →	exodeoxyribonuclease I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,554,138	0	G	A	6.8%	251.6 / 5.7	87	A56A (GCG→GCA)	sbcB	exodeoxyribonuclease I
Reads supporting (aligned to +/- strand):  ref base G (41/40);  new base A (3/3);  total (44/43)
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.85e-01

TGAAACCACCGGCATCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTCGACGATGAA  >  NC_002947/1554000‑1554282                                                                                                                                           |                                                                                                                                                 tGAAACCACCGGCATCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGt                                                                                                                                       <  2:1368771/150‑1 (MQ=255)      ccaccGGCATCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAg                                                                                                                                  >  1:800652/1‑150 (MQ=255)       caccGGCATCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAgg                                                                                                                                 >  2:433202/1‑150 (MQ=255)            gCATCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGCCTTCGCCATAAACGAACTCGTTAACGCAGTCATCCTTGATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTa                                                                                                                            >  1:2122094/1‑150 (MQ=255)             cATCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTAc                                                                                                                           >  1:1579965/1‑150 (MQ=255)              aTCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCa                                                                                                                               <  2:385095/145‑1 (MQ=255)              aTCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCa                                                                                                                               >  1:385095/1‑145 (MQ=255)               tCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTgg                                                                                                                                        >  2:1112894/1‑135 (MQ=255)               tCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTgg                                                                                                                                        <  1:1112894/135‑1 (MQ=255)               tCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTAcgc                                                                                                                         <  2:363873/150‑1 (MQ=255)                cAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTAcgcc                                                                                                                        <  1:560903/150‑1 (MQ=255)                               tCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGa                                                                                                         <  1:433202/150‑1 (MQ=255)                                 ggCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGc                                                                                                       <  1:1949984/150‑1 (MQ=255)                                  gCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCa                                                                                                      <  1:435195/150‑1 (MQ=255)                                   ccGCTGCAGGTGGCCGGCCTGCGCACTGCCTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGANCAg                                                                                                     <  2:1498307/150‑1 (MQ=255)                                    cGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAgg                                                                                                    <  2:2311618/150‑1 (MQ=255)                                            gcggccggccTGCGCACTGCCCGCGACGGTAACGAAAAAGATGGGCGAGTTAGGCTTTCTGGGCGACCGTCGGACGTCATCCTGCCGCACCCGGCGGCCGGCCGGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCCGGGCCTGTGc                                                                                            <  2:2361243/148‑1 (MQ=255)                                            gTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCAACCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGc                                                                                            >  1:2032842/1‑150 (MQ=255)                                               gccggccTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGaa                                                                                         <  2:1556324/150‑1 (MQ=255)                                                    ccTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGa                                                                                    <  1:319693/150‑1 (MQ=255)                                                      tGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGt                                                                                  >  1:778930/1‑150 (MQ=255)                                                         gcACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCa                                                                               >  2:1480424/1‑150 (MQ=255)                                                            cTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGa                                                                            <  1:989583/150‑1 (MQ=255)                                                            cTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTCTGCGAAGCCGAGTTCATGa                                                                            <  2:154420/150‑1 (MQ=255)                                                                 ttCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGa                                                                                    <  1:1462003/137‑1 (MQ=255)                                                                 ttCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGa                                                                                    >  2:1462003/1‑137 (MQ=255)                                                                 ttCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTcgc                                                                       <  2:1545929/150‑1 (MQ=255)                                                                   cGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGt                                                                     <  2:778930/150‑1 (MQ=255)                                                                   cGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGt                                                                     <  1:1480424/150‑1 (MQ=255)                                                                   cGACCTCAACGAAACCGACGAGCCAATCAGCCTTTATTGCCTACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGt                                                                     <  1:1738439/150‑1 (MQ=255)                                                                       cTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCAc                                                                 >  2:1539914/1‑150 (MQ=255)                                                                        tCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGc                                                                   >  2:245433/1‑147 (MQ=255)                                                                        tCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGc                                                                   <  1:245433/147‑1 (MQ=255)                                                                        tCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACg                                                                >  2:400988/1‑150 (MQ=255)                                                                         cAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGc                                                               >  2:2379873/1‑150 (MQ=255)                                                                            cGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCa                                                            >  2:1832290/1‑150 (MQ=255)                                                                              aaaTCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGc                                                          <  1:400988/150‑1 (MQ=255)                                                                               aaTCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCt                                                         >  2:1762924/1‑150 (MQ=255)                                                                               aaTCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCt                                                         >  2:179439/1‑150 (MQ=255)                                                                                   gacgaGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGgcg                                                     >  2:1790213/1‑150 (MQ=255)                                                                                   gacgaGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCAGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGgcg                                                     >  2:1994587/1‑150 (MQ=255)                                                                                        gCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAAcc                                                >  1:1468518/1‑150 (MQ=255)                                                                                        gCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAAcc                                                >  2:2013127/1‑150 (MQ=255)                                                                                         ccAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGccc                                                             <  2:2016707/136‑1 (MQ=255)                                                                                         ccAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGccc                                                             >  1:2016707/1‑136 (MQ=255)                                                                                           cgacAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCGTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCggg                                             >  1:1264683/4‑150 (MQ=255)                                                                                           aaTCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCggg                                             >  2:2064461/1‑150 (MQ=255)                                                                                            aTCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCAGCCTGCCTGGTGACAGGCATTACGCCGCAActgc                                                                                                                >  1:695356/1‑82 (MQ=255)                                                                                            aTCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCAGCCTGCCTGGTGACAGGCATTACGCCGCAActgc                                                                                                                >  1:695183/1‑82 (MQ=255)                                                                                            aTCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCAGCCTGCCTGGTGACAGGCATTACGCCGCAActgc                                                                                                                <  2:695183/82‑1 (MQ=255)                                                                                            aTCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCAGCCTGCCTGGTGACAGGCATTACGCCGCAActgc                                                                                                                <  2:695356/82‑1 (MQ=255)                                                                                              cAGNCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAANCNGAGTTCATGACTCGCNTGCACGCCCAGCTGGCGCAACCGGGCAc                                          <  1:1762924/150‑1 (MQ=255)                                                                                                 ccTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCtg                                       >  1:2313420/1‑150 (MQ=255)                                                                                                  cTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAActg                                                                                                                 >  1:331470/1‑75 (MQ=255)                                                                                                  cTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAActg                                                                                                                 <  2:331470/75‑1 (MQ=255)                                                                                                   tttATTGCCGACCGTCGGCCGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATAACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGAACCGCGTCCCCCCCCAGCTGGCCCAACCGGGCAccactg                                     >  2:1060265/1‑146 (MQ=255)                                                                                                   tttATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGgc                                                      >  2:992139/1‑133 (MQ=255)                                                                                                   tttATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGgc                                                      <  1:992139/133‑1 (MQ=255)                                                                                                   tttATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCtgtg                                     <  1:1549834/150‑1 (MQ=255)                                                                                                      aTTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTgtggtg                                  >  1:1542267/1‑150 (MQ=255)                                                                                                       ttGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGc                                                                                                  >  2:934459/1‑85 (MQ=255)                                                                                                       ttGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGc                                                                                                  <  1:934459/85‑1 (MQ=255)                                                                                                       ttGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGc                                 <  1:1539914/150‑1 (MQ=255)                                                                                                           cGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGc                             >  1:2193250/1‑150 (MQ=255)                                                                                                            gACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGt                                                                                  >  2:769524/1‑96 (MQ=255)                                                                                                            gACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGt                                                                                  <  1:769524/96‑1 (MQ=255)                                                                                                            gACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCt                            >  2:382117/1‑150 (MQ=255)                                                                                                                gTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTacaa                        <  1:179439/150‑1 (MQ=255)                                                                                                                    gacgacATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACAcc                    <  1:1790213/150‑1 (MQ=255)                                                                                                                       gacATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCtgtg                                     <  1:2406459/130‑1 (MQ=255)                                                                                                                       gacATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCtgtg                                     >  2:2406459/1‑130 (MQ=255)                                                                                                                          aTCCTGCCGCGCCCGGCAGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGAAGCAGGGTTTTTGCGAAGCAGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTgcgc              <  1:1994587/150‑1 (MQ=255)                                                                                                                          aTCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCAt                                                                              >  1:1366138/1‑86 (MQ=255)                                                                                                                          aTCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCAt                                                                              <  2:1366138/86‑1 (MQ=255)                                                                                                                            ccTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCCACAACACCCTGCGCtt            >  2:2155504/1‑150 (MQ=255)                                                                                                                             cTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTTGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTc           >  1:436070/1‑150 (MQ=255)                                                                                                                             cTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTc           <  2:436070/150‑1 (MQ=255)                                                                                                                             cTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTc           >  2:1954139/1‑150 (MQ=255)                                                                                                                             cTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTc           >  2:882593/1‑150 (MQ=255)                                                                                                                              tGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTcg          >  2:1833126/1‑150 (MQ=255)                                                                                                                               gCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTcga         <  1:382117/150‑1 (MQ=255)                                                                                                                                ccGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTcgac        <  1:1411289/150‑1 (MQ=255)                                                                                                                                   cACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTCGACGAt     <  1:2013127/150‑1 (MQ=255)                                                                                                                                      ccGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTCGACGATGaa  <  1:1954139/150‑1 (MQ=255)                                                                                                                                      ccGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTCGACGATGaa  >  1:1712582/1‑150 (MQ=255)                                                                                                                                      ccGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTCGACGATGaa  >  2:1182418/1‑150 (MQ=255)                                                                                                                                      ccGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTCGACGATGaa  <  2:1468518/150‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                 TGAAACCACCGGCATCAACCCGCGCTGTGATCGGCCGCTGCAGGTGGCCGGCCTGCGCACTGACTTCGACCTCAACGAAATCGACGAGCCAATCAGCCTTTATTGCCGACCGTCGGACGACATCCTGCCGCACCCGGCGGCCTGCCTGGTGACAGGCATTACGCCGCAACTGCTGGCGGAGCAGGGCCTGTGCGAAGCCGAGTTCATGACTCGCGTGCACGCCCAGCTGGCGCAACCGGGCACCTGTGGTGCCGGCTACAACACCCTGCGCTTCGACGATGAA  >  NC_002947/1554000‑1554282

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