Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,360,643 T→A 8.9% intergenic (+27/‑533) gltD → / → gltF glutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,360,6430TA 8.9% 151.1 / 5.8 45intergenic (+27/‑533)gltD/gltFglutamate synthase, 4Fe‑4S protein, small subunit/periplasmic protein
Reads supporting (aligned to +/- strand):  ref base T (18/23);  new base A (2/2);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.50e-01

GGCGTTTGGTTTCCGTCCACACAACATGGAATGGCTGGCAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAA  >  NC_000913/3360388‑3360876
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
ggCGTTTGGTTTCCGTCCACACAACATGGAATGGCTGGCAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAAtt                                                                                                                                                                                                                        >  2:259860/1‑275 (MQ=255)
  cGTGTGGTTTCCGTCAACACAACATGGAATGGCTGGCAAAACCCAGCGTCGGGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCGGCAACCCGAAAATCTGTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGAACGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTa                                                                                                                                                                                                                       <  4:36608/274‑1 (MQ=255)
               tCCAAACAACAGGGAATGGCTGGCAAAAAACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTgcg                                                                                                                                                                                                          <  1:86773/274‑1 (MQ=255)
                   cacaACATGGAAGGGCTGGCAAAACACAGCTTCGGGCTGGATTCACAAGGCCGTATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAAGCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgtt                                                                                                                                                                                                      <  2:341171/274‑1 (MQ=255)
                                  ctggcAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGCCGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGAAGAAGTTCGCGTTGTTTTATGGCacaagc                                                                                                                                                                                      >  4:519706/1‑270 (MQ=255)
                                                 gtctggCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCGTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                                                                                                        <  2:332436/269‑1 (MQ=255)
                                                    gAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGAcc                                                                                                                                                                     <  1:34159/274‑1 (MQ=255)
                                                          gATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGttt                                                                                                                                                                         >  2:337297/1‑264 (MQ=255)
                                                              cacaAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGAcc                                                                                                                                                                     <  4:337296/264‑1 (MQ=255)
                                                               acaAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          <  4:259859/274‑1 (MQ=255)
                                                                                                            gACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACt                                                                                                                                                                                                                                           >  5:813909/1‑148 (MQ=255)
                                                                                                                   aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCtc                                                                                                                                                                                                                                    <  8:781445/148‑2 (MQ=255)
                                                                                                                   aaCCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTg                                                                                                                                                                                                                                    <  6:813909/148‑1 (MQ=255)
                                                                                                                                         gcgATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                                                                                                                                          <  6:400137/120‑1 (MQ=255)
                                                                                                                                         gcgATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                                                                                                                                          >  5:400137/1‑120 (MQ=255)
                                                                                                                                            atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCaa                                                                                                                                                                                                                           <  8:89625/132‑1 (MQ=255)
                                                                                                                                            atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCaa                                                                                                                                                                                                                           >  7:89625/1‑132 (MQ=255)
                                                                                                                                            atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAg                                                                                                                                                                                                              <  5:923303/145‑1 (MQ=255)
                                                                                                                                            atatCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAg                                                                                                                                                                                                              >  6:923303/1‑145 (MQ=255)
                                                                                                                                                              aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGt                                                          >  4:253415/1‑275 (MQ=255)
                                                                                                                                                              aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGa                                                                                                                                                                                                                                >  6:940360/1‑109 (MQ=255)
                                                                                                                                                              aTCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACACAGTCTGGCGa                                                                                                                                                                                                                                <  5:940360/109‑1 (MQ=255)
                                                                                                                                                                     ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtgtg                                                                                                                                                                                                                                        >  8:291268/1‑94 (MQ=255)
                                                                                                                                                                     ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtgtg                                                                                                                                                                                                                                        <  7:291268/94‑1 (MQ=255)
                                                                                                                                                                     ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTgg                                                                                                                                                                                                                                   <  8:450405/99‑1 (MQ=255)
                                                                                                                                                                     ggtgACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTgg                                                                                                                                                                                                                                   >  7:450405/1‑99 (MQ=255)
                                                                                                                                                                        gACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACctc                                                 <  2:253416/274‑1 (MQ=255)
                                                                                                                                                                         aCCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                                                                                                                                       >  2:436884/1‑155 (MQ=255)
                                                                                                                                                                               aTTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                                                                                                                                          >  8:691560/1‑82 (MQ=255)
                                                                                                                                                                               aTTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACtg                                                                                                                                                                                                                                          <  7:691560/82‑1 (MQ=255)
                                                                                                                                                                                    cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                                                                                                                               >  5:129991/1‑120 (MQ=255)
                                                                                                                                                                                    cGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCt                                                                                                                                                                                               <  6:129991/120‑1 (MQ=255)
                                                                                                                                                                                        ggTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCta                                                                                                                                                              >  2:365869/1‑149 (MQ=255)
                                                                                                                                                                                            gTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          <  4:365868/149‑1 (MQ=255)
                                                                                                                                                                                            gTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCaa                             <  2:519707/274‑1 (MQ=255)
                                                                                                                                                                                             tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                                                                                                                                       <  5:640810/135‑1 (MQ=255)
                                                                                                                                                                                             tAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGa                                                                                                                                                                       >  6:640810/1‑135 (MQ=255)
                                                                                                                                                                                                               tatGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttg                                                                                                                                                                                                     >  6:212446/1‑87 (MQ=255)
                                                                                                                                                                                                               tatGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCgttg                                                                                                                                                                                                     <  5:212446/87‑1 (MQ=255)
                                                                                                                                                                                                                    aCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                                                                                                                              <  6:486211/137‑1 (MQ=255)
                                                                                                                                                                                                                    aCTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAAtttt                                                                                                                                              >  5:486211/1‑137 (MQ=255)
                                                                                                                                                                                                                      tggctggAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTaa  >  2:161052/1‑275 (MQ=255)
                                                                                                                                                                                                                                 ttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          >  7:386492/1‑112 (MQ=255)
                                                                                                                                                                                                                                 ttAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTataat                                                                                                                                                          <  8:386492/112‑1 (MQ=255)
                                                                                                                                                                                                                                                   aaTTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGCCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTTTTCACCAAGTCTGATGAGTGAAGCACAGAAATCAATTTCTCATTTACAATCCGTCGCa                                                                   >  3:11206/1‑181 (MQ=255)
                                                                                                                                                                                                                                                    aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                                                                                                        >  7:507159/1‑79 (MQ=38)
                                                                                                                                                                                                                                                    aTTCGCCAGACACAGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTg                                                                                                                                                                        <  8:507159/79‑1 (MQ=38)
                                                                                                                                                                                                                                                     ttCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTc                                                                                                                                                                                                >  6:110761/1‑54 (MQ=255)
                                                                                                                                                                                                                                                     ttCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTc                                                                                                                                                                                                <  5:110761/54‑1 (MQ=255)
                                                                                                                                                                                                                                                       cGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGAAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTa                                                               <  1:11206/181‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtga                                                                                                               <  6:301873/127‑1 (MQ=255)
                                                                                                                                                                                                                                                             aCTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCtga                                                                                                               >  5:301873/1‑127 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                         
GGCGTTTGGTTTCCGTCCACACAACATGGAATGGCTGGCAAAACACAGCGTCGAGCTGGATTCACAAGGCCGCATCATCGCCCCGGAAGGCAGCGACAACGCCTTCCAGACCAGCAACCCGAAAATCTTTGCTGGCGGCGATATCGTCCGTGGTTCCGATCTGGTGGTGACCGCTATTGCCGAAGGTCGTAAGGCGGCAGACGGTATTATGAACTGGCTGGAAGTTTAAGCGAGGTAACAATGAATTCGCCAGACTGTGTCTGGCGAATCAAATTACTGACGAAGTGCGCGTTGTTTTCTGGCAAATGCGGCATGAGCGTTTGACCTTTCCTATAATAGGGATAATTTTCAATAAATTGTGCTTTGTTCACCAAGTCTGATGAGTGAAGCACATCAAGCAATTTCTCATTTACAATCCGTCTCAATTACTTGTTAACACCTCTGACAATACCGATTCTGCAAGCTATTCAGTCTTTTCATTCCATGTAA  >  NC_000913/3360388‑3360876

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: