Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,742,003	A→T	12.7%	intergenic (‑28/+385)	PP_3305 ← / ← PP_3306	TerC family membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,742,003	0	A	T	12.7%	168.0 / 15.6	63	intergenic (‑28/+385)	PP_3305/PP_3306	TerC family membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/28);  new base T (3/5);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

CCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTA  >  NC_002947/3741854‑3742148                                                                                                                                                       |                                                                                                                                                  ccGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaa                                                                                                                                                   >  1:774274/1‑150 (MQ=255)  cgcgGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGaac                                                                                                                                                  <  1:826884/150‑1 (MQ=255)       tGGAGCACACCTAAGTCAAGCGCGAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTa                                                                                                                                             <  1:492623/150‑1 (MQ=255)        ggAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTAt                                                                                                                                            <  2:774274/150‑1 (MQ=255)          cgatctccTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                          <  1:306405/144‑1 (MQ=255)          aGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGa                                                                                                                                          <  2:815236/150‑1 (MQ=255)          aGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                          <  1:736291/150‑1 (MQ=255)               acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGATATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCg                                                                                                                                     >  1:349887/1‑150 (MQ=255)               acCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCg                                                                                                                                     >  2:392146/1‑150 (MQ=255)                ccTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGa                                                                                                                                    >  2:826513/1‑150 (MQ=255)                 cTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTAATAGATGTACATCTACGAcacatcg                                                                                                                                   >  2:306405/1‑144 (MQ=255)                  tAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACa                                                                                                                                  <  1:935570/150‑1 (MQ=255)                    aGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATg                                                                                                                                <  1:826513/150‑1 (MQ=255)                      tCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtataca                                                                                                                                   <  1:1086822/145‑6 (MQ=255)                      tCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGACtataca                                                                                                                                   >  2:1086822/1‑140 (MQ=255)                          gcgcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTacga                                                                                                                                          >  2:709562/1‑131 (MQ=255)                          gcgcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGa                                                                                                                                          <  1:709562/134‑1 (MQ=255)                          gcgcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTAACAAGAACG‑TCTTGGTAGAGATACGTAATAAGGGCAGACATGAGCCCCTCCTGTATAGGAGTAGTTGAACAAATATGACTCCGACATGGCGGta                                                                                                                          >  2:1097291/1‑150 (MQ=255)                            gcCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGtata                                                                                                                        >  1:1285690/1‑150 (MQ=255)                                 cAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGt                                                                                                                   <  1:861750/150‑1 (MQ=255)                                             aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAAc                                                                                                                      <  2:1115371/135‑1 (MQ=255)                                             aacaGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAAc                                                                                                                      >  1:1115371/1‑135 (MQ=255)                                                   gAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCa                                                                                                 <  1:1150205/150‑1 (MQ=255)                                                          gAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATaa                                                                                          >  2:387711/1‑150 (MQ=255)                                                               ccaaagccaGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAAt                                                                                     >  1:92521/1‑150 (MQ=255)                                                                 aaagccaGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAAttt                                                                                   >  1:513236/1‑150 (MQ=255)                                                                           gTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGcgcg                                                                         >  2:1073737/1‑150 (MQ=255)                                                                                tACCAAG‑ACGTTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAg                                                                    >  1:242220/1‑150 (MQ=255)                                                                                tACCAAG‑ACGTTCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAg                                                                    >  1:674825/1‑150 (MQ=255)                                                                                 actaaGAACG‑TCTTGTTAGAGATACGGACTAAGGGCAGACATTTGCCCCTCTTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGTTTCCGGTGTCACTAATAAAGAATTTATTGAGCGCGAGTAgg                                                                   <  2:1002795/147‑1 (MQ=255)                                                                                    aaGAACG‑TCTTGGTAGAGATACGTACTNAGGGCAGACATGNGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTg                                                                >  1:1355131/1‑150 (MQ=255)                                                                                      gAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCa                                                              >  1:214235/1‑150 (MQ=255)                                                                                          g‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGAc                                                                                                                                         >  1:434061/1‑71 (MQ=38)                                                                                          g‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCATAGTTGTTCAACTACGAc                                                                                                                                         <  2:434061/71‑1 (MQ=37)                                                                                              ttGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATAGCGGTATAACCTTCAAAGATGATGAGCTGGGCCGATAGAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaa                                                       <  1:387711/150‑1 (MQ=255)                                                                                                ggTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGaaaa                                                     <  1:270110/150‑1 (MQ=255)                                                                                                    gagaTACGTACTAAGGGCAGATATGAGCCCCGCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGCGACGCTGAACTGCACTAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACAtt                                                 <  2:349887/150‑1 (MQ=255)                                                                                                    gagaTACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACAtt                                                 <  2:674825/150‑1 (MQ=255)                                                                                                          cGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACaaa                                           <  2:1060910/150‑1 (MQ=255)                                                                                                                aaGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGCCTCCGACATGGCGGTATAACCGTCAGACGCGCCCGGCGTCAAGAAGCTACGACGAATCCGGCGCGATTGGTTTGCAAGTAGAAAACATTTACAAATCTTGt                                     <  1:1097291/150‑1 (MQ=255)                                                                                                                    gCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATc                                                                                <  1:622569/103‑1 (MQ=255)                                                                                                                    gCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATc                                                                                >  2:622569/1‑103 (MQ=255)                                                                                                                    gCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCAt                                 <  2:92521/150‑1 (MQ=255)                                                                                                                          aTGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAg                           >  1:226289/1‑150 (MQ=255)                                                                                                                                   tcctGGATAGTCGTATTTGAACAACTATGCCTCCGAGATCGCGGTATAACCGTCAGAGGGGACCGGCTTCACGAATAACGAATTTATCGAGCGCGAGTACGCTGCAAGTAGAAAACATGTACAAATATTGTCCATGACAAGCTAAAACGc                  <  2:242482/150‑1 (MQ=255)                                                                                                                                    cctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCt                 <  2:1355131/150‑1 (MQ=255)                                                                                                                                     ctGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCtt                <  2:1285690/150‑1 (MQ=255)                                                                                                                                       gTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCtaga              >  1:242482/1‑147 (MQ=255)                                                                                                                                        tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                          >  2:1167528/1‑105 (MQ=255)                                                                                                                                        tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTa                                                          <  1:1167528/105‑1 (MQ=255)                                                                                                                                        tataGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCa             <  2:242220/150‑1 (MQ=255)                                                                                                                                            gTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGa         >  2:444990/1‑150 (MQ=255)                                                                                                                                             tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTc                                                                                                                  <  1:98138/44‑1 (MQ=255)                                                                                                                                             tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTc                                                                                                                  >  2:98138/1‑44 (MQ=255)                                                                                                                                             tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                              <  1:185669/48‑1 (MQ=255)                                                                                                                                             tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCagag                                                                                                              >  2:185669/1‑48 (MQ=255)                                                                                                                                             tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTc                                    <  2:459721/122‑1 (MQ=255)                                                                                                                                             tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTc                                    >  1:459721/1‑122 (MQ=255)                                                                                                                                             tCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGaa        >  2:1180952/1‑150 (MQ=255)                                                                                                                                                  gTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACaaa                                           >  1:8807/1‑110 (MQ=255)                                                                                                                                                  gTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACaaa                                           <  2:8807/110‑1 (MQ=255)                                                                                                                                                  gTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTg               >  1:1104104/1‑138 (MQ=255)                                                                                                                                                  gTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTg               <  2:1104104/138‑1 (MQ=255)                                                                                                                                                  gTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAaact   >  1:911750/1‑150 (MQ=255)                                                                                                                                                   ttGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAACCATTTACAAATCTTGTACATGCCAATCTAAAACGCTTGCAGTGAAAaacta  <  1:392146/150‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  CCGCGGTGGAGCACACCTAAGTCAAGCGCCAACAGTGTTAGCACAACAGCGAAAAAGGAAAGCCAAAGCCAGAGGTCAGTACCAAGAACG‑TCTTGGTAGAGATACGTACTAAGGGCAGACATGAGCCCCTCCTGTATAGTCGTAGTTGAACAACTATGACTCCGACATGGCGGTATAACCGTCAGAGGGGCCCGGCGTCACGAATAAAGAATTTATCGAGCGCGAGTAGGCTGCAAGTAGAAAACATTTACAAATCTTGTCCATGACAAGCTAAAACGCTTGCAGTGAAAAACTA  >  NC_002947/3741854‑3742148

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