Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,939,163	A→C	5.9%	intergenic (‑65/‑114)	PP_1738 ← / → PP_1739	hypothetical protein/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,939,163	0	A	C	5.9%	207.8 / 4.8	68	intergenic (‑65/‑114)	PP_1738/PP_1739	hypothetical protein/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base C (2/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TGGCAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCTGCCG  >  NC_002947/1939014‑1939306                                                                                                                                                      |                                                                                                                                                tGGCAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCaaa                                                                                                                                                 <  2:587336/150‑1 (MQ=255)   gCAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGt                                                                                                                                               >  1:450375/1‑150 (MQ=255)    cAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCGAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGtt                                                                                                                                              <  1:1216453/150‑1 (MQ=255)                tcgacgcTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGt                                                                                                                                  <  1:715570/146‑1 (MQ=255)                 ccGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTg                                                                                                                                 >  1:1385915/1‑150 (MQ=255)                   gcgcTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGc                                                                                                                               <  2:376434/150‑1 (MQ=255)                         aaGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTgg                                                                                                                         >  1:972591/1‑150 (MQ=255)                         aaGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTgg                                                                                                                         <  2:527967/150‑1 (MQ=255)                                cacaGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACtt                                                                                                                  <  2:450375/150‑1 (MQ=255)                                  caGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGc                                                                                                                <  2:1385915/150‑1 (MQ=255)                                       aTGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCt                                                                                                           <  2:972591/150‑1 (MQ=255)                                          aGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCa                                                                                                        >  2:1493531/1‑150 (MQ=255)                                                 tGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAACTTTGCCCGcctggt                                                                                                                                   >  2:1608924/1‑111 (MQ=255)                                                 tGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAACTTTGCCCGcctggt                                                                                                                                   <  1:1608924/116‑6 (MQ=255)                                                 tGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGc                                                                                                 >  1:830820/1‑150 (MQ=255)                                                       gCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCgg                                                                                           <  2:208816/150‑1 (MQ=255)                                                                 gAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATGTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGAttt                                                                                 >  2:1418658/1‑150 (MQ=255)                                                                    gcgcAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGc                                                                              >  2:1437190/1‑150 (MQ=255)                                                                        aaaGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTa                                                                          >  1:426491/1‑150 (MQ=255)                                                                            aTCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGaaa                                                                      >  1:708331/1‑150 (MQ=255)                                                                            aTCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGaaa                                                                      >  2:339558/1‑150 (MQ=255)                                                                             tCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAAc                                                                     >  2:701776/1‑150 (MQ=255)                                                                               gAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   <  2:139255/150‑1 (MQ=255)                                                                                   cATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCAt                                                               <  2:1554357/150‑1 (MQ=255)                                                                                   cATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCA‑GCGGGCAACTTTGCCCGcctg                                                                                                                                     >  2:1250995/1‑76 (MQ=37)                                                                                   cATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCA‑GCGGGCAACTTTGCCCGcctg                                                                                                                                     <  1:1250995/79‑4 (MQ=37)                                                                                       tAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTACGATATGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGc                                                           >  1:515605/1‑150 (MQ=255)                                                                                        aaTCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTg                                                                                                                          <  1:245681/86‑1 (MQ=255)                                                                                        aaTCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTg                                                                                                                          >  2:245681/1‑86 (MQ=255)                                                                                          tCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGAt                                                        <  2:426491/150‑1 (MQ=255)                                                                                            ccccAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTAc                                                                                                                    >  1:1203252/1‑88 (MQ=255)                                                                                            ccccAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTAc                                                                                                                    <  2:1203252/88‑1 (MQ=255)                                                                                             cccAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCAt                                                     >  2:1558260/1‑150 (MQ=255)                                                                                                aaaaGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCAt                                                                                                                             >  1:647796/1‑75 (MQ=255)                                                                                                aaaaGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCAt                                                                                                                             <  2:647796/75‑1 (MQ=255)                                                                                                 aaaGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGAt                                                 <  1:1493531/150‑1 (MQ=255)                                                                                                      aCGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCttt                                            <  1:701776/150‑1 (MQ=255)                                                                                                        gCAAGGTCGAGCAGTGTGCAGGCTATGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCg                                          >  1:703960/1‑150 (MQ=255)                                                                                                         cagggTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGa                                         >  1:1643078/4‑150 (MQ=255)                                                                                                          aaGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAATTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGAt                                        >  2:1028279/1‑150 (MQ=255)                                                                                                          aaGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGAt                                        <  2:1643078/150‑1 (MQ=255)                                                                                                               cGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAATTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCg                                   <  1:1028279/150‑1 (MQ=255)                                                                                                                   cAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTAt                               <  1:1437190/150‑1 (MQ=255)                                                                                                                    aGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTAtt                              <  1:1418658/150‑1 (MQ=255)                                                                                                                     gtgtGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTg                             >  1:202184/1‑150 (MQ=255)                                                                                                                          cAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGaa                        <  1:1558260/150‑1 (MQ=255)                                                                                                                          cAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGaa                        >  2:527170/1‑150 (MQ=255)                                                                                                                           aGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   <  1:1467416/106‑1 (MQ=255)                                                                                                                           aGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   >  2:1467416/1‑106 (MQ=255)                                                                                                                           aGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAAc                       >  1:536704/1‑150 (MQ=255)                                                                                                                             gCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   <  1:1487268/104‑1 (MQ=255)                                                                                                                             gCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACgg                                                                   >  2:1487268/1‑104 (MQ=255)                                                                                                                             gCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAg                     >  1:152958/1‑150 (MQ=255)                                                                                                                              cTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    >  1:651479/1‑150 (MQ=255)                                                                                                                               tAGGGCCACAACCAGGCGGGGAAAGTTGCCTGCTGGGTTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGa                   <  2:158018/150‑1 (MQ=255)                                                                                                                                aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATc                                                                                                            >  1:1419027/1‑60 (MQ=255)                                                                                                                                aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATc                                                                                                            <  2:1419027/60‑1 (MQ=255)                                                                                                                                aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGgtgt                                                                                      >  2:174628/1‑82 (MQ=255)                                                                                                                                aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGgtgt                                                                                      <  1:174628/82‑1 (MQ=255)                                                                                                                                aGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGAt                  <  2:515605/150‑1 (MQ=255)                                                                                                                                  gACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGt                >  2:1144185/1‑150 (MQ=255)                                                                                                                                     cacaACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACGTAGCTGCCGGGGTGTGAtttt                                                                                <  1:1147538/83‑1 (MQ=255)                                                                                                                                     cacaACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGAtttt                                                                                >  2:1147538/1‑83 (MQ=255)                                                                                                                                     cacaACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGtctc             >  2:787524/1‑150 (MQ=255)                                                                                                                                         aCCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTc         >  2:628331/1‑150 (MQ=255)                                                                                                                                            aGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCt      <  1:527170/150‑1 (MQ=255)                                                                                                                                                gggCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAAAAGGGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCTGCCg  >  1:40771/1‑150 (MQ=255)                                                                                                                                                 ggCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    <  1:599220/131‑1 (MQ=255)                                                                                                                                                 ggCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    >  2:599220/1‑131 (MQ=255)                                                                                                                                                    aaaGTTGCCCGCTGCGGTGGCATTTGACGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    <  1:1218742/128‑1 (MQ=255)                                                                                                                                                    aaaGTTGCCCGCTGCGGTGGCATTTGACGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAgg                    >  2:1218742/1‑128 (MQ=255)                                                                                                                                                      |                                                                                                                                                TGGCAGCGTCGATATCCCGCGCTGAAGCGGCCACAGCCATGAGCAAGCTGCCGAGCAGGGCAGCGAGGCGCAAAGATCGAAGCATGTAATCCCCCAAAAGAACGCAAGGTCGAGCAGTGTGCAGGCTAGGACCACAACCAGGCGGGCAAAGTTGCCCGCTGCGGTGGCATTTGGCGTACTTGCTATCTGCAACATAGCTGCCGGGGTGTGATTTTGCCTTAGAAACGGGCATAAGCGATCATGGATGCTTTCGATTACCGCTATTGGAGAACAGGATGTCTCGTTCGCTGCCG  >  NC_002947/1939014‑1939306

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