Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,810,065 T→G 15.3% intergenic (+17/‑146) yniC → / → ydjM 2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,810,0650TG15.3% 277.6 / 37.5 112intergenic (+17/‑146)yniC/ydjM2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (47/47);  new base G (7/10);  total (54/58)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.22e-01

Only 100 of 117 total aligned reads displayed.

GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGT  >  NC_000913/1809967‑1810164
                                                                                                   |                                                                                                   
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtc                                                                                                    >  1:1145029/1‑100 (MQ=255)
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtc                                                                                                    <  2:582664/100‑1 (MQ=255)
   aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCg                                                                                                 <  1:95500/100‑1 (MQ=255)
    aaTGATCCACTTTTTGTATTAGCAGACGTCAAATTTTCATCGCTGACAGAGGTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCg                                                                                                 <  1:639333/99‑1 (MQ=255)
    aaTGATCCACTTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGc                                                                                                <  2:2728239/100‑1 (MQ=255)
      tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACAGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCt                                                                                               >  2:1315799/1‑99 (MQ=255)
      tGATCCACGTTTTGTATTAGCAGACGACAAACTTTCATCGCTGACAGAACTCACTGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTg                                                                                              <  1:1170308/100‑1 (MQ=255)
       gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGc                                                                                             >  1:2461826/1‑100 (MQ=255)
        aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGcc                                                                                            >  1:1999632/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGcccc                                                                                          >  1:963338/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGcccc                                                                                          <  2:2477067/100‑1 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGcccc                                                                                          >  2:1945167/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGccc                                                                                           >  1:2603376/1‑99 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGcccc                                                                                          <  1:2462014/100‑1 (MQ=38)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGcccc                                                                                          <  1:2492888/100‑1 (MQ=38)
           cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                         >  1:1123023/1‑100 (MQ=255)
           cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCt                                                                                         <  2:1996979/100‑1 (MQ=38)
           cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCt                                                                                         <  1:410231/100‑1 (MQ=39)
           cACGTTCTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCt                                                                                         <  1:1449461/100‑1 (MQ=38)
            aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                        <  2:1251943/100‑1 (MQ=38)
             cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAAAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCTgg                                                                                       <  1:1550591/100‑1 (MQ=38)
              gTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAAAGCACCAGGGGCAGCGAGATCACTGcccc                                                                                          <  1:735025/97‑1 (MQ=37)
              gTTTTGTATTACCAGACGTCACACTTCCATCGCTGACAGAACTCACCGCAAGAGACCTTTTCGGTTAATG‑ACCAGGGCCAGCGAGACCACTGCCCCTGGt                                                                                      <  1:977256/100‑1 (MQ=255)
               ttttGTATTAGCATACGTCAAACGTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGt                                                                                      <  1:1732974/99‑1 (MQ=255)
                 ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCACTGGTTCt                                                                                   <  2:544688/100‑1 (MQ=255)
                 ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCTGGtcat                                                                                   <  1:1278621/100‑4 (MQ=37)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                  <  2:117644/100‑1 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                  >  2:64513/1‑100 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                  <  2:963341/100‑1 (MQ=255)
                   gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                                                 >  1:2815302/1‑100 (MQ=255)
                    tATTAGCCGACGTCAAGCTTTCATCGCTGACAGAACTCGCCGCAAAAGCCCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                                                <  1:152057/100‑1 (MQ=255)
                    tATTAGCAGACGTCAAACTTTCAGCGATGACAGATCTCACCGCAAAAGACCTTGCCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                                                 <  1:208866/99‑1 (MQ=255)
                       tAGCAGACGTCAAACGTTCATCGCTGACAGAACTCACCGCAAAATACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATc                                                                              <  1:2141683/99‑1 (MQ=255)
                        aGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTg                                                                            >  2:2602010/1‑100 (MQ=255)
                         gCAGACGTCAACCTTTTATCGCTGACAGACCTCACCGCAAAAGCCCGTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                           <  1:806531/100‑1 (MQ=255)
                          cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                           <  1:1052694/99‑1 (MQ=255)
                            gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                        <  2:1181526/100‑1 (MQ=255)
                            gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                        <  2:2168260/100‑1 (MQ=255)
                             aCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                                       <  2:1123028/100‑1 (MQ=255)
                               gTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                                                     >  1:333184/1‑100 (MQ=255)
                                tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCa                                                                    <  2:600333/100‑1 (MQ=255)
                                 cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                                                   >  1:52822/1‑100 (MQ=255)
                                  aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                                                  >  1:547306/1‑100 (MQ=255)
                                  aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                                                  >  1:2781529/1‑100 (MQ=255)
                                   aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                                                 <  2:720748/100‑1 (MQ=255)
                                    aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                                                <  2:28663/100‑1 (MQ=255)
                                    aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                                                <  2:1474282/100‑1 (MQ=255)
                                      tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                              <  2:1148974/100‑1 (MQ=255)
                                       ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACt                                                             >  1:940393/1‑100 (MQ=255)
                                         cATCGCTGACAGAACTCGCCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                                           >  2:815721/1‑100 (MQ=255)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                          >  2:2787891/1‑100 (MQ=255)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                          >  2:746945/1‑100 (MQ=255)
                                            cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACgg                                                        >  1:489109/1‑100 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                      >  2:939738/1‑100 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                      <  1:2302/100‑1 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGGGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                      <  2:906560/100‑1 (MQ=255)
                                                  cAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                                                   >  1:1809385/1‑99 (MQ=255)
                                                    gAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCg                                                >  1:1845217/1‑100 (MQ=255)
                                                      cgatctcGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                              <  2:718887/94‑1 (MQ=255)
                                                       cTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                             <  1:1262615/100‑1 (MQ=255)
                                                       cTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                             >  2:2613349/1‑100 (MQ=255)
                                                         cACCGCAAAAGACCTTCTCGATTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCagt                                           >  1:1773723/1‑98 (MQ=255)
                                                          aCCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                                          >  1:1732426/1‑100 (MQ=255)
                                                           ccGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                                         >  2:2430079/1‑100 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCagatcg                                        >  1:718884/1‑95 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGACCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAC‑GATTAATccgtcagttga                                       <  2:2552180/100‑11 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATG‑ACCAGAGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                                        <  1:725262/100‑1 (MQ=255)
                                                             gCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTaa                                       >  2:2656162/1‑100 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                      <  2:652761/100‑1 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                      <  1:2787873/100‑1 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                      <  2:1625618/100‑1 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                     <  2:175888/100‑1 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                     >  1:2885931/1‑100 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                     <  2:1159110/100‑1 (MQ=255)
                                                                       ttCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                               >  1:143816/1‑98 (MQ=255)
                                                                        tctcGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt                            >  2:2397999/1‑100 (MQ=255)
                                                                        tctcGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGta                            >  2:742296/1‑99 (MQ=255)
                                                                         ctcGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                           >  2:1084989/1‑100 (MQ=255)
                                                                          tcGGTTAATG‑ACCACGGGCAGTCATCCCGCTGCCCATGGATCGTTATCTGGAATAGCATGCAACGGACGGATAAACCGTCCATCCAAGAGAAACAGCTac                          >  1:1655172/1‑100 (MQ=255)
                                                                           cGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTaca                         >  1:2870539/1‑100 (MQ=255)
                                                                            ggTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                        >  1:1194337/1‑100 (MQ=255)
                                                                             gataaagaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                        >  1:733107/9‑100 (MQ=25)
                                                                             gataaagaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                        >  1:732710/9‑100 (MQ=25)
                                                                                aaTG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                    <  2:1620232/100‑1 (MQ=255)
                                                                                aaTG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                    >  2:1763216/1‑100 (MQ=255)
                                                                                 aTG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                   <  1:2656146/100‑1 (MQ=255)
                                                                                  tG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                  >  2:533109/1‑100 (MQ=255)
                                                                                    aaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacc                 >  2:1933640/2‑100 (MQ=38)
                                                                                     aCCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                <  2:1870610/100‑1 (MQ=255)
                                                                                     aCCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                >  1:1464129/1‑100 (MQ=255)
                                                                                     aCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                >  2:632534/1‑100 (MQ=38)
                                                                                      ccAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcaccac               <  1:2350248/100‑1 (MQ=255)
                                                                                      caaGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAAATGCATTCAACTGACGGATTAATCGTCGATTTAAGAGCCAGAGTTAAACCGTcaccac               >  1:940064/3‑100 (MQ=25)
                                                                                         ggggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTc            <  1:578241/100‑1 (MQ=255)
                                                                                         ggggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTc            >  2:2321662/1‑100 (MQ=255)
                                                                                         ggggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAAAAAGAGTTACACCGTCACCACTTc            >  1:2171857/1‑100 (MQ=255)
                                                                                             cAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTg        <  1:1013236/100‑1 (MQ=255)
                                                                                              aGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGc       <  2:842113/100‑1 (MQ=255)
                                                                                                 gATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACt    >  2:406672/1‑100 (MQ=255)
                                                                                                   tctcGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGt  >  2:1625892/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGT  >  NC_000913/1809967‑1810164

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