Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,897,829	G→A	6.4%	intergenic (‑173/+61)	PP_5172 ← / ← PP_5173	hypothetical protein/RND family transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,897,829	0	G	A	6.4%	182.7 / 2.5	62	intergenic (‑173/+61)	PP_5172/PP_5173	hypothetical protein/RND family transporter
Reads supporting (aligned to +/- strand):  ref base G (33/25);  new base A (2/2);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

GCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATG  >  NC_002947/5897681‑5897977                                                                                                                                                     |                                                                                                                                                      gcCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggt                                                                                                                                                      <  1:1102285/150‑1 (MQ=255) gcCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggt                                                                                                                                                      >  2:1844025/1‑150 (MQ=255)  ccGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGTATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACggtg                                                                                                                                                     <  2:600421/150‑1 (MQ=255)   cGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGTCTTGAGGGCAATAGTAGTCAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTCGGAGCGGGTTTACCCGCGAATACggtgg                                                                                                                                                    <  2:996026/150‑1 (MQ=255)       cAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGa                                                                                                                                                <  2:1282656/150‑1 (MQ=255)            gTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCAc                                                                                                                                           >  1:1858695/1‑150 (MQ=255)                gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                       <  1:1142659/150‑1 (MQ=255)                gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                       >  1:770793/1‑150 (MQ=255)                gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                       <  2:1901717/150‑1 (MQ=255)                gAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATc                                                                                                                                       <  1:1683323/150‑1 (MQ=255)                  aGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGc                                                                                                                                     >  1:971393/1‑150 (MQ=255)                   ggggTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTGTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCa                                                                                                                                    >  2:996679/1‑150 (MQ=255)                        tttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTcgc                                                                                                                               <  1:1319288/150‑1 (MQ=255)                         ttGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTcgcg                                                                                                                              <  1:53218/150‑1 (MQ=25)                          tGCGTCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:1398314/150‑1 (MQ=255)                          tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:1608738/150‑1 (MQ=255)                          tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  1:923561/150‑1 (MQ=255)                          tGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATCCTAACCCTGAGTGTATAGTTATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  2:49692/150‑1 (MQ=255)                               cTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTaaa                                                                                                                        >  1:1253454/1‑150 (MQ=255)                                      gAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTc                                                                                                                 <  2:1858695/150‑1 (MQ=255)                                          gCAATAGTAGTGAGGCGGGATTAGTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCAGTCGCGGGTAAACCCGCtactac                                                                                                             >  2:950617/1‑150 (MQ=255)                                            aaTAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGGTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCa                                                                                                           >  1:287919/1‑150 (MQ=255)                                                 tagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGCATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTGTACCCGCGAATACGGTGGTCGATGCACCATCGCATTAGCGGGTAAAACCGCTACTACCAGGGG‑c                                                                                                     >  1:542051/1‑150 (MQ=255)                                                 tagtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGGGAACACGCCTCTGGTAGGCGCGGGCGTACCCGCGAATACGGCGGTCGATGCACCATCGCATTCGCGGGTAAACACGCTCGTACCAGGGG‑c                                                                                                     >  1:218849/1‑150 (MQ=38)                                                   gtGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCCCTGGTAGGAGCGGGTTTACCCGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGAGG‑cgt                                                                                                   <  1:313630/150‑3 (MQ=14)                                                    tGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTAGATGCACCATCGCATTCGCGGGTAAAACCGCTCCTAACAGGGG‑CCat                                                                                                  >  1:1851704/1‑150 (MQ=255)                                                          ggggTTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGc                                                                                            >  2:1540598/1‑150 (MQ=255)                                                             attattCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCtt                                                                                         <  1:757465/150‑1 (MQ=255)                                                               tattCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGACGTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGGCCCATATggcccctg                                                                                        <  1:1555404/150‑7 (MQ=39)                                                                attCTTACAAATACTGTTTAACAATTTTTATACAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCTAATACGGTGGTCGATGCCCCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAg                                                                                      <  1:878739/150‑1 (MQ=38)                                                                   cTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCg                                                                                   >  1:1440121/1‑150 (MQ=255)                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCgt                                                                                  >  1:1118941/1‑150 (MQ=255)                                                                    ttAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCgt                                                                                  >  2:574021/1‑150 (MQ=255)                                                                            aCTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGc                                                                          <  2:742841/150‑1 (MQ=255)                                                                               gTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGc                                                                                    <  2:1628832/137‑1 (MQ=255)                                                                               gTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGc                                                                                    >  1:1628832/1‑137 (MQ=255)                                                                                           ttATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTg                                                           <  1:1540598/150‑1 (MQ=255)                                                                                            tatataACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGt                                                          <  2:1253454/150‑1 (MQ=255)                                                                                                 aaCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAg                                                     <  1:996679/150‑1 (MQ=255)                                                                                                  aCCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAgg                                                    >  1:689892/1‑150 (MQ=255)                                                                                                        tAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGt                                              >  2:559913/1‑150 (MQ=255)                                                                                                         aaGAGAAAACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTa                                             >  1:1699648/1‑150 (MQ=255)                                                                                                           gagaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAgg                                                                                     <  1:1597465/108‑1 (MQ=255)                                                                                                           gagaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAgg                                                                                     >  2:1597465/1‑108 (MQ=255)                                                                                                           gagaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                           <  1:1822476/150‑1 (MQ=255)                                                                                                           gagaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACa                                           >  2:1671692/1‑150 (MQ=255)                                                                                                             gaACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAgc                                         >  2:1842512/1‑150 (MQ=255)                                                                                                                     ctctGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAg                                 <  1:1844025/150‑1 (MQ=255)                                                                                                                          gTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             >  1:649970/1‑54 (MQ=255)                                                                                                                          gTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCgg                                                                                                                             <  2:649970/54‑1 (MQ=255)                                                                                                                             ggAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCCGGCCACGTCCAGCGCCCGCAGGAAcagca                         <  2:1367863/150‑1 (MQ=255)                                                                                                                                    gTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAg                  >  1:1273346/1‑150 (MQ=255)                                                                                                                                     tttACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGc                 >  1:991789/1‑150 (MQ=255)                                                                                                                                       tACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCaa               <  2:689892/150‑1 (MQ=255)                                                                                                                                         cccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAgg             >  2:1142160/1‑150 (MQ=255)                                                                                                                                          ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTCCAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt            >  1:596458/1‑150 (MQ=25)                                                                                                                                          ccGCGAATGCGATGGTGCATCGACCACCGTATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt            >  1:1821871/1‑150 (MQ=25)                                                                                                                                          ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt            >  2:486203/1‑150 (MQ=255)                                                                                                                                          ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGt            >  1:669507/1‑150 (MQ=255)                                                                                                                                          ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGAAAAGt            >  2:1477250/1‑150 (MQ=255)                                                                                                                                          ccGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGCCCGCAAggg            >  1:858864/1‑149 (MQ=255)                                                                                                                                           cgcgAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGTCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTc           >  2:1128536/1‑150 (MQ=255)                                                                                                                                              gAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑c                                                                                                     <  1:43562/57‑1 (MQ=255)                                                                                                                                              gAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑c                                                                                                     >  2:43562/1‑57 (MQ=255)                                                                                                                                                aTACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCAc      >  1:535158/1‑150 (MQ=255)                                                                                                                                                 tACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCAGTGGCTCGCGGATCTTGTACCAGGCCACGTAAATTTTCGGCTGGAACAGCAGCGTCAGCAAGGGCGCCACg     >  1:1866496/1‑150 (MQ=255)                                                                                                                                                    ggtggtCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACgatg  >  1:683015/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCCGGGCAGAGGTTTGAAGGGGTTTTGCGCCTGGCTTGAGGGCAATAGTAGTGAGGCGGGATTATTCTTAAAAATACTGTTTAAGAATTTTTATATAACCATTTAAGAGAACACGCCTCTGGTAGGAGCGGGTTTACCCGCGAATACGGTGGTCGATGCACCATCGCATTCGCGGGTAAACCCGCTCCTACCAGGGG‑CCATATGGGCCTTCAGGCGTGTTTCTGCGGCTCGCGGATCTTGTACCAGGCCACGTACAGCGCCGGCAGGAACAGCAGCGTCAGCAAGGTCGCCACGATG  >  NC_002947/5897681‑5897977

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