Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,324,838 T→A 5.3% Y31* (TAT→TAA) ‡ yciQ → enhancer of membrane protein expression; putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,324,8380TA 5.3% 255.2 / 4.4 75Y31* (TAT→TAA) ‡yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base A (2/2);  total (39/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

GCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCC  >  NC_000913/1324690‑1324984
                                                                                                                                                     |                                                                                                                                                  
gCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑At                                                                                                                                                    >  1:481914/1‑149 (MQ=255)
 cGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATg                                                                                                                                                   <  2:350862/149‑1 (MQ=255)
  gTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATga                                                                                                                                                  >  1:478196/1‑149 (MQ=255)
      gTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaagaa                                                                                                                                              <  2:243730/149‑1 (MQ=255)
      gTGATGTGAAGCCTTTCTTATAACAAAGTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATgaagaa                                                                                                                                              <  1:267855/149‑1 (MQ=255)
        gATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAgg                                                                                                                                            <  2:213220/149‑1 (MQ=255)
           gtgAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGAt                                                                                                                                         <  2:427322/149‑1 (MQ=255)
             gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTCTGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGTGAtct                                                                                                                                       >  2:255903/1‑149 (MQ=255)
             gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGAtct                                                                                                                                       >  3:418115/1‑149 (MQ=255)
              aaGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGAtctc                                                                                                                                      >  2:243670/1‑149 (MQ=255)
               aGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGAtctcc                                                                                                                                     >  1:350934/1‑148 (MQ=255)
               aGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑CATAAGAAGGGAtctcc                                                                                                                                     <  2:350934/149‑2 (MQ=255)
                    ttCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCtt                                                                                                                                >  3:306479/1‑149 (MQ=255)
                     tCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCttt                                                                                                                               <  3:559836/149‑1 (MQ=255)
                        tataACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGAt                                                                                                                            <  3:77247/149‑1 (MQ=255)
                        tataACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGAt                                                                                                                            <  1:65372/149‑1 (MQ=255)
                        tataACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGAt                                                                                                                            <  1:243670/149‑1 (MQ=255)
                          taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                          >  2:554726/1‑149 (MQ=255)
                           aaCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGtt                                                                                                                         >  1:319978/1‑149 (MQ=255)
                            aCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGCAAGCTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGAACCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTg                                                                                                                        >  1:399860/1‑149 (MQ=255)
                              aaaTTGCGTACCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTATGTTAGGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTCaa                                                                                                                      <  2:555059/149‑1 (MQ=255)
                                aTTGCGTTCCACGGATGGAGGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGc                                                                                                                    <  3:376758/149‑1 (MQ=255)
                                aTTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGc                                                                                                                    >  3:256237/1‑149 (MQ=255)
                                   gCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAc                                                                                                                 <  2:454143/149‑1 (MQ=255)
                                    cGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                                                                                                <  4:612822/149‑1 (MQ=255)
                                      ttCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCActct                                                                                                              <  1:462022/149‑1 (MQ=255)
                                          aCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCaa                                                                                                          >  4:377455/1‑149 (MQ=255)
                                             gATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                                                                                                       >  3:318559/1‑149 (MQ=255)
                                                ggAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGAt                                                                                                    <  4:14308/149‑1 (MQ=255)
                                                 gAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATg                                                                                                   <  2:478196/149‑1 (MQ=255)
                                                  aaGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                  >  3:17681/1‑149 (MQ=255)
                                                  aaGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                  <  2:180363/149‑1 (MQ=255)
                                                  aaGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                  >  1:192875/1‑149 (MQ=255)
                                                   aGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATggg                                                                                                 <  4:431754/149‑1 (MQ=255)
                                                      cTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATAATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTg                                                                                              >  4:622138/1‑149 (MQ=255)
                                                          ggCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                                                                                                       <  4:111451/120‑1 (MQ=255)
                                                          ggCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                                                                                                       >  3:111451/1‑120 (MQ=255)
                                                                    tttCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa                                                                                >  3:81203/1‑149 (MQ=255)
                                                                    tttCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa                                                                                <  4:522788/149‑1 (MQ=255)
                                                                         cTGCGTTTTGCTGTTGAAAGCAGGGCTGTATGTATCATCTCTAAGTTATGCAGAAGACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGGTCTTGGATGTTAAAGCCACTCTCAAACCAGTTGGAGTCCTCGACATCAAAGaaaaaa                                                                           <  2:204001/149‑3 (MQ=255)
                                                                         cTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACCGAGATCCCTTCTA‑ATGAAGAAGGAATCTCGCTCTTTGATGTTGAAGTCACTCTGCAACAAGATGGGGTGCTCGACATCAAAGAAAata                                                                           >  4:56634/1‑149 (MQ=255)
                                                                                   cTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCa                                                                 <  3:377455/149‑1 (MQ=255)
                                                                                    tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                                <  2:77927/149‑1 (MQ=255)
                                                                                            gCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGaaa                                                        <  2:488156/149‑1 (MQ=255)
                                                                                              aGGGCTGTTAGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAACGAAAATATTCATTTTCAGGCGCGAAAtg                                                      >  4:438691/1‑148 (MQ=255)
                                                                                                ggCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCAT‑AAGAAGGGAtc                                                                                                                                        <  1:286645/65‑1 (MQ=25)
                                                                                                ggCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCAT‑AAGAAGGGAtc                                                                                                                                        >  2:286645/1‑65 (MQ=25)
                                                                                                 gNTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATcagc                                                   >  3:255992/1‑149 (MQ=255)
                                                                                                        gTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTaa                                            >  4:366238/1‑149 (MQ=255)
                                                                                                         tATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAg                                           <  3:366238/149‑1 (MQ=255)
                                                                                                         tATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAg                                           >  4:36132/1‑149 (MQ=255)
                                                                                                          atcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGc                                          >  2:69758/1‑149 (MQ=255)
                                                                                                          atcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGc                                          <  1:69758/149‑1 (MQ=255)
                                                                                                               ctctAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGc                                     >  2:809/1‑149 (MQ=255)
                                                                                                                tctAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCt                                    >  1:665137/1‑149 (MQ=255)
                                                                                                                tctAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCt                                    >  2:638581/1‑149 (MQ=255)
                                                                                                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                                  >  2:221742/1‑149 (MQ=255)
                                                                                                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                                  >  1:540324/1‑149 (MQ=255)
                                                                                                                   aaGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCtttt                                 <  1:609062/149‑1 (MQ=255)
                                                                                                                    aGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTa                                <  2:523556/149‑1 (MQ=255)
                                                                                                                      ttATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATc                              <  3:632833/149‑1 (MQ=255)
                                                                                                                         tGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTg                           >  1:675957/1‑149 (MQ=255)
                                                                                                                          gCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGa                          >  1:77349/1‑149 (MQ=255)
                                                                                                                                  cacGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCAc                  <  4:255992/149‑1 (MQ=255)
                                                                                                                                     ggAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCGGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGAc               <  4:418115/149‑1 (MQ=255)
                                                                                                                                      gagaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACt              >  3:612678/1‑149 (MQ=255)
                                                                                                                                       agaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTa             <  4:306479/149‑1 (MQ=255)
                                                                                                                                       agaTCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTa             <  2:481914/149‑1 (MQ=255)
                                                                                                                                          tCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg          >  4:600217/1‑149 (MQ=255)
                                                                                                                                           cccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGa         >  2:181497/1‑149 (MQ=255)
                                                                                                                                            ccTTCTT‑CATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGAt        >  1:23366/1‑149 (MQ=255)
                                                                                                                                            ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGt                            <  1:431522/129‑1 (MQ=255)
                                                                                                                                            ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGt                            >  2:431522/1‑129 (MQ=255)
                                                                                                                                            ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGAt        >  1:655013/1‑149 (MQ=255)
                                                                                                                                            ccTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTGTATCGTGATTTACCACGACTATGGAt        >  4:571238/1‑149 (MQ=255)
                                                                                                                                                  t‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGcc  >  4:286038/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTT‑ATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCC  >  NC_000913/1324690‑1324984

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑