Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,324,837 A→C 5.3% Y31S (TAT→TCT) ‡ 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,8370AC 5.3% 254.8 / 7.2 76Y31S (TAT→TCT) ‡yciQenhancer of membrane protein expression; putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (38/34);  new base C (2/2);  total (40/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 = 5.25e-01

GCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCC  >  NC_000913/1324690‑1324984
                                                                                                                                                    |                                                                                                                                                   
gCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTAt                                                                                                                                                    >  1:481914/1‑149 (MQ=255)
 cGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATg                                                                                                                                                   <  2:350862/149‑1 (MQ=255)
  gTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATga                                                                                                                                                  >  1:478196/1‑149 (MQ=255)
      gTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATgaagaa                                                                                                                                              <  2:243730/149‑1 (MQ=255)
      gTGATGTGAAGCCTTTCTTATAACAAAGTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATgaagaa                                                                                                                                              <  1:267855/149‑1 (MQ=255)
        gATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAgg                                                                                                                                            <  2:213220/149‑1 (MQ=255)
           gtgAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGAt                                                                                                                                         <  2:427322/149‑1 (MQ=255)
             gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTCTGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGTGAtct                                                                                                                                       >  2:255903/1‑149 (MQ=255)
             gAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtct                                                                                                                                       >  3:418115/1‑149 (MQ=255)
              aaGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGAtctc                                                                                                                                      >  2:243670/1‑149 (MQ=255)
               aGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGAtctcc                                                                                                                                     >  1:350934/1‑148 (MQ=255)
               aGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGAtctcc                                                                                                                                     <  2:350934/149‑2 (MQ=255)
                    ttCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCtt                                                                                                                                >  3:306479/1‑149 (MQ=255)
                     tCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCttt                                                                                                                               <  3:559836/149‑1 (MQ=255)
                        tataACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGAt                                                                                                                            <  3:77247/149‑1 (MQ=255)
                        tataACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGAt                                                                                                                            <  1:65372/149‑1 (MQ=255)
                        tataACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGAt                                                                                                                            <  1:243670/149‑1 (MQ=255)
                          taACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGt                                                                                                                          >  2:554726/1‑149 (MQ=255)
                           aaCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGtt                                                                                                                         >  1:319978/1‑149 (MQ=255)
                            aCAAATTGCGTTCCACGGATGGAAGACTATGGCAGGCAAGCTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGAACCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTg                                                                                                                        >  1:399860/1‑149 (MQ=255)
                              aaaTTGCGTACCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTATGTTAGGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTCaa                                                                                                                      <  2:555059/149‑1 (MQ=255)
                                aTTGCGTTCCACGGATGGAGGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGc                                                                                                                    <  3:376758/149‑1 (MQ=255)
                                aTTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGc                                                                                                                    >  3:256237/1‑149 (MQ=255)
                                   gCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAc                                                                                                                 <  2:454143/149‑1 (MQ=255)
                                    cGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCAct                                                                                                                <  4:612822/149‑1 (MQ=255)
                                      ttCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCActct                                                                                                              <  1:462022/149‑1 (MQ=255)
                                          aCGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCaa                                                                                                          >  4:377455/1‑149 (MQ=255)
                                             gATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCa                                                                                                       >  3:318559/1‑149 (MQ=255)
                                                ggAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGAt                                                                                                    <  4:14308/149‑1 (MQ=255)
                                                 gAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATg                                                                                                   <  2:478196/149‑1 (MQ=255)
                                                  aaGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                  >  3:17681/1‑149 (MQ=255)
                                                  aaGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                  <  2:180363/149‑1 (MQ=255)
                                                  aaGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATgg                                                                                                  >  1:192875/1‑149 (MQ=255)
                                                   aGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATggg                                                                                                 <  4:431754/149‑1 (MQ=255)
                                                      cTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATAATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTg                                                                                              >  4:622138/1‑149 (MQ=255)
                                                          ggCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                                                                                                       <  4:111451/120‑1 (MQ=255)
                                                          ggCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGa                                                                                                                       >  3:111451/1‑120 (MQ=255)
                                                                    tttCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa                                                                                >  3:81203/1‑149 (MQ=255)
                                                                    tttCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGa                                                                                <  4:522788/149‑1 (MQ=255)
                                                                         cTGCGTTTTGCTGTTGAAAGCAGGGCTGTATGTATCATCTCTAAGTTATGCAGAAGACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGGTCTTGGATGTTAAAGCCACTCTCAAACCAGTTGGAGTCCTCGACATCAAAGaaaaaa                                                                           <  2:204001/149‑3 (MQ=255)
                                                                         cTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACCGAGATCCCTTCTAATGAAGAAGGAATCTCGCTCTTTGATGTTGAAGTCACTCTGCAACAAGATGGGGTGCTCGACATCAAAGAAAata                                                                           >  4:56634/1‑149 (MQ=255)
                                                                                   cTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCa                                                                 <  3:377455/149‑1 (MQ=255)
                                                                                    tGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAg                                                                <  2:77927/149‑1 (MQ=255)
                                                                                            gCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGaaa                                                        <  2:488156/149‑1 (MQ=255)
                                                                                              aGGGCTGTTAGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAACGAAAATATTCATTTTCAGGCGCGAAAtg                                                      >  4:438691/1‑148 (MQ=255)
                                                                                                ggCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGAtc                                                                                                                                        <  1:286645/65‑1 (MQ=25)
                                                                                                ggCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTCATAAGAAGGGAtc                                                                                                                                        >  2:286645/1‑65 (MQ=25)
                                                                                                 gNTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATcagc                                                   >  3:255992/1‑149 (MQ=255)
                                                                                                        gTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTaa                                            >  4:366238/1‑149 (MQ=255)
                                                                                                         tATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAg                                           <  3:366238/149‑1 (MQ=255)
                                                                                                         tATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAg                                           >  4:36132/1‑149 (MQ=255)
                                                                                                          atcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGc                                          >  2:69758/1‑149 (MQ=255)
                                                                                                          atcatcTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGc                                          <  1:69758/149‑1 (MQ=255)
                                                                                                               ctctAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGc                                     >  2:809/1‑149 (MQ=255)
                                                                                                                tctAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCt                                    >  1:665137/1‑149 (MQ=255)
                                                                                                                tctAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCt                                    >  2:638581/1‑149 (MQ=255)
                                                                                                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                                  >  2:221742/1‑149 (MQ=255)
                                                                                                                  tAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCttt                                  >  1:540324/1‑149 (MQ=255)
                                                                                                                   aaGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCtttt                                 <  1:609062/149‑1 (MQ=255)
                                                                                                                    aGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTa                                <  2:523556/149‑1 (MQ=255)
                                                                                                                      ttATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATc                              <  3:632833/149‑1 (MQ=255)
                                                                                                                         tGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTg                           >  1:675957/1‑149 (MQ=255)
                                                                                                                          gCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGa                          >  1:77349/1‑149 (MQ=255)
                                                                                                                                  cacGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCAc                  <  4:255992/149‑1 (MQ=255)
                                                                                                                                     ggAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCGGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGAc               <  4:418115/149‑1 (MQ=255)
                                                                                                                                      gagaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACt              >  3:612678/1‑149 (MQ=255)
                                                                                                                                       agaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTa             <  4:306479/149‑1 (MQ=255)
                                                                                                                                       agaTCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTa             <  2:481914/149‑1 (MQ=255)
                                                                                                                                          tCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATgg          >  4:600217/1‑149 (MQ=255)
                                                                                                                                           cccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGa         >  2:181497/1‑149 (MQ=255)
                                                                                                                                            ccTTCTTCATAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGAt        >  1:23366/1‑149 (MQ=255)
                                                                                                                                            ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGt                            <  1:431522/129‑1 (MQ=255)
                                                                                                                                            ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGt                            >  2:431522/1‑129 (MQ=255)
                                                                                                                                            ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGAt        >  1:655013/1‑149 (MQ=255)
                                                                                                                                            ccTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTGTATCGTGATTTACCACGACTATGGAt        >  4:571238/1‑149 (MQ=255)
                                                                                                                                                  tATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGcc  >  4:286038/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GCGTAGGTGATGTGAAGCCTTTCTTATAACAAATTGCGTTCCACGGATGGAAGACTATGGCAGGGAAGTTTCGCTGCATTTTGCTGTTGATAGCAGGGCTGTTTGTATCATCTCTAAGTTATGCAGAAAACACGGAGATCCCTTCTTATGAAGAAGGGATCTCGCTCTTTGATGTTGAAGCCACTCTGCAACCAGATGGGGTGCTCGACATCAAAGAAAATATTCATTTTCAGGCGCGAAATCAGCAGATTAAGCACGGCTTTTATCGTGATTTACCACGACTATGGATGCAGCC  >  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: