Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,659,640 G→A 5.5% E524E (GAG→GAA ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,659,6400GA 5.5% 242.8 / 2.9 73E524E (GAG→GAAynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (34/35);  new base A (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.32e-01

ATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTC  >  NC_000913/1659501‑1659781
                                                                                                                                           |                                                                                                                                             
aTCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGa                                                                                                                                       <  1:109120/148‑1 (MQ=255)
  cAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGata                                                                                                                                     >  7:262165/1‑148 (MQ=255)
    gCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtatt                                                                                                                                   <  7:6432/148‑1 (MQ=255)
    gCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtatt                                                                                                                                   <  1:194782/148‑1 (MQ=255)
    gCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtatt                                                                                                                                   <  7:305154/148‑1 (MQ=255)
    gCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtatt                                                                                                                                   <  1:40663/148‑1 (MQ=255)
     cATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGc                                                                                                                                              >  2:38144/1‑136 (MQ=255)
     cATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGc                                                                                                                                              <  1:38144/136‑1 (MQ=255)
      aTTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtattat                                                                                                                                 <  7:292403/148‑1 (MQ=255)
       ttCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATattatt                                                                                                                                <  6:262280/148‑1 (MQ=255)
        tCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTc                                                                                                                               >  2:271465/1‑148 (MQ=255)
               tCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGa                                                                                                                        <  3:215777/148‑1 (MQ=255)
               tCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATAGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGa                                                                                                                        >  2:35448/1‑148 (MQ=255)
                caacaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtatt                                                                                                                                   >  2:5333/1‑136 (MQ=255)
                caacaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtatt                                                                                                                                   <  1:5333/136‑1 (MQ=255)
                caacaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATattatt                                                                                                                                >  1:37700/1‑139 (MQ=255)
                caacaaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATattatt                                                                                                                                <  2:37700/139‑1 (MQ=255)
                   caaAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTAc                                                                                                                    >  6:286568/1‑148 (MQ=255)
                      aaCGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGcc                                                                                                                 >  3:18735/1‑148 (MQ=255)
                            tGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAAc                                                                                                           <  1:271465/148‑1 (MQ=255)
                                aTTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAACAGGAAGatat                                                                                                                                    <  7:74542/119‑1 (MQ=255)
                                aTTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAACAGGAAGatat                                                                                                                                    <  7:79204/119‑1 (MQ=255)
                                aTTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAACAGGAAGatat                                                                                                                                    >  8:74542/1‑119 (MQ=255)
                                aTTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAACAGGAAGatat                                                                                                                                    >  8:79204/1‑119 (MQ=255)
                                      cAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATAtc                                                                                                 >  5:258947/1‑148 (MQ=255)
                                        ggACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATAtctc                                                                                               <  7:331348/148‑1 (MQ=255)
                                          aCGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCAt                                                                                             <  8:262165/148‑1 (MQ=255)
                                                cgaaatgcgaaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTa                                                                                                                            <  5:51175/111‑1 (MQ=255)
                                                cgaaatgcgaaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTa                                                                                                                            >  6:51175/1‑111 (MQ=255)
                                                       cgaaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCt                                                                                >  5:5661/1‑148 (MQ=255)
                                                         aaatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGt                                                                              >  6:107705/1‑148 (MQ=255)
                                                          aatgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCAGTTTCCTCCAGCCTGTc                                                                             >  6:17097/1‑148 (MQ=255)
                                                           atgATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAAAATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTAATCCTAACGACTACGCCGGTAACATGGGATATCTCATTATCCACCAGCCTCTCa                                                                            >  8:205950/1‑148 (MQ=255)
                                                                       aTCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGaaaaa                                                                 <  1:35448/147‑1 (MQ=255)
                                                                         cGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTc                                                              >  5:180097/1‑148 (MQ=255)
                                                                         cGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTc                                                              >  1:117891/1‑148 (MQ=255)
                                                                           aaaaCTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGTATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGa                                                            <  5:113318/148‑1 (MQ=255)
                                                                            aaaCTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGaa                                                           >  5:299885/1‑148 (MQ=255)
                                                                              aCTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtattat                                                                                                                                 >  1:87487/1‑76 (MQ=255)
                                                                              aCTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGAtattat                                                                                                                                 <  2:87487/76‑1 (MQ=255)
                                                                                   aTGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCaaa                                                     <  6:5661/147‑1 (MQ=255)
                                                                                      aCCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCg                                                  <  4:18735/147‑1 (MQ=255)
                                                                                      aCCTCATCGGAAAAATATGCCAACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATCGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGa                                                 >  4:208632/1‑148 (MQ=255)
                                                                                            tCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACt                                           >  5:182968/1‑148 (MQ=255)
                                                                                              ggCAAAATATGCCGACATTCTGCTGCCATAACTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTg                                          <  3:38561/147‑1 (MQ=255)
                                                                                              ggCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATggg                                                                                                      >  1:131411/1‑87 (MQ=255)
                                                                                              ggCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATggg                                                                                                      <  2:131411/87‑1 (MQ=255)
                                                                                                 aaaaTATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATc                                      >  2:276850/1‑148 (MQ=255)
                                                                                                    atatGCCGACATTCTGCTGCCAGACCTGAGGACCGTTGAGCAGTAAGCTATTATTCCTAACGACTACGCCGTTAACATGGGTTATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTg                                   <  7:232051/148‑1 (MQ=255)
                                                                                                     tatGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGa                                  <  3:105942/148‑1 (MQ=255)
                                                                                                       tGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAg                                 <  6:299885/147‑1 (MQ=255)
                                                                                                          cGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTg                               >  2:183047/1‑146 (MQ=255)
                                                                                                           gACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAg                            >  1:117203/1‑148 (MQ=255)
                                                                                                            aCATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGt                           >  1:109212/1‑148 (MQ=255)
                                                                                                             cATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTc                          >  4:40311/1‑148 (MQ=255)
                                                                                                                tCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTcgc                        >  7:173750/1‑147 (MQ=255)
                                                                                                                tCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTcg                         <  8:146416/146‑1 (MQ=255)
                                                                                                                tCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTcg                         >  7:146416/1‑146 (MQ=255)
                                                                                                                 ctgctgCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTcgcg                       >  8:290659/1‑147 (MQ=255)
                                                                                                                    ctgcCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGaa                     <  2:117203/146‑1 (MQ=255)
                                                                                                                    ctgcCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGa                      <  7:290659/145‑1 (MQ=255)
                                                                                                                        cAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGAAACGTCt               <  6:180097/148‑1 (MQ=255)
                                                                                                                           aCCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGAGCCTGAGTGAAGTCGCGAAACGTCTTgg            <  6:182968/148‑1 (MQ=255)
                                                                                                                            ccTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGAAACGTCTTGGa           <  5:286568/148‑1 (MQ=255)
                                                                                                                                  gaCCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGa                                        >  4:97216/1‑113 (MQ=255)
                                                                                                                                  gaCCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGa                                        <  3:97216/113‑1 (MQ=255)
                                                                                                                                   aCCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACg    >  3:96531/1‑148 (MQ=255)
                                                                                                                                     cGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTc  <  6:258947/148‑1 (MQ=255)
                                                                                                                                      gTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAg                            <  2:100826/121‑1 (MQ=255)
                                                                                                                                      gTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAg                            >  1:100826/1‑121 (MQ=255)
                                                                                                                                      gTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTc  >  6:152308/1‑147 (MQ=255)
                                                                                                                                          aGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATcc                                     >  8:167781/1‑108 (MQ=255)
                                                                                                                                          aGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATcc                                     <  7:167781/108‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                             
ATCAGCATTCTGACATCAACAAAACGCATGAAATTCTGCAGGACGAATCGAAATGCGAAATGATTGTGGTCATCGAAAACTTTATGACCTCATCGGCAAAATATGCCGACATTCTGCTGCCAGACCTGATGACCGTTGAGCAGGAAGATATTATTCCTAACGACTACGCCGGTAACATGGGATATCTCATTTTCCTCCAGCCTGTCACCAGCGAAAAATTCGAACGCAAACCGATTTACTGGATCCTGAGTGAAGTCGCGAAACGTCTTGGACCTGACGTC  >  NC_000913/1659501‑1659781

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