Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,851,737	T→A	7.7%	intergenic (+84/+35)	apt → / ← recR	adenine phosphoribosyltransferase/recombination protein RecR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,851,737	0	T	A	7.7%	189.1 / 7.3	65	intergenic (+84/+35)	apt/recR	adenine phosphoribosyltransferase/recombination protein RecR
Reads supporting (aligned to +/- strand):  ref base T (31/29);  new base A (3/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.89e-01

GCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGC  >  NC_002947/4851600‑4851884                                                                                                                                          |                                                                                                                                                    gCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATgg                                                                                                                                         <  2:589222/150‑1 (MQ=255)           ggCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               >  2:924524/2‑150 (MQ=255)                aCCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTAGTGGGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGccc                                                                                                                          <  2:568989/150‑1 (MQ=255)                   ttCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCaa                                                                                                                       <  2:488196/150‑1 (MQ=255)                    tCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAACCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAAt                                                                                                                      <  1:137536/150‑1 (MQ=255)                    tCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAAt                                                                                                                      <  1:437404/150‑1 (MQ=255)                      tGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                    <  1:924524/150‑1 (MQ=255)                          tGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCa                                                                                                                <  1:442527/150‑1 (MQ=255)                              cGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTCCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCGGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGc                                                                                                            >  1:581227/1‑150 (MQ=255)                                     tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATgg                                                                                                     >  2:50816/1‑150 (MQ=255)                                     tCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATgg                                                                                                     <  1:298550/150‑1 (MQ=255)                                       ccTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCCGAGCGAAATGGGc                                                                                                   >  2:129056/1‑150 (MQ=255)                                       ccTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGGAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGc                                                                                                   >  1:613655/1‑150 (MQ=255)                                        cTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGGATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGcc                                                                                                  >  1:720564/1‑150 (MQ=255)                                         tCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGTCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGAAGCTTAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                                                                                 <  1:699317/150‑1 (MQ=39)                                            gCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                             <  1:56168/103‑2 (MQ=255)                                            gCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTAAACCcgcgc                                                                                                                                             >  2:56168/1‑102 (MQ=255)                                              gAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGcc                                                                                            >  1:1027285/1‑150 (MQ=255)                                                 tATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGGAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATGAGAGCGAAATGGGCCGGCGCCCTg                                                                                         >  1:193868/1‑150 (MQ=255)                                                  aTTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCaga                                                                                                              >  1:294629/1‑128 (MQ=255)                                                  aTTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCaga                                                                                                              <  2:294629/128‑1 (MQ=255)                                                       gTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                    >  2:856854/1‑117 (MQ=255)                                                       gTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGa                                                                                                                    <  1:856854/117‑1 (MQ=255)                                                        tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGgcc                                                                                                                                      >  2:411957/1‑98 (MQ=255)                                                        tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGgcc                                                                                                                                      <  1:411957/98‑1 (MQ=255)                                                        tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               >  1:105309/1‑105 (MQ=255)                                                        tGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               <  2:105309/105‑1 (MQ=255)                                                          ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               >  2:579222/1‑103 (MQ=255)                                                          ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGc                                                                                                                               <  1:579222/103‑1 (MQ=255)                                                          ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCAt                                                                                >  2:680839/1‑150 (MQ=255)                                                          ggCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCAt                                                                                <  2:1027285/150‑1 (MQ=255)                                                           gCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATg                                                                               <  1:1135200/150‑1 (MQ=255)                                                            tttttGGAGCTGCGTTGCAGCCCATTCGCGGGTAACCCCAGTCGTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATgg                                                                              <  1:616286/149‑1 (MQ=39)                                                            cTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATgg                                                                              <  2:505860/150‑1 (MQ=255)                                                            cTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATgg                                                                              <  1:904234/150‑1 (MQ=255)                                                             ttttGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATggg                                                                             >  1:830052/1‑150 (MQ=255)                                                                   ggCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAggg                                                                       >  1:312932/1‑150 (MQ=255)                                                                         tttGCAGCCCATTAGATTATAAACCCGCTCCTACGAAAGCCGCGTATGTGTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATTAGAGAGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTAccgc                                                                 >  1:947769/1‑150 (MQ=38)                                                                              aGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcc                                                            <  2:520801/150‑1 (MQ=255)                                                                               gcgcATTCGCGGCTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGAGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGACCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                           >  2:550370/4‑150 (MQ=255)                                                                               gCCCATTCGCGGGTAACCCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                           <  1:1049899/150‑1 (MQ=255)                                                                               gCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTcca                                                           >  1:507827/1‑150 (MQ=255)                                                                                cccATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTGTGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTccac                                                          >  2:697409/1‑150 (MQ=255)                                                                                       gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                                                                                 <  2:1142262/104‑1 (MQ=255)                                                                                       gcgGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCg                                                                                                 >  1:1142262/1‑104 (MQ=255)                                                                                           gTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCt                                               >  2:1059901/1‑150 (MQ=255)                                                                                           gTAAACCCGCTCCTACGAAAGCCGCGTATGATTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCt                                               >  2:536108/1‑150 (MQ=255)                                                                                                cccGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgccg                                          <  2:610179/150‑1 (MQ=255)                                                                                                 ccGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTcgccgc                                         >  1:516623/1‑150 (MQ=255)                                                                                                        tACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg                                  <  2:613655/150‑1 (MQ=255)                                                                                                        tACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTgg                                  >  2:821507/1‑150 (MQ=255)                                                                                                         aCGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGc                                 >  1:1017640/1‑150 (MQ=255)                                                                                                           gAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGcac                               >  1:728903/1‑150 (MQ=255)                                                                                                                         ttttGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGc                 >  2:17991/1‑150 (MQ=255)                                                                                                                           ttGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCga               <  2:312932/150‑1 (MQ=255)                                                                                                                           ttGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCga               <  1:821507/150‑1 (MQ=255)                                                                                                                             gTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga             <  1:284858/150‑1 (MQ=255)                                                                                                                             ccattcgcgGGTAAACCCGCGCATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCaa                    <  2:438102/138‑1 (MQ=255)                                                                                                                             ccattcgcgGGTAAACCCGCGCATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCaa                    >  1:438102/6‑143 (MQ=255)                                                                                                                                ttcgcgGGTAAACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGc          >  1:954318/3‑150 (MQ=255)                                                                                                                                ggCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGc          >  2:389392/1‑150 (MQ=255)                                                                                                                                      ggTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCaa    <  2:1017640/150‑1 (MQ=255)                                                                                                                                       gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGGGGCCTAGTGGCACACCATGGGCAATGCGAGACGCCACCAAg   >  1:968310/1‑150 (MQ=255)                                                                                                                                       gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga             <  2:249535/140‑1 (MQ=255)                                                                                                                                       gTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCgaga             >  1:249535/1‑140 (MQ=255)                                                                                                                                        tttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTg                                   >  1:1160833/1‑117 (MQ=255)                                                                                                                                        tttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTg                                   <  2:1160833/117‑1 (MQ=255)                                                                                                                                        tttACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGc  <  1:17991/150‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                    GCGGCCGGGGTGCCTACCTTCTGCCTGACCGAGTTTTCCCTCAGCGAGTATTGAGTGGGCTTTTGGGGCTGCTTTGCAGCCCATTCGCGGGTAAACCCGCTCCTACGAAAGCCGCGTATGTTTTGTAGGCGCGGGTTTACCCGCGAATGGGCCGCCAAGCGGCCCCAATGAATCAGAGCGAAATGGGCCGGCGCCCTGCAAAGGCATGGGCCAGGGTACCGCCGTCCACCATCTCCAGCTCGCCGCCCAGTGGCACACCATGGGCAATGCGAGACGCCACCAAGC  >  NC_002947/4851600‑4851884

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