Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,577,653	T→A	6.9%	L23Q (CTG→CAG)	PP_3158 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,577,653	0	T	A	6.9%	179.5 / 3.3	58	L23Q (CTG→CAG)	PP_3158	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/28);  new base A (2/2);  total (28/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

CACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTGCTGAT  >  NC_002947/3577505‑3577794                                                                                                                                                        |                                                                                                                                              cACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCt                                                                                                                                               <  1:186146/149‑1 (MQ=255)    cGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCtgtg                                                                                                                                            >  1:109937/1‑149 (MQ=255)       aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                                         >  7:179686/1‑149 (MQ=255)       aCTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTgcgc                                                                                                                                         >  6:20488/1‑149 (MQ=255)           tttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTg                                                                                                                                      >  7:105352/1‑148 (MQ=255)            ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGt                                                                                                                                     <  8:179686/148‑1 (MQ=255)            ttCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTa                                                                                                                                    >  6:160556/1‑149 (MQ=255)             tCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTAc                                                                                                                                   >  3:244080/1‑149 (MQ=255)               ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCTACCCTGTTCCGATGGTTGCAAGGTACTCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                                  <  8:295010/148‑1 (MQ=255)               ccATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACg                                                                                                                                  >  8:270317/1‑148 (MQ=255)                cATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGCGCACAGGCATGCCAGGGCACCTGtcagta                                                                                                                                <  6:219781/149‑6 (MQ=38)                 aTTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACgcag                                                                                                                               <  1:266494/149‑1 (MQ=255)                       gggCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGcagcagcag                                                                                                                         <  2:109937/149‑1 (MQ=255)                          ctctATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGc                                                                                                                      >  4:114196/1‑149 (MQ=255)                             tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGcc                                                                                                                     <  6:256792/147‑1 (MQ=255)                             tATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAAGTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                                   >  1:87417/1‑149 (MQ=255)                               tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCg                                                                                                                    <  2:271396/146‑1 (MQ=255)                               tGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGt                                                                                                                   <  7:39925/147‑1 (MQ=255)                                gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                                        <  4:137270/144‑5 (MQ=18)                                gCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTGCTGCGTACAGGCGCACAGGCATGCCAGGGCACCTGTCAgtagcaggcgg                                                                                                                        >  3:137270/1‑140 (MQ=18)                                     gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATgg                                                                                                            <  4:244080/148‑1 (MQ=255)                                     gCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                           <  6:61465/149‑1 (MQ=255)                                          gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                           <  4:185055/144‑1 (MQ=255)                                          gCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGt                                                                                                           >  3:185055/1‑144 (MQ=255)                                                cTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcg                                                                                                <  7:270317/149‑1 (MQ=255)                                                     tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCGTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                           >  7:218030/1‑149 (MQ=255)                                                     tCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaac                                                                                           >  4:183299/1‑149 (MQ=255)                                                       ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                                         >  6:242441/1‑149 (MQ=255)                                                       ccGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacg                                                                                          <  5:160556/148‑1 (MQ=255)                                                          gcgcAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacggca                                                                                       >  1:379792/1‑148 (MQ=255)                                                              aGCCGCCCAATATCGACCCGTAGCCGAAGGTAGCAAGCTACGCGAATTTCCTACCCCGCCTGCTACTG‑‑‑ACAGTTGCCCTGGCATGCCTGTGCGCCTTTACGCATCAGCAGGTCCGTGACATGGTCAACCAGCTCGGCAACGTCAAAgcc                                                                                  <  7:166119/149‑1 (MQ=255)                                                                      aaTATCGACCCGGTGCCGATGGTTGCAAGTTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCt                                                                          >  6:215466/1‑149 (MQ=255)                                                                         atCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGtt                                                                       >  6:374055/1‑149 (MQ=255)                                                                            gACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCa                                                                    <  8:105352/149‑1 (MQ=255)                                                                              cccGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGTCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGa                                                                  <  2:379792/149‑1 (MQ=255)                                                                                    gCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCg                                                            <  7:72046/149‑1 (MQ=255)                                                                                             tGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTcggcaacgg                                                                                         <  2:281141/111‑1 (MQ=255)                                                                                               cAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTg                                                 <  3:114196/149‑1 (MQ=255)                                                                                                aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                                <  5:215466/149‑1 (MQ=255)                                                                                                aaGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTgg                                                >  2:102131/1‑149 (MQ=255)                                                                                                      cgcgAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGc                                          >  1:1461/1‑149 (MQ=255)                                                                                                          aaTTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTc                                      <  7:120340/149‑1 (MQ=255)                                                                                                                 tccCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                               >  6:15830/1‑149 (MQ=255)                                                                                                                   ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCTCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGTCCGCCTGGCCACGCTGTCGATgcgc                               <  8:138593/147‑1 (MQ=255)                                                                                                                   ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATgcgc                               <  8:218030/147‑1 (MQ=255)                                                                                                                   ccccGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGcgcg                              >  4:365268/1‑148 (MQ=255)                                                                                                                         cTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCt                       >  1:93426/1‑149 (MQ=255)                                                                                                                            ccaCTG‑‑‑ACCCGGGCCCTTGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGGGACATGGGCAACCACCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa                    <  5:144422/147‑1 (MQ=255)                                                                                                                             tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa                    >  2:128631/1‑148 (MQ=255)                                                                                                                             tACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGACTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACAGCCTGGCCACGCTGTCGATGCGCGACAACCTGCa                    >  4:93652/1‑148 (MQ=255)                                                                                                                                  cgtaCAGGCGCACAGGCATGCCAGGGCACCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAg                 >  7:43058/4‑149 (MQ=255)                                                                                                                                     aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                               >  3:69522/1‑84 (MQ=255)                                                                                                                                     aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGc                                                                               <  4:69522/84‑1 (MQ=255)                                                                                                                                     aCAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGc                <  2:1461/147‑1 (MQ=255)                                                                                                                                      cAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGGCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTg             <  8:376579/149‑1 (MQ=255)                                                                                                                                            cccTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTAGctc       >  7:140944/1‑148 (MQ=255)                                                                                                                                            cccTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACctg       <  5:15830/149‑1 (MQ=255)                                                                                                                                                 gCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTGCtgat  >  8:314229/1‑149 (MQ=255)                                                                                                                                                       cTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGCCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCa                                                               >  6:333367/1‑82 (MQ=255)                                                                                                                                                       cTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGCCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCa                                                               <  5:333367/82‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                              CACCGCACTTTTTCCCATTTTCGGGCTCTATGCTTGGCCAGGCCGCGCTGACTCCCGGCGCAGCCGCCCAATATCGACCCGGTGCCGATGGTTGCAAGCTACGCGAATTTCCTCCCCCGCCTGCTACTG‑‑‑ACAGGTGCCCTGGCATGCCTGTGCGCCTGTACGCAGCAGCAGGGCCGTGACATGGTCAACCAGCTCGGCAACGGCAAACCCAGCGAGCTGTTCCAGACCAGCGTCGACCGCCTGGCCACGCTGTCGATGCGCGACAACCTGCAAAGCCTGTACCTGCTGAT  >  NC_002947/3577505‑3577794

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