Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	436,702	G→A	7.9%	G38G (GGG→GGA)	pgpA →	phosphatidylglycerophosphatase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	436,702	0	G	A	7.9%	190.1 / 5.0	63	G38G (GGG→GGA)	pgpA	phosphatidylglycerophosphatase A
Reads supporting (aligned to +/- strand):  ref base G (33/25);  new base A (2/3);  total (35/28)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

GTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATGG  >  NC_000913/436553‑436851                                                                                                                                                      |                                                                                                                                                       gTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTggg                                                                                                                                                        >  1:641665/1‑150 (MQ=255)  tGACGGCGAACCTGTTACATTAGACTGGAAAGGAAATGACCATTTTGCCACGCCATAAAGATTTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGAGGCATCTACTTGCTGTCGGATTCGGAACTGGATTAAGGCCGATCGTTCCTCGGa                                                                                                                                                       >  2:789979/1‑150 (MQ=255)          aaCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGAT‑cgggc                                                                                                                                             <  1:953/150‑3 (MQ=255)           aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCt                                                                                                                                             <  2:1365614/150‑1 (MQ=255)           aCCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGATCGGGCt                                                                                                                                             <  2:858171/150‑1 (MQ=255)            ccTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCC‑ATCGTCCCAGGAACGATCGGGCtt                                                                                                                                            <  1:650620/150‑2 (MQ=255)               gTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCt                                                                                                                                         <  2:670942/150‑1 (MQ=255)                 tACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTgg                                                                                                                                       >  2:723645/1‑150 (MQ=255)                  aCATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGc                                                                                                                                      >  2:567089/1‑150 (MQ=255)                    aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAg                                                                                                                                    <  2:14317/150‑1 (MQ=255)                    aTTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAg                                                                                                                                    <  2:1426574/150‑1 (MQ=255)                       aGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGa                                                                                                                                 >  1:205639/1‑150 (MQ=255)                          cTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTc                                                                                                                              >  1:1325229/1‑150 (MQ=255)                                  gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑g                                                                                                                                                 >  1:1012577/1‑123 (MQ=255)                                  gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑g                                                                                                                                                 <  2:1012577/123‑1 (MQ=255)                                  gATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTgg                                                                                                                      <  2:223072/150‑1 (MQ=255)                                      tGACCTTTTTGCGACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGc                                                                                                                                   <  1:955562/133‑1 (MQ=255)                                      tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGc                                                                                                                                   >  2:955562/1‑133 (MQ=255)                                      tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGt                                                                                                                           >  2:659949/1‑141 (MQ=255)                                      tGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGt                                                                                                                           <  1:659949/141‑1 (MQ=255)                                            ttttGCCACGCCATAATGATGTCGCGAAATGTCGCCTGAAGATGAGTAATCGGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCGCGTTCGTTCCTGGGTCGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCtgatga                                                                                                            <  2:991045/150‑1 (MQ=255)                                              ttGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGAcc                                                                                                          >  2:958991/1‑150 (MQ=255)                                                 ccACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCttt                                                                                                       >  1:510815/1‑150 (MQ=255)                                                  cACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCtttt                                                                                                      >  2:1396252/1‑150 (MQ=255)                                                     gCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGc                                                                                                   <  1:723645/150‑1 (MQ=255)                                                        aTAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                                >  2:461715/1‑150 (MQ=255)                                                        aTAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCt                                                                                                >  1:259525/1‑150 (MQ=255)                                                           aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGc                                                                                             >  1:650273/1‑150 (MQ=255)                                                           aaGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGGCGGATTCGGAAGTGGATTAAGCCCGAGCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGc                                                                                             <  1:1252892/150‑1 (MQ=255)                                                            aGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGGTCTGGTATCTGATGACCTTTTTGCCCTGGCa                                                                                            <  2:1289428/150‑1 (MQ=255)                                                                gTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGctc                                                                                        >  1:1087486/1‑150 (MQ=255)                                                                     gAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTc                                                                                   <  1:93280/150‑1 (MQ=255)                                                                          gTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtgg                                                                              >  1:953538/1‑150 (MQ=255)                                                                             gCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtgg                                                                           >  1:669892/1‑150 (MQ=255)                                                                             gCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCtggtgg                                                                           >  2:114157/1‑150 (MQ=255)                                                                              ccTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTATTTGCCCTGGCAGCTCTACTCGCTggtggt                                                                          >  1:1356610/1‑150 (MQ=255)                                                                                tGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGa                                                                        >  1:776567/1‑150 (MQ=255)                                                                                    gATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                                    <  2:510815/150‑1 (MQ=255)                                                                                    gATGAGTAATCCGTGGCATCTACTTGCTGTCAGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCt                                                                    <  1:1385744/150‑1 (MQ=255)                                                                                      tGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGAt                                                                       <  2:206344/145‑1 (MQ=255)                                                                                      tGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGAt                                                                       >  1:206344/1‑145 (MQ=255)                                                                                         gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGa                                                               >  1:115409/1‑150 (MQ=255)                                                                                         gTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGa                                                               >  1:787611/1‑150 (MQ=255)                                                                                           aaTCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATc                                                             <  1:114157/150‑1 (MQ=255)                                                                                             tCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTg                                                           <  2:865042/150‑1 (MQ=255)                                                                                             tCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTg                                                           <  1:567089/150‑1 (MQ=255)                                                                                                     aTCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCg                                                   <  1:1382370/150‑1 (MQ=255)                                                                                                        tACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCt                                                >  1:1144659/1‑150 (MQ=255)                                                                                                            tGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCt                                            >  1:274089/1‑150 (MQ=255)                                                                                                              cTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCttt                                          >  1:162978/1‑150 (MQ=255)                                                                                                                gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                                        >  1:350698/1‑150 (MQ=255)                                                                                                                gTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGt                                        <  1:958991/150‑1 (MQ=255)                                                                                                                          gAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTg                                                                                                    <  1:794575/80‑1 (MQ=255)                                                                                                                          gAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTg                                                                                                    >  2:794575/1‑80 (MQ=255)                                                                                                                            aGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCg                            >  2:156866/1‑150 (MQ=255)                                                                                                                               ggATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGaaa                         <  2:274089/150‑1 (MQ=255)                                                                                                                                gATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAg                        <  2:259525/150‑1 (MQ=255)                                                                                                                                      gCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg                  >  2:433100/1‑150 (MQ=255)                                                                                                                                       gcccaTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATgg                  >  1:555406/5‑150 (MQ=255)                                                                                                                                       cccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg                 <  2:1144659/150‑1 (MQ=255)                                                                                                                                       cccGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTAATCAAACGGCGAAAGACATggg                 >  1:454955/1‑150 (MQ=255)                                                                                                                                        cccaTCGTCCCAGGAACGATCGGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATggg                 >  2:363304/4‑150 (MQ=255)                                                                                                                                          gATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTTTCATCAAACGGCGAAAGACATGGgtgt              >  1:413306/1‑150 (MQ=255)                                                                                                                                          gATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGgtgt              >  1:188198/1‑150 (MQ=255)                                                                                                                                                    gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGa                           <  1:16084/127‑1 (MQ=255)                                                                                                                                                    gggACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGa                           >  2:16084/1‑127 (MQ=255)                                                                                                                                                      gACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATgg  <  2:162978/150‑1 (MQ=255)                                                                                                                                                      gACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATgg  <  1:1413827/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                       GTGACGGCGAACCTGTTACATTAGACTGGAAAGGATATGACCATTTTGCCACGCCATAAAGATGTCGCGAAAAGTCGCCTGAAGATGAGTAATCCGTGGCATCTACTTGCTGTCGGATTCGGAAGTGGATTAAGCCCGATCGTTCCTGGGACGAT‑GGGCTCGCTGGCAGCGATTCCGTTCTGGTATCTGATGACCTTTTTGCCCTGGCAGCTCTACTCGCTGGTGGTGATGCTGGGGATCTGTATCGGCGTCTATCTTTGTCATCAAACGGCGAAAGACATGGGTGTGCACGATCATGG  >  NC_000913/436553‑436851

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