Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	665,094	G→C	5.1%	P33A (CCT→GCT)	rlpA ←	septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	665,094	0	G	C	5.1%	325.4 / 6.0	97	P33A (CCT→GCT)	rlpA	septal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (51/40);  new base C (3/2);  total (54/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

AGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGAC  >  NC_000913/664946‑665243                                                                                                                                                      |                                                                                                                                                      aGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGATAATTTCAACTATAGGGCCGTTACATACCGCAGGc                                                                                                                                                      >  1:311833/1‑150 (MQ=255)     ccGCCTGGCTAACTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTAAAGGGCCGTTACATACCGCAGGCTGCg                                                                                                                                                  <  2:1568234/150‑1 (MQ=255)      cGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCgg                                                                                                                                                 >  2:1413327/1‑150 (MQ=255)      cGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCgg                                                                                                                                                 >  2:642719/1‑150 (MQ=255)       gCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGt                                                                                                                                                <  1:694939/150‑1 (MQ=255)          tGGCTAAATCGAGACGGATGCTGCACGATTTTGTAGCTTTTACCGTCGCGCGTGCAATCCTGATATGCCGTCGCGTGCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTatg                                                                                                                                             <  1:1235692/150‑3 (MQ=255)           ggCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTatgt                                                                                                                                            <  1:1052926/150‑4 (MQ=255)             cTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTAAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACtt                                                                                                                                          >  1:863845/1‑150 (MQ=255)             cTAAATCGAGACGGATCCTGCACGAATTTGTAGCTTTTAACGTCGCCCAGGTAATCCTGATATCCCGTCGCGTACAGTAGTTCGAAACGCGGGTCCG‑CCCCCCTAATTACAACTATAGGGCCGTAACATACCGCAGGCTGCGGTacacat                                                                                                                                          >  1:893886/1‑148 (MQ=255)                 aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CGCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACCCTTACCg                                                                                                                                      >  2:88176/1‑150 (MQ=255)                 aTCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCg                                                                                                                                      <  2:927570/150‑1 (MQ=255)                   cGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTc                                                                                                                                    >  2:439075/1‑150 (MQ=255)                   cGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTc                                                                                                                                    >  2:1405908/1‑150 (MQ=255)                          gATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGACGTCGCTTTCACTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCAGTTACATACCGCAGGCTGCGGTCCACTTACCGTCTGTTGCt                                                                                                                             >  1:70533/1‑150 (MQ=255)                             ccTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGAc                                                                                                                          <  2:1419343/150‑1 (MQ=255)                              cTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGAcc                                                                                                                         <  1:1676255/150‑1 (MQ=255)                                   cGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatcat                                                                                                                    >  2:237372/1‑150 (MQ=255)                                   cGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatcat                                                                                                                    <  1:1093925/150‑1 (MQ=255)                                   cGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACcatcaa                                                                                                                    >  2:451862/1‑149 (MQ=255)                                   cGATTTTGTAGCTTTAACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGTCTGCGGTACACTTACCGTCTGTTGCTGACcatcat                                                                                                                    <  1:235178/150‑1 (MQ=255)                                       tttGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCt                                                                                                                <  1:750976/150‑1 (MQ=255)                                        ttGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑GCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTATGTTTCTGACCATCATAGCtt                                                                                                               >  1:1493461/1‑150 (MQ=255)                                         tGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg                                                                                                              <  1:1462077/150‑1 (MQ=255)                                         tGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTg                                                                                                              >  2:145670/1‑150 (MQ=255)                                          gTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGt                                                                                                             <  1:642719/150‑1 (MQ=255)                                            aGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAc                                                                                                           >  1:1160528/1‑150 (MQ=255)                                             gCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACa                                                                                                          >  2:541415/1‑150 (MQ=255)                                              cTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTTATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACAt                                                                                                         >  1:1079165/1‑150 (MQ=255)                                              cTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTAAGTGGTTCGAAACGCGGGTCCG‑CGCCGCTAATTTAAACTATAGGGCCGTTACATACCGCAGGCTGCGGTAAACTTAACGTCTGATGCTGACCATCATCGCTTGttaat                                                                                                         >  1:1058052/1‑146 (MQ=255)                                                  tACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcg                                                                                                     <  2:1023608/150‑1 (MQ=255)                                                    ccGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgcg                                                                                                   <  2:779171/150‑1 (MQ=255)                                                    ccGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCcgcg                                                                                                   >  2:1266633/1‑150 (MQ=255)                                                       tCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGc                                                                                                >  1:444055/1‑150 (MQ=255)                                                         gcgcTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGCGCAt                                                                                              >  1:80710/1‑150 (MQ=255)                                                          cgcTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCAtt                                                                                             >  1:358145/1‑150 (MQ=255)                                                           gcTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTc                                                                                            >  2:378617/1‑150 (MQ=255)                                                           gcTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTc                                                                                            >  2:1161115/1‑150 (MQ=255)                                                            cTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTcc                                                                                           <  2:164048/150‑1 (MQ=255)                                                            cTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTcc                                                                                           >  2:1443915/1‑150 (MQ=255)                                                            cTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTcc                                                                                           <  2:1264504/150‑1 (MQ=255)                                                            cTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCAGTTGCAGACCATCCTCGCTTGTACATGCCGCGAGCATTcc                                                                                           >  2:1283353/1‑150 (MQ=255)                                                              ggTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTg                                                                                         <  1:221758/150‑1 (MQ=255)                                                              ggTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTg                                                                                         >  2:1223411/1‑150 (MQ=255)                                                                tAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGcc                                                                                       <  2:587575/150‑1 (MQ=255)                                                                  aTCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc                                                                                     >  1:362202/1‑150 (MQ=255)                                                                  aTCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc                                                                                     <  1:857359/150‑1 (MQ=255)                                                                  aTCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCcgc                                                                                     >  1:1496157/1‑150 (MQ=255)                                                                     cTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGAt                                                                                  >  1:42835/1‑150 (MQ=255)                                                                       gATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGc                                                                                >  2:588105/1‑150 (MQ=255)                                                                        aTTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCa                                                                               >  2:1497330/1‑150 (MQ=255)                                                                         tttGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAg                                                                              >  1:1236716/1‑150 (MQ=255)                                                                         tttGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAg                                                                              <  1:428410/150‑1 (MQ=255)                                                                           agacgacgcgTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATATCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTAACCGTCTGTTGCTGACCATCATCGCTTGTACAAGCCGCGAGCATTCCTGCCGCGATGCAGAt                                                                            <  1:1667824/144‑1 (MQ=255)                                                                             ccGTTGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATcc                                                                          <  1:614087/150‑1 (MQ=255)                                                                             ccGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATcc                                                                          >  1:404707/1‑150 (MQ=255)                                                                               gTCGCGTTCAGTGGTGCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCg                                                                        <  2:30923/150‑1 (MQ=255)                                                                                gcgcgTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCCGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGa                                                                       <  1:1255862/149‑1 (MQ=255)                                                                                 cgcgTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAg                                                                      >  2:1294596/1‑150 (MQ=255)                                                                                 cgcgTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAg                                                                      <  1:1087035/150‑1 (MQ=255)                                                                                  gcgTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGc                                                                     >  1:965528/1‑150 (MQ=255)                                                                                   cgTTCAGTGGTTCGAAACGCGGGTGCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGcc                                                                    >  1:317162/1‑150 (MQ=255)                                                                                     ttCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCAc                                                                  <  2:1236716/150‑1 (MQ=255)                                                                                       cAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTg                                                                >  2:599665/1‑150 (MQ=255)                                                                                            gTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTAc                                                           <  1:1112822/150‑1 (MQ=255)                                                                                             ttCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCGTCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACg                                                          >  1:423814/1‑150 (MQ=255)                                                                                                  aaCGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGc                                                                                <  1:1198488/123‑1 (MQ=255)                                                                                                  aaCGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGc                                                                                >  2:1198488/1‑123 (MQ=255)                                                                                                        ggTCCGCCCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAc                                                <  1:548436/150‑1 (MQ=255)                                                                                                        ggTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAcc                                               <  2:1366996/150‑1 (MQ=255)                                                                                                        ggTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAcc                                               <  1:1262153/150‑1 (MQ=255)                                                                                                        ggTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCAcc                                               <  1:1098608/150‑1 (MQ=255)                                                                                                         gTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACct                                              <  2:1106306/150‑1 (MQ=255)                                                                                                           ccG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACctct                                            >  1:1352704/1‑150 (MQ=255)                                                                                                             g‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTTCATGCCGCGAGCATTCCTGCATCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTa                                          <  1:88920/150‑1 (MQ=255)                                                                                                             g‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCATGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTAAGCATTGCGCACCTCTTa                                          >  2:1109739/1‑150 (MQ=255)                                                                                                               ccccGCTAATTTCAACTATAGGGCCGTTTCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTac                                         <  2:1693883/150‑1 (MQ=255)                                                                                                                cccGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTaca                                        <  1:588105/150‑1 (MQ=255)                                                                                                                    cTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCt                                    <  1:1413327/150‑1 (MQ=255)                                                                                                                      aaTTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCttt                                  <  1:1497330/150‑1 (MQ=255)                                                                                                                      aaTTTCAACTATAGGGCCGTTACATACCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGAGCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCTCCTCTTACACGCttt                                  >  1:1266784/1‑150 (MQ=255)                                                                                                                        tttCAACTATATGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGCCCATCATCGCTTGTACATTCCTCGAGCATTCCTGCCTCGATGCTGATCACGAGCCACTGCTTACGCATTGCGCACCTCATACACGCtatta                                >  2:1034671/1‑149 (MQ=255)                                                                                                                        tttCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTc                                <  2:584382/150‑1 (MQ=255)                                                                                                                         ttCAACTATAGGTCCGTTACTTCCCGCAGGCTGCGGTACACTTACCGTCGGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCAGCCCCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCg                               <  1:88176/150‑1 (MQ=255)                                                                                                                           cAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGac                             >  2:1378866/1‑150 (MQ=255)                                                                                                                            aaCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATTCAGATCCCGAGCCACTGCTTACGCTTTGAGCACCTCTTACACGCTTTTCGaaa                            >  2:1420473/1‑148 (MQ=255)                                                                                                                            aaCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGaca                            <  2:863845/150‑1 (MQ=255)                                                                                                                            aaCTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGaca                            >  2:1413912/1‑150 (MQ=255)                                                                                                                                ataGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAt                        <  1:281849/150‑1 (MQ=255)                                                                                                                                   gggCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACAtttt                     >  1:19003/1‑150 (MQ=255)                                                                                                                                      ccGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCt                  >  1:841318/1‑150 (MQ=255)                                                                                                                                       cGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCtg                 >  2:518570/1‑150 (MQ=255)                                                                                                                                           aCATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTAAATGCCGCGAGCATTCCTGCCGCGATGGAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTggg             <  1:1283353/150‑1 (MQ=255)                                                                                                                                            tgtACCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGAGGgt            >  1:411316/3‑150 (MQ=255)                                                                                                                                             gtACCGCAGCCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATACCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtg           >  2:647033/2‑150 (MQ=255)                                                                                                                                              tACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGgtgt          >  2:385385/1‑150 (MQ=255)                                                                                                                                                ccGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTgg        >  2:1305057/1‑150 (MQ=255)                                                                                                                                                  gCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGAt      >  2:1692139/1‑150 (MQ=255)                                                                                                                                                     ggCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAAAATTTTCCTGTGGGTGTGGATTGa   >  1:1382885/1‑150 (MQ=255)                                                                                                                                                     ggCTGCGGTACACTTACCGTCTGTTGCTGACCATCATAGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGATTACGCATTGCGCACATCTTACACGCTATTCGACAACATTTTCCTGTGGGTGTGGATTGa   >  1:911481/1‑150 (MQ=255)                                                                                                                                                      gCTGCGGTACACTTACCGTCTGTTGCTGACCATCCTCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGAc  <  2:423814/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                      AGTCCCGCCTGGCTAAATCGAGACGGATCCTGCACGATTTTGTAGCTTTTACCGTCGCGCTGGTAATCCTGATTTGCCGTCGCGTTCAGTGGTTCGAAACGCGGGTCCG‑CCCCGCTAATTTCAACTATAGGGCCGTTACATACCGCAGGCTGCGGTACACTTACCGTCTGTTGCTGACCATCATCGCTTGTACATGCCGCGAGCATTCCTGCCGCGATGCAGATCCCGAGCCACTGCTTACGCATTGCGCACCTCTTACACGCTTTTCGACAACATTTTCCTGTGGGTGTGGATTGAC  >  NC_000913/664946‑665243

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