Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,264,059	C→G	8.3%	intergenic (+188/‑35)	trpF → / → accD	N‑(5'‑phosphoribosyl)anthranilate isomerase/acetyl‑CoA carboxylase carboxyltransferase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,264,059	0	C	G	8.3%	173.6 / 5.8	60	intergenic (+188/‑35)	trpF/accD	N‑(5'‑phosphoribosyl)anthranilate isomerase/acetyl‑CoA carboxylase carboxyltransferase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (24/30);  new base G (2/3);  total (26/34)
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.19e-01

CTTTCGCAGCCCTGCGGGCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAGTGCCCGTCCTGCGAGGCCGTGCTGT  >  NC_002947/2263911‑2264208                                                                                                                                                          |                                                                                                                                                       cTTTCGCAGCCCTGCGGGCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGcc                                                                                                                                                       <  1:919681/150‑1 (MQ=255) cTTTCGCAGCCCTGCGGGCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGcc                                                                                                                                                       >  2:479375/1‑150 (MQ=255)   cacggcctcgcaggacgggcacttgtgccacaggccctcaggcaccgagctctttttcacctcggaacgcatgatcgaagggatcagtttgtcgactaaccagttgctcatgctttctttctccagtattGGCAGGGC‑‑‑‑‑CGGTGACGCCGt                                                                                                                                                     <  2:259953‑M1/20‑1 (MQ=255)    tCGCAGCCCTGCGGGCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTc                                                                                                                                                    <  1:387903/150‑1 (MQ=255)      gCAGCCCTGCGGGCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTccc                                                                                                                                                  >  1:843650/1‑150 (MQ=255)              aggacgggcacttgtgccacaggccctcaggcaccgagctctttttcacctcggaacgcatgatcgaagggatcagtttgtcgactaaccagttgctcatgctttctttctccagtattGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGc                                                                                                                                         >  1:259953‑M1/120‑150 (MQ=255)                  gCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑Tg                                                                                                                                          <  1:1105387/145‑1 (MQ=255)                  gCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑Tg                                                                                                                                          >  2:1105387/1‑145 (MQ=255)                      cATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTg                                                                                                                                 <  1:631410/150‑1 (MQ=255)                        tGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGga                                                                                                                               >  1:521386/1‑150 (MQ=255)                        tGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGga                                                                                                                               <  1:657541/150‑1 (MQ=255)                           gtgccacaggccctcaggcaccgagctctttttcacctcggaacgcatgatcgaagggatcagtttgtcgactaaccagttgctcatgctttctttctccagtattGGCAGGGCAGGGCCGGTGACGCCGGCCCCCGCC‑TGCCAATACTg                                                                                                                                 >  2:1006748‑M1/107‑150 (MQ=255)                           cAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGagaa                                                                                                                            >  2:1124192/1‑150 (MQ=255)                           cAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGagaa                                                                                                                            >  1:297791/1‑150 (MQ=255)                            aGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGagaaa                                                                                                                           >  1:491415/1‑150 (MQ=255)                                   gtgGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCg                                                                                                                                                      >  1:355314/1‑117 (MQ=255)                                   gtgGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCg                                                                                                                                                      <  2:355314/117‑1 (MQ=255)                                    tgGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCAt                                                                                                                   >  1:19480/1‑150 (MQ=255)                                      gcggcAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGa                                                                                                                 <  1:168879/150‑1 (MQ=255)                                         gcAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCa                                                                                                              >  1:508140/1‑150 (MQ=255)                                          cAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCaa                                                                                                             <  2:843650/150‑1 (MQ=255)                                           aaaaaaaCACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAAc                                                                                                            >  2:131334/1‑150 (MQ=255)                                                 acacATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGa                                                                                                                 <  1:261255/139‑1 (MQ=255)                                                 acacATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGa                                                                                                                 >  2:261255/1‑139 (MQ=255)                                                    caTTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTc                                                                                                   <  1:742492/150‑1 (MQ=255)                                                             tGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAg                                                                                                     >  1:925593/1‑139 (MQ=255)                                                             tGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAg                                                                                                     <  2:925593/139‑1 (MQ=255)                                                               gAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGAt                                                                                        <  1:509473/150‑1 (MQ=255)                                                                 gTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATcc                                                                                      <  2:882314/150‑1 (MQ=255)                                                                    gcgcgcTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGGTAGTCTACAAACTGATCCCtt                                                                                   <  2:457609/150‑1 (MQ=255)                                                                         cTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATc                                                                              <  2:70878/150‑1 (MQ=255)                                                                           ggAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGc                                                                                                                                                        >  2:297347/1‑75 (MQ=255)                                                                           ggAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGc                                                                                                                                                        <  1:297347/75‑1 (MQ=255)                                                                           ggAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCAt                                                                            <  2:297791/150‑1 (MQ=255)                                                                            gAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATg                                                                           >  1:829558/1‑150 (MQ=255)                                                                            gAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATg                                                                           >  1:414556/1‑150 (MQ=255)                                                                                  ccGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAAGGC‑‑‑‑‑CGGTTACTTCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTcc                                                                     <  1:131334/150‑1 (MQ=255)                                                                                       ggcCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCGTCGATCATGCGTTCCGAGGt                                                                <  2:521386/150‑1 (MQ=255)                                                                                          cGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCa                                                                                                                    >  1:892194/1‑95 (MQ=255)                                                                                          cGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCa                                                                                                                    <  2:892194/95‑1 (MQ=255)                                                                                               caccaTCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAaga                                                        <  2:885454/150‑1 (MQ=255)                                                                                                   taaccagttgctcatgctttctttctccagtattGGCAGG‑C‑‑‑‑‑GGGGGACGGCGTCACCGGCCCTGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTc                                                    <  2:998250‑M1/116‑1 (MQ=255)                                                                                                   aTCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTc                                                    >  1:1057473/1‑150 (MQ=255)                                                                                                      ttttCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGt                                                 <  2:886216/149‑1 (MQ=255)                                                                                                       tttCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGtt                                                                       <  1:56879/127‑1 (MQ=255)                                                                                                       tttCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGtt                                                                       >  2:56879/1‑127 (MQ=255)                                                                                                       tttCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACAGCGTTGCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTg                                                >  1:749240/1‑150 (MQ=255)                                                                                                          cATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCt                                             <  1:271167/150‑1 (MQ=255)                                                                                                               aaGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAggg                                        <  2:1021441/150‑1 (MQ=255)                                                                                                                 tgctttctttctccagtattGGCAGG‑C‑‑‑‑‑GGGGGACGGCGTCACCGGCCCTGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGcc                                      >  1:998250‑M1/21‑150 (MQ=255)                                                                                                                  aTTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCt                                     <  1:479375/150‑1 (MQ=255)                                                                                                                     ttGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCtgtg                                  <  2:414556/150‑1 (MQ=255)                                                                                                                           tCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTc                                                                                                   <  1:782259/79‑1 (MQ=255)                                                                                                                           tCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTc                                                                                                   >  2:782259/1‑79 (MQ=255)                                                                                                                           tCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACaa                            >  1:799165/1‑150 (MQ=255)                                                                                                                           tCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACaa                            <  2:491415/150‑1 (MQ=255)                                                                                                                           tCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACaa                            >  1:794753/1‑150 (MQ=255)                                                                                                                           gcAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCAAGAAAAAAAAAAAAAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACaa                            <  2:709570/149‑1 (MQ=255)                                                                                                                            tccagtattGGCAGG‑C‑‑‑‑‑GGGGGACGGCGTCACCGGCCCTGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAg                           <  1:196306‑M1/141‑1 (MQ=255)                                                                                                                             ccagtattGGCAGG‑C‑‑‑‑‑GGGGGACGGCGTCACCGGCCCTGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAGt                          <  1:697262‑M1/142‑1 (MQ=255)                                                                                                                               gggcagggcagggc‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAGTGc                        <  2:749240/150‑1 (MQ=255)                                                                                                                                 gcagggcagggc‑‑‑‑‑CGGTGACGCCGTCCCACGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGGGAAACAGAGCTCGGTGCCTGGGGGCCTGTGGCACACGTGccc                      >  1:367833/1‑150 (MQ=255)                                                                                                                                      GCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAGTGCCCGTCCt                 <  1:679776‑M1/150‑1 (MQ=255)                                                                                                                                        AGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAGTGCCCGTCCTGc               >  1:141362‑M1/1‑150 (MQ=255)                                                                                                                                                          ccGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAGTGCCCGTCCTGCGAGGCCGTGCTGt  <  1:1124192/150‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                       CTTTCGCAGCCCTGCGGGCCGCATGCCAGCACCTGTGGCGGCAAAAAAACACATTAGCGGTGGAGTGGCGCGCTGGAAGCCCCGGCGGCCGGTCCACCATCGTTTCATAAAAGATTTTGAGCTCAAGGGCAGGGCAGGGC‑‑‑‑‑CGGTGACGCCGTCCCCCGCC‑TGCCAATACTGGAGAAAGAAAGCATGAGCAACTGGTTAGTCGACAAACTGATCCCTTCGATCATGCGTTCCGAGGTGAAAAAGAGCTCGGTGCCTGAGGGCCTGTGGCACAAGTGCCCGTCCTGCGAGGCCGTGCTGT  >  NC_002947/2263911‑2264208

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