Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,104,697	A→G	10.5%	F63L (TTC→CTC)	PP_1879 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,104,697	0	A	G	10.5%	144.2 / 7.3	57	F63L (TTC→CTC)	PP_1879	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (26/25);  new base G (4/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-01

TCAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCG  >  NC_002947/2104549‑2104839                                                                                                                                                     |                                                                                                                                                tCAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGTTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGaa                                                                                                                                                 >  5:100103/1‑149 (MQ=255) tCAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGaa                                                                                                                                                 >  3:122150/1‑149 (MQ=255)     tAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTcc                                                                                                                                             <  2:124625/149‑1 (MQ=255)       aGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCa                                                                                                                                            <  7:101419/148‑1 (MQ=255)        gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCA‑cc                                                                                                                                         >  1:353708/1‑149 (MQ=255)        gCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑cc                                                                                                                                         >  8:10188/1‑149 (MQ=255)             tCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTCCA‑CCCgtgg                                                                                                                                    <  2:353708/149‑4 (MQ=255)               gACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGaa                                                                                                                                                 <  5:117892/135‑1 (MQ=255)               gACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGaa                                                                                                                                                 >  6:117892/1‑135 (MQ=255)               gACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAggg                                                                                                                                  <  7:10188/149‑1 (MQ=255)                 ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTcacc                                                                                                                                           >  5:294490/1‑136 (MQ=255)                 ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTcacc                                                                                                                                           <  6:294490/139‑4 (MQ=255)                 ccGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTACAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAggggg                                                                                                                                <  7:164284/149‑1 (MQ=255)                         aaTATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTAAGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTTCCCACGGGTGGAAATCAA‑CACGCAGGGGGGAATAGGt                                                                                                                        >  1:80283/1‑149 (MQ=255)                              gAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAGTTACA‑CCCgtg                                                                                                                                     >  3:145284/1‑129 (MQ=255)                                   gCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGaa                                                                                                               <  6:264506/148‑1 (MQ=255)                                      ctTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTc                                                                                                            >  7:36982/1‑148 (MQ=255)                                                         gCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTccc                                                                                          <  8:36982/147‑1 (MQ=255)                                                           tGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCATCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCGGGAACTCGACCGTGTTGGGTTTCCCCCCt                                                                                      >  4:211937/1‑149 (MQ=255)                                                            gAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCtt                                                                                     >  3:290186/1‑149 (MQ=255)                                                             aaGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCtt                                                                                     >  3:224667/1‑148 (MQ=255)                                                             aaGGGCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTc                                                                                    <  6:40775/149‑1 (MQ=255)                                                                  cttcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACAt                                                                               >  5:74109/1‑149 (MQ=255)                                                                    tcttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTACCCCCCTGCAACATc                                                                              >  6:12536/1‑148 (MQ=255)                                                                     cttGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGt                                                                            >  4:293365/1‑149 (MQ=255)                                                                         tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTccac                                                                        >  7:325504/1‑149 (MQ=255)                                                                         tCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTcca                                                                         <  1:173724/148‑1 (MQ=255)                                                                                       aTTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACgg                                                          <  1:335020/149‑1 (MQ=255)                                                                                       aTTCCAGTAGAGGAAAAGGTATTCCTCATCTAATTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGTGGGGAATAGCTTGCCCGGAACTCGACCGAGTTGGGTTTCCACCCTTCATCATAGTCCACAACGCTTGGTACgt                                                          >  1:218178/1‑148 (MQ=255)                                                                                            aGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTggg                                                     <  4:224667/149‑1 (MQ=255)                                                                                             gTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTgggg                                                    <  2:114486/149‑1 (MQ=255)                                                                                              tAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTACA‑CACTCAGGGGGGAATAGGTTGACCGGAACTCGACCGTGTTGGGTTTCAACCCTTCAACATAGTCAACCAGGCTTGTTACGGTTggggg                                                   >  5:312149/1‑149 (MQ=255)                                                                                               agagGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGCTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCTCCATCGTCCACCCCGCTTGGTCCGGTTcgggg                                                   >  7:332346/1‑148 (MQ=255)                                                                                                gagGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGa                                                 >  2:87003/1‑149 (MQ=255)                                                                                                   gAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGAc                                                >  1:54661/1‑147 (MQ=255)                                                                                                   gAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACaa                                              <  4:122150/149‑1 (MQ=255)                                                                                                    aaaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACaa                                              <  6:100103/148‑1 (MQ=255)                                                                                                     aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGtt                                                                                                   >  8:324662/1‑94 (MQ=255)                                                                                                     aaaGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGtt                                                                                                   <  7:324662/94‑1 (MQ=255)                                                                                                              ccTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTAAGGTTGGTGGGACAAGAGAGGACGAg                                   >  4:238571/1‑149 (MQ=255)                                                                                                               cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGTGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGt                                  <  6:312149/149‑1 (MQ=255)                                                                                                               cTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGt                                  >  5:77202/1‑149 (MQ=255)                                                                                                                  aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCATCCCGCAGGGGGGAATAGGTTGCCCGGAACTCGAc                                                                                                         <  1:242741/76‑1 (MQ=255)                                                                                                                  aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAggggg                                                                                                                                <  4:15239/52‑1 (MQ=255)                                                                                                                  aTCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTc                                 <  1:87003/147‑1 (MQ=255)                                                                                                                    cTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGcc                                      >  6:165392/1‑140 (MQ=255)                                                                                                                    cTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGcc                                      <  5:165392/140‑1 (MQ=255)                                                                                                                         tCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAg                         <  3:293365/148‑1 (MQ=255)                                                                                                                            gggTGCGAGTTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCAACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGt                     >  5:359366/1‑149 (MQ=255)                                                                                                                                  gAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGt               >  5:180581/1‑149 (MQ=255)                                                                                                                                   aaTTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGtt              >  1:309729/1‑149 (MQ=255)                                                                                                                                        gCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTTGAGTTCTGc          <  6:359366/148‑1 (MQ=255)                                                                                                                                        gCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGc          <  6:77202/148‑1 (MQ=255)                                                                                                                                         cctgcGGGTGGAGTTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCgg        >  7:128105/5‑149 (MQ=255)                                                                                                                                         cctgcGGGTGGAGTTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCgg        >  4:343838/5‑149 (MQ=255)                                                                                                                                         cccACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg         <  4:290186/148‑1 (MQ=255)                                                                                                                                         cccACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCg         <  2:309729/148‑1 (MQ=255)                                                                                                                                          ccACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGt       >  1:138613/1‑149 (MQ=255)                                                                                                                                            aCGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCGAGTTCAGGGGTTGGAGTTCTGCGGTTa     >  5:237470/1‑149 (MQ=255)                                                                                                                                              gggTGGAGTTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTAcc   <  8:128105/149‑1 (MQ=255)                                                                                                                                                 tGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCg  >  7:260241/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                TCAATAAGCTTGTCGACCGGTACCAATATGAATAGCTCTTCAGGTTCGATAGGCTTGCTGAAGGTCTTCTTGTCGACCAGTACGTCATTCCAGTAGAGGAAAAGGTATTCCTCATCTATTTCCGGGTGCGAATTTGCCCACGGGTGGAACTCCA‑CCCGCAGGGGGGAATAGGTTGCCCGGAACTCGACCGTGTTGGGTTTCCCCCCTTCAACATCGTCCACCACGCTTGGTACGGTTGGGGGGACAAGAGAGGCCGAGTCCAGTTCAGGGGTTGGAGTTCTGCGGTTACCG  >  NC_002947/2104549‑2104839

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