Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,090,056	C→G	6.5%	S36C (TCC→TGC)	gudP →	glucarate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,090,056	0	C	G	6.5%	194.7 / 3.6	61	S36C (TCC→TGC)	gudP	glucarate transporter
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/2);  total (31/30)
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.49e-01

CCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTT  >  NC_002947/4089915‑4090203                                                                                                                                               |                                                                                                                                                      cccGCACAAAAATAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTg                                                                                                                                                <  2:261228/149‑1 (MQ=255)  ccGCACAAAAATAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTG‑‑c                                                                                                                                             >  6:85761/1‑149 (MQ=255)      acaAAAATAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTG‑‑cgata                                                                                                                                         <  5:85761/149‑5 (MQ=255)      acaAAAAAAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTGGAGCCTG‑‑cgata                                                                                                                                         <  6:55649/149‑5 (MQ=255)         aaaaTAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGa                                                                                                                                      >  4:188960/1‑149 (MQ=255)            aTAACTAGTGGGAGTATTGCCAGCATGGCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGAcc                                                                                                                                    >  1:14601/1‑149 (MQ=255)             tAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTc                                                                                                                                  <  6:178244/149‑1 (MQ=255)              aaCTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCg                                                                                                                                 <  5:170383/149‑1 (MQ=255)                 tAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCa                                                                                                                              >  2:263740/1‑149 (MQ=255)                   gTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCAtt                                                                                                                            >  2:169245/1‑149 (MQ=255)                   gTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTACTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCAtt                                                                                                                            >  2:206008/1‑149 (MQ=255)                        aGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGc                                                                                                                       <  4:242974/149‑1 (MQ=255)                        aGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGc                                                                                                                       >  1:245038/1‑149 (MQ=255)                         gTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACGTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGcc                                                                                                                      >  3:178990/1‑149 (MQ=255)                          tATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCa                                                                                                                              <  4:211744/140‑1 (MQ=255)                          tATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCa                                                                                                                              >  3:211744/1‑140 (MQ=255)                          tATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCg                                                                                                                     >  4:295605/1‑149 (MQ=255)                               ccAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAAcc                                                                                                                <  7:220981/149‑1 (MQ=255)                                  gCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTc                                                                                                             >  4:53050/1‑149 (MQ=255)                                  gCAT‑GCAAGCGACTAAGAAAACGCATGTGCGATACCTGATCCTGTTCATGCTGTTTCTGGTGACAACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTc                                                                                                             >  7:83007/1‑149 (MQ=255)                                          atcGACTAAGAAAACGCATGTTCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTtata                                                                                                      <  4:202335/147‑1 (MQ=255)                                                 aaGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCgg                                                                                                                                <  1:306624/116‑1 (MQ=255)                                                 aaGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCgg                                                                                                                                >  2:306624/1‑116 (MQ=255)                                                   gAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATATTTTCCg                                                                                             >  5:276093/1‑149 (MQ=255)                                                        cGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTc                                                                                        <  4:178990/149‑1 (MQ=255)                                                            tgtgCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGAt                                                                                    <  3:295605/149‑1 (MQ=255)                                                            tgtgCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGAt                                                                                    <  2:14601/149‑1 (MQ=255)                                                                cgcTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGc                                                                                <  1:169245/149‑1 (MQ=255)                                                                   tACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCATG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCATTCGGATGTGCTTa                                                                             >  7:125962/1‑149 (MQ=255)                                                                         aTCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTGCACAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACg                                                                           <  1:263740/147‑1 (MQ=255)                                                                            cTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGTGCTTACGtggctg                                                                     <  8:83007/148‑1 (MQ=255)                                                                             tGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGtggctggc                                                                   >  1:193337/1‑149 (MQ=255)                                                                               ttCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGc                                                                                            >  3:228674/1‑122 (MQ=255)                                                                               ttCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGc                                                                                            <  4:228674/122‑1 (MQ=255)                                                                                tCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGcc                                                                                                                      >  5:146733/1‑95 (MQ=255)                                                                                tCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGcc                                                                                                                      <  6:146733/95‑1 (MQ=255)                                                                                    gCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGAtt                                                             >  6:49075/1‑148 (MQ=255)                                                                                           cTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGcc                                                                                                                      <  6:59528/84‑1 (MQ=255)                                                                                           cTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGcc                                                                                                                      >  5:59528/1‑84 (MQ=255)                                                                                              gTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCtggtgg                                                     <  5:49075/146‑1 (MQ=255)                                                                                                gACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CACAAAGACCTCGTCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGtggctggc                                                >  5:228136/1‑149 (MQ=255)                                                                                                  ccACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCAGCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTa                                              >  1:234527/1‑149 (MQ=255)                                                                                                    aCGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACt                                            >  1:159050/1‑149 (MQ=255)                                                                                                    aCGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACt                                            <  2:159050/149‑1 (MQ=255)                                                                                                          aaCTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATc                                       <  3:136851/148‑1 (MQ=255)                                                                                                              aTGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTc                                  <  3:188960/149‑1 (MQ=255)                                                                                                               tGCCGACCGCGCGACCATTGCTAGCGCAGGCTCCAGCCTG‑‑CAGAAAGTCCTCGTCATTGTCGCCGTAACGCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTc                                  <  8:125962/148‑1 (MQ=255)                                                                                                                 ccGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGt                               >  3:167422/1‑149 (MQ=255)                                                                                                                           gACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGtggctg                                                  <  8:308799/120‑1 (MQ=255)                                                                                                                           gACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGtggctg                                                  >  7:308799/1‑120 (MQ=255)                                                                                                                              cATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTc                   >  5:288622/1‑148 (MQ=255)                                                                                                                              cATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTc                   >  5:169048/1‑148 (MQ=255)                                                                                                                                ttGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTAc                <  2:234527/149‑1 (MQ=255)                                                                                                                                ttGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTAc                <  3:53050/149‑1 (MQ=255)                                                                                                                                    tATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACg                     <  7:194264/140‑1 (MQ=255)                                                                                                                                    tATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACg                     >  8:194264/1‑140 (MQ=255)                                                                                                                                    tATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCt            >  8:178108/1‑149 (MQ=255)                                                                                                                                    tATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCt            >  4:271334/1‑149 (MQ=255)                                                                                                                                      tCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCtt           <  2:193337/148‑1 (MQ=255)                                                                                                                                           ggCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCAt     <  4:167422/149‑1 (MQ=255)                                                                                                                                              tCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCtt  >  2:195760/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                      CCCGCACAAAAATAACTAGTGGGAGTATTGCCAGCAT‑GCAAGCGACTAAGAAAACGCATGTGCGCTACCTGATCCTGTTCATGCTGTTTCTGGTGACCACGATCAACTATGCCGACCGCGCGACCATTGCTATCGCAGGCTCCAGCCTG‑‑CAGAAAGACCTCGGCATTGACGCCGTAACCCTCGGTTATATCTTTTCCGCCTTCGGATGGGCTTACGTGGCTGGCCAGATTCCTGGTGGCTGGCTACTCGATCGCTTCGGTTCCAAGAACGTCTACGCCTTCAGCATCTT  >  NC_002947/4089915‑4090203

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