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.

GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTA  >  NC_000913/1809967‑1810165
                                                                                                    |                                                                                                   
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                     >  1:1145029/1‑100 (MQ=255)
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                     <  2:582664/100‑1 (MQ=255)
gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                     >  2:2136115/1‑100 (MQ=255)
 cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtct                                                                                                    <  2:2587497/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)
     aTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCt                                                                                                >  2:1053822/1‑100 (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                                                                                           <  2:2477067/100‑1 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                           >  1:963338/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                           >  2:1945167/1‑100 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                            >  1:2603376/1‑99 (MQ=255)
          ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGcccc                                                                                           <  1:2462014/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)
                 ttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTGGtcat                                                                                    <  1:1278621/100‑4 (MQ=37)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   <  2:117644/100‑1 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   <  2:963341/100‑1 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   >  2:64513/1‑100 (MQ=255)
                  tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                                   <  1:1525461/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:2168260/100‑1 (MQ=255)
                            gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                         <  2:1181526/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:2781529/1‑100 (MQ=255)
                                  aaaCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAtt                                                                   >  1:547306/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)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                           >  2:2787891/1‑100 (MQ=255)
                                          aTCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGAc                                                           >  2:746945/1‑100 (MQ=255)
                                           tCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACg                                                          >  1:878665/1‑100 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                       <  1:2302/100‑1 (MQ=255)
                                              cTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGAt                                                       >  2:939738/1‑100 (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                                              >  2:2613349/1‑100 (MQ=255)
                                                       cTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                              <  1:1262615/100‑1 (MQ=255)
                                                         cACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAt                                            >  2:802880/1‑100 (MQ=255)
                                                          aCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAAtt                                           >  1:1732426/1‑100 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCagatcg                                         >  1:718884/1‑95 (MQ=255)
                                                            cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGACCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGATTAATccgtcagttga                                        <  2:2552180/100‑11 (MQ=255)
                                                              cAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAag                                       <  2:301158/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                                      <  2:1159110/100‑1 (MQ=255)
                                                               aaaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                      >  1:2885931/1‑100 (MQ=255)
                                                                       ttCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                                >  1:143816/1‑98 (MQ=255)
                                                                        tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtt                             >  2:2397999/1‑100 (MQ=255)
                                                                        tctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGta                             >  2:742296/1‑99 (MQ=255)
                                                                         ctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGtgt                            <  1:353219/100‑3 (MQ=255)
                                                                          tcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                            <  1:695570/99‑1 (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:1763216/1‑100 (MQ=255)
                                                                                aaTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                     <  2:1620232/100‑1 (MQ=255)
                                                                                 aTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                    <  1:2656146/100‑1 (MQ=255)
                                                                                 aTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                    <  1:18319/100‑1 (MQ=255)
                                                                                  tGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac                   >  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)
                                                                                                    |                                                                                                   
GCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCACTGTA  >  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: