Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,645,902	A→T	7.0%	L32Q (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,902	0	A	T	7.0%	194.3 / 4.6	71	L32Q (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 = 4.63e-01

GTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGAC  >  NC_002947/3645757‑3646049                                                                                                                                                     |                                                                                                                                                    gTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCcag                                                                                                                                                    >  6:48448/1‑147 (MQ=255)   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                                                                                                                                             >  8:47081/1‑148 (MQ=255)       ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTg                                                                                                                                             >  1:352608/1‑148 (MQ=255)       ccAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGCACAGTTTTGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTg                                                                                                                                             <  1:333794/148‑1 (MQ=255)        cAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGaa                                                                                                                                           >  8:228684/1‑149 (MQ=255)         aGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaag                                                                                                                                          <  1:128358/149‑2 (MQ=255)         aGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCTGCTGGGTGaag                                                                                                                                          <  7:72207/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:52598/1‑149 (MQ=255)                                        aGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAg                                                                                                           <  7:47081/149‑1 (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                                                                         >  7:321098/1‑135 (MQ=255)                                                                                           gTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAgg                                                                         <  8:321098/135‑1 (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                                                 <  5:47509/130‑1 (MQ=255)                                                                                                                        gACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAg                                                 >  6:47509/1‑130 (MQ=255)                                                                                                                               ccACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTg                          <  8:271195/146‑1 (MQ=255)                                                                                                                               ccACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTg                          >  7:271195/1‑146 (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)                                                                                                                                                     |                                                                                                                                                    GTATTCCCAGCCGTGATACCCGCAGGCGATATTGCCCTCAGCACTGACGAAGCCTTTGGACAGTTTGGCGGTGCGGTGGCAG‑‑‑CAGCGGTCCCGCATTGCCACCGGGGTACCATCGGGACGCTTCCACAACACCAGCGCTTCACCCAGCAGGGTGAATGGCTTGGGGCCGTCTTCGAGGGCTGTCATAGGGAGCAGGGCATACCAGAAACGGCGCAATACAGGTTGTTTGGTAACCAACATAATCAGGATTCCTCACAGCGTTGAGCTTGTTATGGTGTTCAGCGGCTGCCGAC  >  NC_002947/3645757‑3646049

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