Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,356,815	A→C	6.3%	V649G (GTT→GGT)	infB ←	translation initiation factor IF‑2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,356,815	0	A	C	6.3%	196.6 / 4.6	64	V649G (GTT→GGT)	infB	translation initiation factor IF‑2
Reads supporting (aligned to +/- strand):  ref base A (31/29);  new base C (2/2);  total (33/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

TGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTCGGCAACCAC  >  NC_002947/5356667‑5356963                                                                                                                                                     |                                                                                                                                                     tGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGaa                                                                                                                                                      >  6:59901/1‑149 (MQ=255) tGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGaa                                                                                                                                                      >  2:5320/1‑149 (MQ=255) tGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGaa                                                                                                                                                      >  2:22198/1‑149 (MQ=255)        gAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGtt                                                                                                                                               >  8:164652/1‑149 (MQ=255)        gAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGtt                                                                                                                                               >  1:242541/1‑149 (MQ=255)         aaGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTg                                                                                                                                              <  1:106376/149‑1 (MQ=255)          aGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGtt                                                                                                                                               >  8:132549/1‑147 (MQ=255)           gCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGTTCTCAAGACCGATGtgcgt                                                                                                                                            <  8:79271/149‑5 (MQ=255)           gCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                            <  5:131653/149‑1 (MQ=255)              aGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                            <  3:249891/146‑1 (MQ=255)              aGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAg                                                                                                                                            >  4:249891/1‑146 (MQ=255)               gcgcCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGtt                                                                                                                                        >  5:295165/1‑149 (MQ=255)                cgcCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACGAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGttt                                                                                                                                       >  2:223497/1‑149 (MQ=255)                  ccAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCTCTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTATTGAGAACGATGTTGAGGGTTttc                                                                                                                                     <  6:208376/149‑1 (MQ=255)                           cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAgg                                                                                                                                           >  1:50199/1‑134 (MQ=255)                           cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAgg                                                                                                                                           <  2:50199/134‑1 (MQ=255)                           cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTtcttct                                                                                                                                 <  6:137556/144‑1 (MQ=255)                           cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTtcttct                                                                                                                                 >  5:137556/1‑144 (MQ=255)                           cATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTcc                                                                                                                            >  3:112358/1‑149 (MQ=255)                                ctctCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGAcc                                                                                                                       >  1:172332/1‑149 (MQ=255)                                 tctcGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGAcc                                                                                                                       >  3:156558/1‑148 (MQ=255)                                   tcGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCAt                                                                                                                    >  5:118966/1‑149 (MQ=255)                                     ggTTATACCAACGACGCCACCACCAATCACGCGCACTTGCAATTAGTCGTTGCCCACGCCGCCGACCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                                   <  4:13842/148‑1 (MQ=255)                                     ggTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                                   <  2:172332/148‑1 (MQ=255)                                       tGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGGCCCATGGtc                                                                                                                >  2:22235/1‑149 (MQ=255)                                         aTACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGtct                                                                                                               <  4:112358/148‑1 (MQ=255)                                         aTACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGtct                                                                                                               <  7:132549/148‑1 (MQ=255)                                             caccGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCg                                                                                                             >  7:250961/1‑146 (MQ=255)                                               ccGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGGCCCATGGTCTCGAAGAt                                                                                                        <  1:22235/149‑1 (MQ=255)                                                  aCGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGtt                                                                                                     <  6:118966/149‑1 (MQ=255)                                                    gCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCg                                                                                                             >  6:13536/1‑139 (MQ=255)                                                    gCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCg                                                                                                             <  5:13536/139‑1 (MQ=255)                                                      caccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTc                                                                                                  >  3:95915/1‑148 (MQ=255)                                                       accaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTc                                                                                                  >  7:259197/1‑147 (MQ=255)                                                       accaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCa                                                                                                <  5:59901/149‑1 (MQ=255)                                                        ccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTttct                                                                                                                                    <  4:243316/112‑1 (MQ=255)                                                        ccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                                   <  8:231448/129‑1 (MQ=255)                                                        ccaccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATg                                                                                                                   >  7:231448/1‑129 (MQ=255)                                                           ccaATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCtt                                                                                            <  2:242541/149‑1 (MQ=255)                                                                 aCGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCgg                                                                                       >  7:176429/1‑148 (MQ=255)                                                                     gcacttgcacttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGa                                                                                  <  6:295165/149‑1 (MQ=255)                                                                       acttgcacttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAg                                                                                 >  5:321316/1‑148 (MQ=255)                                                                             acttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGc                                                                          >  4:73767/1‑149 (MQ=255)                                                                               ttCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGACGTGAGCACGGGc                                                                          >  5:131370/1‑147 (MQ=255)                                                                                 cgtcgtTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGc                                                                      >  8:238592/1‑149 (MQ=255)                                                                                      ttGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGAc                                                                 <  1:223497/149‑1 (MQ=255)                                                                                           cAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTcgc                                                            <  4:95915/149‑1 (MQ=255)                                                                                             ggCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGtt                                                                                                     <  8:293314/106‑1 (MQ=255)                                                                                             ggCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGtt                                                                                                     >  7:293314/1‑106 (MQ=255)                                                                                                      gCGAACCCTGCAGCGCTTCCAGACAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGATTGCCGCCGTGAGCACGGGCCAGCTTGACATCGCGATACTTGCCt                                                 >  4:158723/1‑149 (MQ=255)                                                                                                         aaCCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCATGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGa                                              <  8:176429/149‑1 (MQ=255)                                                                                                             cTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGAAgg                                           >  7:294956/1‑148 (MQ=255)                                                                                                                 aGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAg                                      <  7:164652/149‑1 (MQ=255)                                                                                                                    gcTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGc                                   <  1:106651/149‑1 (MQ=255)                                                                                                                        ccAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATATTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACt                               >  6:103052/1‑149 (MQ=255)                                                                                                                           gagaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGcc                                                  <  5:127208/127‑1 (MQ=255)                                                                                                                           gagaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGcc                                                  >  6:127208/1‑127 (MQ=255)                                                                                                                            agaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCAc                           >  2:6111/1‑149 (MQ=255)                                                                                                                            agaACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCAc                           <  6:321316/149‑1 (MQ=255)                                                                                                                                   tcaacaTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGa                                              >  1:218788/4‑123 (MQ=255)                                                                                                                                   tcaacaTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGa                                              <  2:218788/120‑1 (MQ=255)                                                                                                                                      acaTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttctt                 <  7:238592/149‑1 (MQ=255)                                                                                                                                      acaTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttctt                 <  1:6111/149‑1 (MQ=255)                                                                                                                                       caTCGTTCTCAAGACCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttctt                 >  3:55650/1‑148 (MQ=255)                                                                                                                                       caTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttcttc                >  7:228703/1‑149 (MQ=255)                                                                                                                                       caTCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTttcttc                <  1:22198/149‑1 (MQ=255)                                                                                                                                            gTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAAATTCACGGGCTTTCTTCtcgtc           >  6:182978/1‑149 (MQ=255)                                                                                                                                                     aCGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTCGGCAACCAc  <  8:250961/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     TGCATTCGAAGCCAGCGCCAGGTTGGCATCGCTCTCGGTGATACCACCGACGCCACCACCAATCACGCGCACTTGCACTTCGTCGTTGCCCAGGCCGCCGAGCGAACCCTGCAGCGCTTCCAGAGAACCACGCACATCGGTCTTGAGAACGATGTTGAGGGTTTTCTTCTCTTCCTGACCCATGGTCTCGAAGATGTTTTCCAGCTTGCCGGCGTGAGCACGGGCCAGCTTGACCTCGCGATACTTGCCTTGACGGAACAGGGCAACTTCACGGGCTTTCTTCTCGTCGGCAACCAC  >  NC_002947/5356667‑5356963

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