Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,788,441	G→C	7.6%	Y79* (TAC→TAG)	PP_3350 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,788,441	0	G	C	7.6%	182.6 / 5.2	62	Y79* (TAC→TAG)	PP_3350	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/27);  new base C (3/2);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.70e-01

CGGGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGC  >  NC_002947/3788298‑3788585                                                                                                                                                  |                                                                                                                                                 cGGGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCgg                                                                                                                                            <  2:3547/150‑1 (MQ=255)    gTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGCTACCGGGcc                                                                                                                                         >  2:541386/1‑150 (MQ=255)     tATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGCTACCGGGccg                                                                                                                                        <  1:541386/150‑2 (MQ=255)        tctcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAg                                                                                                                                     <  2:1277229/150‑1 (MQ=255)         ctcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGc                                                                                                                                    >  1:1031727/1‑150 (MQ=255)          tcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGTGT‑GACATATCAGCCCGCGCATGCCAAAGAAGGGTTTGTCCGAAGCAT‑CCGCCTGGTAGCGGGCCTCAGcc                                                                                                                                   >  1:798474/1‑150 (MQ=255)          tcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGcc                                                                                                                                   >  2:717922/1‑150 (MQ=255)          tcGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGcc                                                                                                                                   >  2:12145/1‑150 (MQ=255)             aCCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGCTACCGGGccgat                                                                                                                                      <  2:436344/144‑4 (MQ=255)             aCCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGCTACCGGGccgat                                                                                                                                      >  1:436344/1‑141 (MQ=255)               cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAg                                                                                                                              <  1:1051528/150‑1 (MQ=255)               cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAg                                                                                                                              >  1:704903/1‑150 (MQ=255)               cTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAg                                                                                                                              <  1:403677/150‑1 (MQ=255)                tGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGt                                                                                                                             <  1:12145/150‑1 (MQ=255)                         gTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCATCAGCCCGGAGAAGCAGccgg                                                                                                                    >  2:80356/1‑148 (MQ=255)                             gCACGGCATCGTCGACCAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCTCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGCCCCGGTAGCGGCCCTCAGCCCCGAGTTGCACCCAGGCCg                                                                                                                >  2:128038/1‑150 (MQ=255)                                 ggCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGActc                                                                                                            <  1:634412/150‑1 (MQ=255)                                           aaCAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGa                                                                                                  >  2:861854/1‑150 (MQ=255)                                             cAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCATACCGTCTTCAAGC‑TGCTTCTGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATg                                                                                                >  2:1106720/1‑150 (MQ=255)                                              aGGTGCAGGTCGGCATGGGCCTGCAGGCGTTGCTGCAGCTAGGAATCGTCGT‑TCAGCCCGCGCATGCGAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATgg                                                                                               >  2:376334/1‑150 (MQ=255)                                              aGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATgg                                                                                               <  1:868872/150‑1 (MQ=255)                                                     ggTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCa                                                                                        >  2:1076434/1‑150 (MQ=255)                                                      gTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAg                                                                                       >  2:754604/1‑150 (MQ=255)                                                      gTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGACGT‑TCAGCCCGCGAAGGCCAAAGAAGGGTTGGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAg                                                                                       >  1:532670/1‑150 (MQ=255)                                                           aaTGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCAt                                                                                  <  2:651803/149‑1 (MQ=255)                                                               ggCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGaa                                                                              >  1:724695/1‑150 (MQ=255)                                                                gCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGcccatgccg                                                                             >  2:1292529/1‑143 (MQ=255)                                                                     cAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGCAAGGGGc                                                                        >  2:399935/1‑150 (MQ=255)                                                                      aGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGa                                                                                                                               >  2:489370/1‑94 (MQ=255)                                                                      aGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGa                                                                                                                               <  1:489370/94‑1 (MQ=255)                                                                        ccGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCCGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGt                                                                     <  2:821014/150‑1 (MQ=255)                                                                         aacggctgcAGATAGGAATCGCAGC‑TCATCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTg                                                                    <  2:331781/146‑1 (MQ=255)                                                                          gTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTgg                                                                   >  1:446690/1‑150 (MQ=255)                                                                              ctgcAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCtt                                                               <  1:156111/150‑1 (MQ=255)                                                                                gcAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGc                                                                                                                                          <  1:1174197/73‑1 (MQ=255)                                                                                gcAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGc                                                                                                                                          >  2:1174197/1‑73 (MQ=255)                                                                                   gATAGGAATCGTCGT‑TCAGCCTGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGTGGGTCTGCGGCCCGCGTTGCAGCCAGGCCGACTATGACATTCGATAGTGGCGCAGGGCATTGAGAGTGGCAGTGGCCTTGGCgg                                                          >  2:86745/1‑150 (MQ=255)                                                                                     tAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCgggg                                                        >  1:566660/1‑150 (MQ=255)                                                                                        gAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcg                                                     >  1:961444/1‑150 (MQ=255)                                                                                        gAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACCCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcg                                                     >  2:379089/1‑150 (MQ=255)                                                                                            cgtcgt‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCga                                                 <  1:168515/150‑1 (MQ=255)                                                                                                     aGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGGCAGCCGATGGTAGCGCAGGCCATCGAAAGAGATCACGCGCTGGGCGGGGTCGTCGAGAAAGCgg                                         <  1:1106720/150‑1 (MQ=255)                                                                                                     aGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCgg                                         <  1:273840/150‑1 (MQ=255)                                                                                                       cccGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTACGCGGGGTCGTCGAGAAAGCGGTa                                       <  2:798474/150‑1 (MQ=255)                                                                                                          gcgcAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTc                                    >  2:1406746/1‑150 (MQ=255)                                                                                                          gcgcAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTc                                    >  2:22697/1‑150 (MQ=255)                                                                                                             cAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTc                                 <  1:1406746/150‑1 (MQ=255)                                                                                                                  cAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCggg                                                         >  1:689382/1‑121 (MQ=255)                                                                                                                  cAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCggg                                                         <  2:689382/121‑1 (MQ=255)                                                                                                                   aaaGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCa                           <  1:717922/150‑1 (MQ=255)                                                                                                                        agGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGACCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAg                      <  1:80356/150‑1 (MQ=255)                                                                                                                            tttGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGTCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGc                  <  1:86745/150‑1 (MQ=255)                                                                                                                            tttGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGAGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                   <  1:425571/117‑1 (MQ=255)                                                                                                                            tttGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGAGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGgtcgtc                                                   >  2:425571/1‑117 (MQ=255)                                                                                                                               gTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGcccc               <  2:446690/150‑1 (MQ=255)                                                                                                                                tCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGa                                                                                                  <  1:912828/66‑1 (MQ=255)                                                                                                                                tCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGa                                                                                                  >  2:912828/1‑66 (MQ=255)                                                                                                                                   gAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCCGGc           >  2:1337032/1‑150 (MQ=255)                                                                                                                                    aaGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGCAGCGCCGGGcc          >  1:770814/1‑150 (MQ=255)                                                                                                                                            ggCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCg                                          <  2:82467/111‑1 (MQ=255)                                                                                                                                            ggCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCg                                          >  1:82467/1‑111 (MQ=255)                                                                                                                                             gCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCccggc  <  1:376334/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                 CGGGTATTCTCGACCTGGACGAAGGTACGCACGGCATCGTCGAACAGGTGCAGGTCGGCATGGGCCTGCAGCCGTTGCTGCAGATAGGAATCGTCGT‑TCAGCCCGCGCAGGCCAAAGAAGGGTTTGTCCGAAGCAT‑CGGCCCGGTAGCGGGCCTCAGCCCCGAGTTGCAGCCAGGCCGACTCTGACAGCCGATGGTAGCGCAGGCCATCGAAAGGGGCAGTGGCCTTGGCGGGGTCGTCGAGAAAGCGGTAGTCCTCGATCCAGCGAGGGGCCCCTGGCCTCCCCGGC  >  NC_002947/3788298‑3788585

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