Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,810,066 C→A 15.2% intergenic (+18/‑145) 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,0660CA15.2% 313.9 / 34.4 112intergenic (+18/‑145)yniC/ydjM2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base C (47/48);  new base A (7/10);  total (54/58)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 119 total aligned reads displayed.

GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTA  >  NC_000913/1809967‑1810165
                                                                                                    |                                                                                                   
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtc                                                                                                     >  1:1145029/1‑100 (MQ=255)
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtc                                                                                                     <  2:582664/100‑1 (MQ=255)
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtc                                                                                                     >  2:2136115/1‑100 (MQ=255)
 cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtct                                                                                                    <  2:2587497/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)
     aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCt                                                                                                >  2:1053822/1‑100 (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                                                                                           <  2:2477067/100‑1 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGcccc                                                                                           >  1:963338/1‑100 (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)
           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‑ACCAGGGGCAGCGAGATCACTGCCCCTGGtcat                                                                                    <  1:1278621/100‑4 (MQ=37)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   <  2:117644/100‑1 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   <  2:963341/100‑1 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   >  2:64513/1‑100 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   <  1:1525461/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:2168260/100‑1 (MQ=255)
                            gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                         <  2:1181526/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:2781529/1‑100 (MQ=255)
                                  aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                                                   >  1:547306/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)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                           >  2:2787891/1‑100 (MQ=255)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                           >  2:746945/1‑100 (MQ=255)
                                           tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                          >  1:878665/1‑100 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                       <  1:2302/100‑1 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                       >  2:939738/1‑100 (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                                              >  2:2613349/1‑100 (MQ=255)
                                                       cTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                              <  1:1262615/100‑1 (MQ=255)
                                                         cACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                                            >  2:802880/1‑100 (MQ=255)
                                                          aCCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                                           >  1:1732426/1‑100 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCagatcg                                         >  1:718884/1‑95 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGACCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAC‑GATTAATccgtcagttga                                        <  2:2552180/100‑11 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                       <  2:301158/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                                      <  2:1159110/100‑1 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                      >  1:2885931/1‑100 (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‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtgt                            <  1:353219/100‑3 (MQ=255)
                                                                          tcGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                            <  1:695570/99‑1 (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:1763216/1‑100 (MQ=255)
                                                                                aaTG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                     <  2:1620232/100‑1 (MQ=255)
                                                                                 aTG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                    <  1:2656146/100‑1 (MQ=255)
                                                                                 aTG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                    <  1:18319/100‑1 (MQ=255)
                                                                                  tG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                   >  2:533109/1‑100 (MQ=255)
                                                                                    aaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacc                  >  2:1933640/2‑100 (MQ=38)
                                                                                     aCCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                 >  1:1464129/1‑100 (MQ=255)
                                                                                     aCCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                 <  2:1870610/100‑1 (MQ=255)
                                                                                     aCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                 >  2:632534/1‑100 (MQ=38)
                                                                                      ccAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcaccac                <  1:2350248/100‑1 (MQ=255)
                                                                                      ccAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTCACACCTTcacca                 >  1:234599/1‑99 (MQ=38)
                                                                                         ggggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTc             >  2:2321662/1‑100 (MQ=255)
                                                                                         ggggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTc             <  1:578241/100‑1 (MQ=255)
                                                                                          gttgattGATCTCGCTGTCTCTGGTGCTTTATCTGTAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca                 <  1:2762518/89‑1 (MQ=255)
                                                                                             cAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTg         <  1:1013236/100‑1 (MQ=255)
                                                                                              aGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGc        <  2:842113/100‑1 (MQ=255)
                                                                                                tGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAc      <  2:643638/100‑1 (MQ=255)
                                                                                                 gATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACt     >  2:406672/1‑100 (MQ=255)
                                                                                                   tctcGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGt   >  2:1625892/1‑100 (MQ=255)
                                                                                                    ctcGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGta  <  2:828783/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTA  >  NC_000913/1809967‑1810165

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: ‑