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.

GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGT  >  NC_000913/1809967‑1810164
                                                                                                   |                                                                                                   
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                    >  1:1145029/1‑100 (MQ=255)
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                    <  2:582664/100‑1 (MQ=255)
   aaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCg                                                                                                 <  1:95500/100‑1 (MQ=255)
    aaTGATCCACTTTTTGTATTAGCAGACGTCAAATTTTCATCGCTGACAGAGGTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCg                                                                                                 <  1:639333/99‑1 (MQ=255)
    aaTGATCCACTTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGc                                                                                                <  2:2728239/100‑1 (MQ=255)
      tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACAGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCt                                                                                               >  2:1315799/1‑99 (MQ=255)
      tGATCCACGTTTTGTATTAGCAGACGACAAACTTTCATCGCTGACAGAACTCACTGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTg                                                                                              <  1:1170308/100‑1 (MQ=255)
       gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGc                                                                                             >  1:2461826/1‑100 (MQ=255)
        aTCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcc                                                                                            >  1:1999632/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                          >  1:963338/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                          <  2:2477067/100‑1 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                          >  2:1945167/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                           >  1:2603376/1‑99 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                          <  1:2462014/100‑1 (MQ=38)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                          <  1:2492888/100‑1 (MQ=38)
           cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                         >  1:1123023/1‑100 (MQ=255)
           cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCt                                                                                         <  2:1996979/100‑1 (MQ=38)
           cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCt                                                                                         <  1:410231/100‑1 (MQ=39)
           cACGTTCTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCt                                                                                         <  1:1449461/100‑1 (MQ=38)
            aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                        <  2:1251943/100‑1 (MQ=38)
             cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAAAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTgg                                                                                       <  1:1550591/100‑1 (MQ=38)
              gTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAAAGCACCAGGGGCAGCGAGATCACTGcccc                                                                                          <  1:735025/97‑1 (MQ=37)
              gTTTTGTATTACCAGACGTCACACTTCCATCGCTGACAGAACTCACCGCAAGAGACCTTTTCGGTTAATGACCAGGGCCAGCGAGACCACTGCCCCTGGt                                                                                      <  1:977256/100‑1 (MQ=255)
               ttttGTATTAGCATACGTCAAACGTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGt                                                                                      <  1:1732974/99‑1 (MQ=255)
                 ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCACTGGTTCt                                                                                   <  2:544688/100‑1 (MQ=255)
                 ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTGGtcat                                                                                   <  1:1278621/100‑4 (MQ=37)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                  <  2:117644/100‑1 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                  >  2:64513/1‑100 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                  <  2:963341/100‑1 (MQ=255)
                   gTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                                                 >  1:2815302/1‑100 (MQ=255)
                    tATTAGCCGACGTCAAGCTTTCATCGCTGACAGAACTCGCCGCAAAAGCCCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                                                <  1:152057/100‑1 (MQ=255)
                    tATTAGCAGACGTCAAACTTTCAGCGATGACAGATCTCACCGCAAAAGACCTTGCCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCttt                                                                                 <  1:208866/99‑1 (MQ=255)
                       tAGCAGACGTCAAACGTTCATCGCTGACAGAACTCACCGCAAAATACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATc                                                                              <  1:2141683/99‑1 (MQ=255)
                        aGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTg                                                                            >  2:2602010/1‑100 (MQ=255)
                         gCAGACGTCAACCTTTTATCGCTGACAGACCTCACCGCAAAAGCCCGTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                           <  1:806531/100‑1 (MQ=255)
                          cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                           <  1:1052694/99‑1 (MQ=255)
                            gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                        <  2:1181526/100‑1 (MQ=255)
                            gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                        <  2:2168260/100‑1 (MQ=255)
                             aCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAtt                                                                       <  2:1123028/100‑1 (MQ=255)
                               gTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGc                                                                     >  1:333184/1‑100 (MQ=255)
                                tCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCa                                                                    <  2:600333/100‑1 (MQ=255)
                                 cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                                                   >  1:52822/1‑100 (MQ=255)
                                  aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                                                  >  1:547306/1‑100 (MQ=255)
                                  aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                                                  >  1:2781529/1‑100 (MQ=255)
                                   aaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTc                                                                 <  2:720748/100‑1 (MQ=255)
                                    aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                                                <  2:28663/100‑1 (MQ=255)
                                    aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                                                <  2:1474282/100‑1 (MQ=255)
                                      tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                              <  2:1148974/100‑1 (MQ=255)
                                       ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACt                                                             >  1:940393/1‑100 (MQ=255)
                                         cATCGCTGACAGAACTCGCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                                           >  2:815721/1‑100 (MQ=255)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                          >  2:2787891/1‑100 (MQ=255)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                          >  2:746945/1‑100 (MQ=255)
                                            cGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACgg                                                        >  1:489109/1‑100 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                      >  2:939738/1‑100 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                      <  1:2302/100‑1 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGGGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                      <  2:906560/100‑1 (MQ=255)
                                                  cAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTaa                                                   >  1:1809385/1‑99 (MQ=255)
                                                    gAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCg                                                >  1:1845217/1‑100 (MQ=255)
                                                      cgatctcGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                              <  2:718887/94‑1 (MQ=255)
                                                       cTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                             <  1:1262615/100‑1 (MQ=255)
                                                       cTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                             >  2:2613349/1‑100 (MQ=255)
                                                         cACCGCAAAAGACCTTCTCGATTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCagt                                           >  1:1773723/1‑98 (MQ=255)
                                                          aCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                                          >  1:1732426/1‑100 (MQ=255)
                                                           ccGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAttt                                         >  2:2430079/1‑100 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCagatcg                                        >  1:718884/1‑95 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGACCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGATTAATccgtcagttga                                       <  2:2552180/100‑11 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATGACCAGAGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                                        <  1:725262/100‑1 (MQ=255)
                                                             gCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTaa                                       >  2:2656162/1‑100 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                      <  2:652761/100‑1 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                      <  1:2787873/100‑1 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                      <  2:1625618/100‑1 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                     <  2:175888/100‑1 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                     >  1:2885931/1‑100 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                     <  2:1159110/100‑1 (MQ=255)
                                                                       ttCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                               >  1:143816/1‑98 (MQ=255)
                                                                        tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt                            >  2:2397999/1‑100 (MQ=255)
                                                                        tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGta                            >  2:742296/1‑99 (MQ=255)
                                                                         ctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                           >  2:1084989/1‑100 (MQ=255)
                                                                          tcGGTTAATGACCACGGGCAGTCATCCCGCTGCCCATGGATCGTTATCTGGAATAGCATGCAACGGACGGATAAACCGTCCATCCAAGAGAAACAGCTac                          >  1:1655172/1‑100 (MQ=255)
                                                                           cGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTaca                         >  1:2870539/1‑100 (MQ=255)
                                                                            ggTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                        >  1:1194337/1‑100 (MQ=255)
                                                                             gataaagaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                        >  1:733107/9‑100 (MQ=25)
                                                                             gataaagaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                        >  1:732710/9‑100 (MQ=25)
                                                                                aaTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                    <  2:1620232/100‑1 (MQ=255)
                                                                                aaTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                    >  2:1763216/1‑100 (MQ=255)
                                                                                 aTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                   <  1:2656146/100‑1 (MQ=255)
                                                                                  tGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                  >  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)
                                                                                                   |                                                                                                   
GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGT  >  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: