Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,645,905	A→T	7.0%	L31Q (CTG→CAG)	PP_3212 ←	Rieske 2Fe‑2S family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,645,905	0	A	T	7.0%	210.0 / 5.1	71	L31Q (CTG→CAG)	PP_3212	Rieske 2Fe‑2S family protein
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base T (2/3);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.23e-01

ATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGC  >  NC_002947/3645759‑3646053                                                                                                                                                      |                                                                                                                                                     aTTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAgg                                                                                                                                                    >  1:45519/1‑149 (MQ=255)  ttCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAggg                                                                                                                                                   <  1:24552/149‑1 (MQ=255)  ttCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAgg                                                                                                                                                    <  2:5140/148‑1 (MQ=255)    cccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGt                                                                                                                                                  <  5:48448/148‑1 (MQ=255)    cccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGt                                                                                                                                                  <  5:12507/148‑1 (MQ=255)     ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCCTCGGGCCGCTTCCACACCACCAGCGCTTCACCCAGCAGGGTg                                                                                                                                                 >  4:229013/1‑148 (MQ=255)     ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTg                                                                                                                                                 >  1:352608/1‑148 (MQ=255)     ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTg                                                                                                                                                 >  8:47081/1‑148 (MQ=255)     ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGCACAGTTTTGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTg                                                                                                                                                 <  1:333794/148‑1 (MQ=255)      cAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGaa                                                                                                                                               >  8:228684/1‑149 (MQ=255)       aGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaag                                                                                                                                              <  7:72207/149‑2 (MQ=255)       aGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaag                                                                                                                                              <  1:128358/149‑2 (MQ=255)       aGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAAt                                                                                                                                              >  4:227566/1‑149 (MQ=255)        gCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaagc                                                                                                                                             <  6:261059/149‑3 (MQ=255)           gTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCt                                                                                                                                          <  1:22380/149‑1 (MQ=255)           gTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCt                                                                                                                                          >  2:140578/1‑149 (MQ=255)             gATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTg                                                                                                                                        >  2:302777/1‑149 (MQ=255)             gATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCAATGCCACCGGCGTACCATCCGGACGCTTCCACAACACCAGCGCTTCACCCAGAAGCGTGAAAGGCTTg                                                                                                                                        >  2:302782/1‑149 (MQ=255)               tACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAGTGAGAGCGGTCCCTCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACTCCAGCGCTTCACCCAGCAGGGTGAATGTCtt                                                                                                                                         >  1:330091/1‑149 (MQ=255)                aCCCGCAGGGGATATTGCCCTGAGCACTGACGAAGCCTTTTGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTggg                                                                                                                                      <  8:182638/148‑1 (MQ=255)                  ccGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGcc                                                                                                                                   <  2:177754/149‑1 (MQ=255)                     cAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTc                                                                                                                                <  3:393443/149‑1 (MQ=255)                       ggCGATATTGCCCCCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCtt                                                                                                                              <  1:140578/149‑1 (MQ=255)                           atatTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAg                                                                                                                          <  2:45519/149‑1 (MQ=255)                                gCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCt                                                                                                                      <  1:81370/148‑1 (MQ=255)                                      aGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAg                                                                                                               <  7:47081/149‑1 (MQ=255)                                      aGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAg                                                                                                               >  7:52598/1‑149 (MQ=255)                                       gCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACAATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAgg                                                                                                              >  1:154361/1‑149 (MQ=255)                                        cACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAgg                                                                                                              >  2:149980/1‑148 (MQ=255)                                              cGAAGCCTTTGGAAAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCATCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAgg                                                                                                       >  2:215097/1‑149 (MQ=255)                                               gAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAggg                                                                                                      >  2:302590/1‑149 (MQ=255)                                                         gACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGc                                                                                                     <  3:302762/140‑1 (MQ=255)                                                         gACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGc                                                                                                     >  4:302762/1‑140 (MQ=255)                                                         gACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGa                                                                                             >  1:157494/1‑148 (MQ=255)                                                            tgTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCTCATTGCCACCTGGGTACGATCGGGACGCTTGCACAACAGCAGCGCTTCACCCAGCAGGGTTATTGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATAGCAGAAAc                                                                                          <  8:131855/147‑1 (MQ=255)                                                              tttGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGgc                                                                                       <  7:45090/149‑1 (MQ=255)                                                              tttGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGGCATAGGGAGCAGGGCATACCAGAAAc                                                                                          >  7:131855/1‑146 (MQ=255)                                                                tGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGgcg                                                                                      <  7:88375/148‑1 (MQ=255)                                                                    ggtgcggtgGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATa                                                                                 >  5:206605/1‑149 (MQ=255)                                                                        cggtgGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTggg                                                                                                                                      <  7:70961/92‑1 (MQ=255)                                                                        cggtgGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTggg                                                                                                                                      >  8:70961/1‑92 (MQ=255)                                                                             gcag‑‑‑caTCGGTCCCGCATTGCCACCGGGGTACCAGCGGGACGCTTCCACATCACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGgttgtt                                                                        <  4:100788/149‑1 (MQ=255)                                                                             gcag‑‑‑caGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGgttgtt                                                                        <  2:157494/149‑1 (MQ=255)                                                                              cag‑‑‑cagCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGttt                                                                       <  1:158999/149‑1 (MQ=255)                                                                              cag‑‑‑cagCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGttt                                                                       <  2:352608/149‑1 (MQ=255)                                                                               ag‑‑‑cagcGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTg                                                                      <  7:228684/149‑1 (MQ=255)                                                                                      gcGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTa                                                                   >  8:17675/1‑148 (MQ=255)                                                                                         gTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAgg                                                                             <  8:321098/135‑1 (MQ=255)                                                                                         gTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAgg                                                                             >  7:321098/1‑135 (MQ=255)                                                                                          tCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCaa                                                              <  3:227566/149‑1 (MQ=255)                                                                                              gCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACAt                                                           >  8:49115/1‑148 (MQ=255)                                                                                                  tGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAAt                                                        >  6:194662/1‑147 (MQ=255)                                                                                                   gCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCa                                                      <  8:52598/148‑1 (MQ=255)                                                                                                   gCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAg                                                     >  5:128135/1‑149 (MQ=255)                                                                                                    ccACCGGGGTACCATCGGGACTCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGTTTGTTTGGTAACCAACATAATCat                                                     <  1:302782/148‑2 (MQ=255)                                                                                                    ccACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAg                                                     <  1:302777/148‑1 (MQ=255)                                                                                                        cGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTc                                                <  7:17675/149‑1 (MQ=255)                                                                                                            gTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTca                                            <  1:179754/149‑1 (MQ=255)                                                                                                                    gggACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGtt                                     <  1:302590/148‑1 (MQ=255)                                                                                                                    gggACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGtt                                     >  5:43179/1‑148 (MQ=255)                                                                                                                      gACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAg                                                     >  6:47509/1‑130 (MQ=255)                                                                                                                      gACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAg                                                     <  5:47509/130‑1 (MQ=255)                                                                                                                             ccACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTg                              >  7:271195/1‑146 (MQ=255)                                                                                                                             ccACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTg                              <  8:271195/146‑1 (MQ=255)                                                                                                                             ccACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGt                             <  4:181341/147‑1 (MQ=255)                                                                                                                                   cacCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGgtgt                     >  7:151635/1‑149 (MQ=255)                                                                                                                                      cAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGtt                    >  7:134172/1‑147 (MQ=255)                                                                                                                                      cAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGtt                    >  6:169619/1‑147 (MQ=255)                                                                                                                                       agccattCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAg                 >  7:213226/6‑149 (MQ=255)                                                                                                                                           cTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATCCAGGTTGTTTGGTTTCCACAATAATCAGGATTCCTCACAGCGTTGCGCTTGTTATGGTGTTCAGCTGc             >  5:230707/1‑149 (MQ=255)                                                                                                                                           attCACCCTGCTGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGc             >  4:350498/2‑149 (MQ=255)                                                                                                                                                  cagcagGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGag      >  8:362106/1‑148 (MQ=255)                                                                                                                                                  cagcagGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGAc      <  6:206605/149‑1 (MQ=255)                                                                                                                                                      agGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGc  >  6:357450/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     ATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGACCTGC  >  NC_002947/3645759‑3646053

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