Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,461,890	A→C	5.9%	I126I (ATA→ATC) ‡	PP_3955 →	permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,461,890	0	A	C	5.9%	253.1 / 4.4	79	I126I (ATA→ATC) ‡	PP_3955	permease
Reads supporting (aligned to +/- strand):  ref base A (32/42);  new base C (3/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCTTGCGGGAATGGCGACGGCTGCATACCAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAGGGCTGTATGGCA  >  NC_002947/4461743‑4462029                                                                                                                                                    |                                                                                                                                            ccTTGCGGGAATGGCGACGGCTGCATACCAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGCCCACCAGTGACCTGATCGTGGGGGTGTTGCTTGTTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAg                                                                                                                                            >  6:89836/1‑149 (MQ=255)     gCGGGAATGGCGACGGCTTCATACCAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGTCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCg                                                                                                                                          <  5:12196/147‑1 (MQ=255)         gAATGGCGACGGCTGCATACCAGTCGTACTGCGAAGGGCATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGccc                                                                                                                                    <  8:256603/149‑1 (MQ=255)              gCGACGGCTGCATACCAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCcatca                                                                                                                                <  3:214252/148‑1 (MQ=255)                 aCGGCTGCATACCAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTg                                                                                                                            <  2:281845/149‑1 (MQ=255)                            ccAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGttt                                                                                                                  <  7:201587/148‑1 (MQ=255)                            ccAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGttt                                                                                                                  <  7:196001/148‑1 (MQ=255)                             cAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCt                                                                                                                >  3:58546/1‑149 (MQ=255)                              aGTGGTACTTCGAACGGGATCTGATTCAACGCTCCGGTGATCTGAACACCAGTGACCTGATCGTTGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTg                                                                                                               <  6:190019/149‑1 (MQ=255)                               gTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTgg                                                                                                              >  7:222318/1‑149 (MQ=255)                                  gTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGCGAATAGCGTTGCCCATCATTTGCGACCTGTTTCTGGcc                                                                                                            >  8:79774/1‑148 (MQ=255)                                      ttCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCcatc                                                                                                                                 <  7:42489/123‑1 (MQ=255)                                      ttCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCcatc                                                                                                                                 >  8:42489/1‑123 (MQ=255)                                      ttCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATgg                                                                                                       >  2:166754/1‑149 (MQ=255)                                        cGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGt                                                                                                     >  4:63507/1‑149 (MQ=255)                                        cGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGt                                                                                                     >  5:108287/1‑149 (MQ=255)                                             gggATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCtt                                                                                                >  4:277807/1‑149 (MQ=255)                                              ggATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTg                                                                                               <  6:312405/149‑1 (MQ=255)                                              ggATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTGGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTg                                                                                               <  7:221096/149‑1 (MQ=255)                                               gATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTgg                                                                                              >  6:223036/1‑149 (MQ=255)                                                 tCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGc                                                                                             >  3:92274/1‑148 (MQ=255)                                                    gATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCAttt                                                                                                                             >  6:90741/1‑113 (MQ=255)                                                    gATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCAttt                                                                                                                             <  5:90741/113‑1 (MQ=255)                                                    gATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCaa                                                                                          >  2:185426/1‑148 (MQ=255)                                                     aTTCAACGCTCCGGTGATCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGc                                                                                             <  6:280225/144‑1 (MQ=255)                                                     aTTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCaa                                                                                          <  4:58546/147‑1 (MQ=255)                                                     aTTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCaa                                                                                          <  3:277807/147‑1 (MQ=255)                                                          aCGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGc                                                                                   <  8:222318/149‑1 (MQ=255)                                                          aCGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGc                                                                                   <  8:105809/149‑1 (MQ=255)                                                            gCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTg                                                                                                   <  5:332964/131‑1 (MQ=255)                                                            gCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTg                                                                                                   >  6:332964/1‑131 (MQ=255)                                                            gCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGc                                                                                   >  5:41523/1‑147 (MQ=255)                                                                 ggTGAGCTGACCACCAGTGACCTGATCGTGGGGCTCTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGACCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCg                                                                             >  5:42262/1‑148 (MQ=255)                                                                   tGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAAc                                                                          >  6:233091/1‑149 (MQ=255)                                                                       cTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTcgc                                                                      <  1:205749/149‑1 (MQ=255)                                                                         gACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTcgcg                                                                     >  7:285663/1‑148 (MQ=255)                                                                          accaccAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCgcgc                                                                    <  3:300971/148‑1 (MQ=255)                                                                              ccaGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATc                                                                 >  5:177484/1‑147 (MQ=255)                                                                              ccaGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGa                                                               <  1:166754/149‑1 (MQ=255)                                                                               caGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAg                                                              <  2:223559/149‑1 (MQ=255)                                                                                    aCCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTa                                                          <  5:193149/148‑1 (MQ=255)                                                                                     ccTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAAc                                                                          <  7:1709/131‑1 (MQ=255)                                                                                     ccTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAAc                                                                          >  8:1709/1‑131 (MQ=255)                                                                                             tGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGAt                                                >  2:111617/1‑149 (MQ=255)                                                                                              ggggCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATc                                               <  6:41523/149‑1 (MQ=255)                                                                                               gggCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCCGTAGCCTCGAtg                                               <  3:63507/148‑2 (MQ=255)                                                                                                 gCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCAttt                                                                                                                             <  2:168354/68‑1 (MQ=255)                                                                                                 gCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCAttt                                                                                                                             >  1:168354/1‑68 (MQ=255)                                                                                                 gCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTAATGGGCTTAGCGTTGCCCATCATTTGGGCCCTGTTTCTGGCCTATGGTTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCCGTAGCCTCGAtgcgc                                            <  6:42262/149‑5 (MQ=255)                                                                                                 gCTGTTGCTTATTGGTTTGGGCTTTGAAGCCGCGCCACGCGTCATGGGCATACCGTGGCCCATCTTTTGCGCCCTGTTTCTGGCCTAGGGGTGGCTTGGCCAACACTTGCCAGGTGAACTCGCGCATCGAGGCTACGGCCTCGATCAGa                                            <  1:242633/149‑1 (MQ=255)                                                                                                   tGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAgg                                                             <  7:355526/130‑1 (MQ=255)                                                                                                   tGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAgg                                                             >  8:355526/1‑130 (MQ=255)                                                                                                     ttGCTTCTTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTATGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATc                                          <  5:89836/147‑1 (MQ=255)                                                                                                       gCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTc                                       <  1:111617/148‑1 (MQ=255)                                                                                                        cTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCaa                                     >  6:262709/1‑149 (MQ=255)                                                                                                        cTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCaa                                     >  8:67454/1‑149 (MQ=255)                                                                                                         ttATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCACCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCaa                                     >  3:211052/1‑148 (MQ=255)                                                                                                           aTTGGGTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAAt                                    <  2:83355/147‑1 (MQ=255)                                                                                                               gtttggtCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGtt                               <  6:177484/148‑1 (MQ=255)                                                                                                               gtttggtCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGtt                               <  7:79774/148‑1 (MQ=255)                                                                                                                    gtCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATgg                                                                                                       >  2:320575/1‑71 (MQ=255)                                                                                                                    gtCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATgg                                                                                                       <  1:320575/71‑1 (MQ=255)                                                                                                                    gtCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCtt                         <  4:92274/149‑1 (MQ=255)                                                                                                                    gtCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCtt                         <  6:169126/149‑1 (MQ=255)                                                                                                                     tCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTc                        <  1:185426/149‑1 (MQ=255)                                                                                                                     tCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTc                        >  4:269981/1‑149 (MQ=255)                                                                                                                       tttGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCTACACTTGCCAGGCGCACTCGCGCATCGAGGCTAGGGCCTCGTTCAGTTCGTCATGCAGTTGGCCTTCgg                      >  3:274858/1‑149 (MQ=255)                                                                                                                                gcgcgACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCa                                              >  4:28294/1‑116 (MQ=255)                                                                                                                                gcgcgACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCa                                              <  3:28294/116‑1 (MQ=255)                                                                                                                                 cgcgACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAggg            <  5:233091/149‑1 (MQ=255)                                                                                                                                       aaatgaTGGGCAACGCTATGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAGGGCTGTAt      >  8:59376/6‑149 (MQ=255)                                                                                                                                         atgaTGGGCAACGCTATGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGa                                                               >  7:165199/4‑90 (MQ=255)                                                                                                                                         atgaTGGGCAACGCTATGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGa                                                               <  8:165199/87‑1 (MQ=255)                                                                                                                                          tCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGc                                                                              <  2:20444/74‑1 (MQ=255)                                                                                                                                          tCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGc                                                                              >  1:20444/1‑74 (MQ=255)                                                                                                                                           gaTGGGCAACGCTATGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAgg             <  1:3141/137‑1 (MQ=255)                                                                                                                                           gaTGGGCAACGCTATGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAgg             >  2:3141/2‑138 (MQ=255)                                                                                                                                                 cATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAGGGCTGTATGGca  <  7:229481/143‑1 (MQ=255)                                                                                                                                                 cATAGCGTTGCCCATCATTTGCGCCCTCTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAGGGCTGTATGGca  >  8:229481/1‑143 (MQ=255)                                                                                                                                                    |                                                                                                                                            CCTTGCGGGAATGGCGACGGCTGCATACCAGTGGTACTTCGAAGGGGATCTGATTCAACGCTCCGGTGAGCTGACCACCAGTGACCTGATCGTGGGGCTGTTGCTTATTGGTTTGGTCTTTGAAGCCGCGCGACGCGTCATGGGCATAGCGTTGCCCATCATTTGCGCCCTGTTTCTGGCCTATGGGTTGCTTGGCCAACACTTGCCAGGCGAACTCGCGCATCGAGGCTACGGCCTCGATCAGATCGTCAATCAGTTGGCCTTCGGCACCGAAGGGCTGTATGGCA  >  NC_002947/4461743‑4462029

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