Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	58,437	T→A	10.2%	intergenic (‑348/+340)	PP_5424 ← / ← PP_0050	hypothetical protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	58,437	0	T	A	10.2%	163.1 / 8.7	59	intergenic (‑348/+340)	PP_5424/PP_0050	hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/25);  new base A (4/2);  total (32/27)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.74e-01

ACCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAAT  >  NC_002947/58289‑58585                                                                                                                                                     |                                                                                                                                                      tccTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCt                                                                                                                                                       <  7:149154/148‑1 (MQ=255)  ccTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTg                                                                                                                                                      <  3:265362/149‑1 (MQ=255)   cTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCt                                                                                                                                                       >  8:149154/1‑147 (MQ=255)            gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTg                                                                                                                                             >  5:13761/1‑148 (MQ=255)            gccggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTg                                                                                                                                             <  6:13761/148‑1 (MQ=255)              cggccATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAag                                                                                                                                          <  7:16258/149‑1 (MQ=255)                      gTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcg                                                                                                                                  >  1:168530/1‑149 (MQ=255)                       tCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGAcgg                                                                                                                                 <  2:150470/149‑1 (MQ=255)                             aTCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCa                                                                                                                           >  3:193531/1‑149 (MQ=255)                                       gattgaCACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCg                                                                                                                 >  2:167623/1‑149 (MQ=255)                                         ttgaCACGCGTGCCGATTTGGATAGAATTGCAATAGAGGGCTGCAACTGCGCATTTTTCGCGCCTGGTTACCACAATGACCACCAGCAATATTGGCAGAAGAGGAAGCTGATTCCTATGAAGAAGTGACTGCGGGCATAGCCAGACGGa                                                                                                               >  6:154172/1‑149 (MQ=255)                                            acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                           >  5:235327/1‑118 (MQ=255)                                            acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGaa                                                                                                                                           <  6:235327/118‑1 (MQ=255)                                            acacGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTc                                                                                                            >  1:27933/1‑149 (MQ=255)                                              acGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTg                                                                                                          <  8:152214/149‑1 (MQ=255)                                                gcgTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGaa                                                                                                        >  6:273022/1‑149 (MQ=255)                                                    gCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCATTTTTGGCAGCAGAGGAAGCTGGTTGCTCTGAAGTAGTGCCGGCGGGCATTGCCAGGCTGATTCTGAAcc‑c                                                                                                    >  2:265937/1‑148 (MQ=255)                                                            ggATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGc                                                                                            <  2:27933/148‑1 (MQ=255)                                                             gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGc                                                                                            >  5:72980/1‑147 (MQ=255)                                                             gATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCa                                                                                          <  5:144089/149‑1 (MQ=255)                                                                 gCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGAc                                                                                       <  5:273022/148‑1 (MQ=255)                                                                  ccTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGAc                                                                                       >  8:227244/1‑147 (MQ=255)                                                                   cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                           <  5:250333/95‑3 (MQ=255)                                                                   cTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAACCAGCTTCCTCTgct                                                                                                                                           >  6:250333/1‑93 (MQ=255)                                                                    ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                                    >  6:274644/1‑148 (MQ=255)                                                                    ttGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                                    >  8:87790/1‑148 (MQ=255)                                                                      gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAgg                                                                                   <  2:168530/147‑1 (MQ=255)                                                                      gCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                                    <  7:87790/146‑1 (MQ=255)                                                                             aGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAg                                                                           >  3:104735/1‑148 (MQ=255)                                                                             aGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAg                                                                           <  5:22532/148‑1 (MQ=255)                                                                                  tGCAACTGCGCATTCGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAtt                                                                                                             <  1:54880/110‑1 (MQ=255)                                                                                  tGCAACTGCGCATTCGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGAtt                                                                                                             >  2:54880/1‑110 (MQ=255)                                                                                      aCTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCg                                                                 <  4:104735/149‑1 (MQ=255)                                                                                        tGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                                    <  8:281484/128‑1 (MQ=255)                                                                                        tGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAg                                                                                    >  7:281484/1‑128 (MQ=255)                                                                                           gcATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                                     <  5:143613/61‑1 (MQ=255)                                                                                           gcATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTgg                                                                                                                                                     >  6:143613/1‑61 (MQ=255)                                                                                             aTTTGTCGCGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGCGCATTGCCAGTCGGATTCTGAACC‑AAAGGTCGCCACACGAGGAGCTAAAGGATAGCTCCGACTActc                                                          >  2:90571/1‑149 (MQ=255)                                                                                              tttGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTActcc                                                         <  4:193531/149‑1 (MQ=255)                                                                                                  tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                                     <  6:72980/149‑1 (MQ=255)                                                                                                  tCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCAt                                                     >  4:336597/1‑149 (MQ=255)                                                                                                           ggTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTAcc                                             <  5:274644/148‑1 (MQ=255)                                                                                                             tACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGa                                           >  5:129162/1‑148 (MQ=255)                                                                                                                    aTGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATccc                                                  >  4:340089/1‑134 (MQ=255)                                                                                                                    aTGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATccc                                                  <  3:340089/134‑1 (MQ=255)                                                                                                                           cccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCgtcgg                            <  3:336597/149‑1 (MQ=255)                                                                                                                            ccGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCCCGCgtcggt                           >  1:263307/1‑149 (MQ=255)                                                                                                                             cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                                   >  2:322056/1‑140 (MQ=255)                                                                                                                             cGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCAc                                   <  1:322056/140‑1 (MQ=255)                                                                                                                                      tacttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGc                 >  8:173846/6‑149 (MQ=255)                                                                                                                                      ggCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGc                 >  8:66877/1‑149 (MQ=255)                                                                                                                                        cttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGc                                                                  >  4:221898/4‑98 (MQ=255)                                                                                                                                        cttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGc                                                                  <  3:221898/95‑1 (MQ=255)                                                                                                                                         ttcagAGGAACCAGCTTCCTCTGAAGTAGTGACGGCGGCCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGAGACTACTCCTCATACCTTACCGAAAAGTCAAGCGTCGGTAGTAAACGGCCCt              >  5:97768/3‑149 (MQ=255)                                                                                                                                         agcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACCTGAAGGTAGCCAGACTAGGAGCTAAAGGATAGCTGCGACTACAACGAATCCCTTAACGCAAATTAACGCGGCGGTCGTTAACGGccc               >  1:107067/1‑149 (MQ=255)                                                                                                                                         agcagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCt              <  6:129162/149‑1 (MQ=255)                                                                                                                                           cagAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCttt            >  4:62937/1‑149 (MQ=255)                                                                                                                                              agGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTc           >  1:235294/1‑147 (MQ=255)                                                                                                                                               ggAAGCTGTTTACTCTGAATTATTGACTGCGGGCATTGCCAGTCGGATTCTGAACC‑CAATGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTATCCAg        <  5:154172/149‑1 (MQ=255)                                                                                                                                                 ggtcTGGTTCCTCTGAAGTAGTGAGGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCGCTTTCCAGGc      <  6:249627/146‑1 (MQ=255)                                                                                                                                                    cTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCaa   >  8:97535/1‑149 (MQ=255)                                                                                                                                                     tGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAAt  >  5:325101/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      ACCTCCGTGAGGCCGGCCATCGTCATGCATCAGATCTTGATTGACACGCGTGCCGCTTTGGATAGCCTTGCAATTGAGGGCTGCAACTGCGCATTTGTCGGGCCTGGGTACCAGCATGAGCCCCCGCAATTTTGGCAGCAGAGGAAGCTGGTTCCTCTGAAGTAGTGACGGCGGGCATTGCCAGTCGGATTCTGAACC‑CAAGGTCGCCAGACTAGGAGCTAAAGGATCGCTGCGACTACTCCTCATCCCTTACCGAAAAGTCACGCGTCGGTCGTAAACGGCCCTTTCCAGGCCAAT  >  NC_002947/58289‑58585

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