Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,810,062 T→C 15.1% intergenic (+14/‑149) 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,0620TC15.1% 333.8 / 49.2 148intergenic (+14/‑149)yniC/ydjM2‑deoxyglucose‑6‑P phosphatase/inner membrane protein regulated by LexA
Reads supporting (aligned to +/- strand):  ref base T (64/61);  new base C (12/10);  total (77/71)
Fisher's exact test for biased strand distribution p-value = 8.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01

Only 100 of 150 total aligned reads displayed.

AAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAC  >  NC_000913/1809963‑1810161
                                                                                                    |                                                                                                   
    gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtc                                                                                                 >  2:880128/1‑100 (MQ=255)
     cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtct                                                                                                >  1:1151039/1‑100 (MQ=255)
     cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGAtct                                                                                                >  1:94452/1‑100 (MQ=255)
      agaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAAACATTCTCGATTAATG‑ACCAGGGGCAGTGAtct                                                                                                >  1:895153/3‑99 (MQ=255)
          tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTg                                                                                           <  2:1544196/100‑1 (MQ=255)
          tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTg                                                                                           <  1:3232309/100‑1 (MQ=255)
           gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGc                                                                                          >  1:3678686/1‑100 (MQ=255)
           gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCACCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGc                                                                                          <  1:2749044/100‑1 (MQ=255)
             tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGccc                                                                                        >  1:3513650/1‑100 (MQ=255)
             tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGccc                                                                                        >  2:1807484/1‑100 (MQ=255)
              ccACGTTTTTTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGcccc                                                                                       <  2:2982852/100‑1 (MQ=255)
              ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGcccc                                                                                       <  1:1807499/100‑1 (MQ=255)
              ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAAGGGCAGTGATCTCGATGcccc                                                                                       >  1:2456162/1‑100 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑TCCAGTGGCAGTGATCTCGCTGcccc                                                                                       <  2:2445254/99‑1 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                      >  2:419207/1‑100 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                      >  1:3795131/1‑100 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAAAGAACCAGGGGCAGCGAGATCACTGcccc                                                                                       >  2:2116549/1‑100 (MQ=25)
                aCGTTTTGTATTAGCAGACGTCCCACTCTCATCGCTGATAGAACTCCCCGCAAAAGACGTTCTCGGTTACTG‑ACCAGGGGCAGCGCGATCACTGccccca                                                                                     <  2:2128178/100‑3 (MQ=255)
                aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTg                                                                                     >  2:2971694/1‑100 (MQ=255)
                aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                     <  1:1994608/100‑1 (MQ=38)
                 gatcttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTgn                                                                                    <  1:462391/96‑2 (MQ=255)
                 cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCTCACGACCTTCTCGGTTAATG‑ACCAGCGGCAGCGAGATCACTGCCCCTg                                                                                     <  2:784068/99‑1 (MQ=37)
                 cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTgg                                                                                    <  2:3679863/100‑1 (MQ=255)
                 cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAAAGAACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                     <  1:2116568/100‑1 (MQ=25)
                  tttttGTATTATCATACGTCAACCTGTTATCGTTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑TCCAGGGGCAGCGAGATCACTGCCCCt                                                                                      <  2:3198082/96‑1 (MQ=21)
                   ttttGTGTTAGAAGAAGTCAAACTTTCATCGCTGGCAGAACTCACCGCAAAAGTCCTTCTCGGTTATTG‑ACCAGGGTTAGCGATATCACTGCCCCTgg                                                                                    <  2:385174/98‑1 (MQ=255)
                   ttttGTATTAGCAGACTTCAAATTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGCGAGATCACTGCCCCTGGtc                                                                                  <  1:694566/100‑2 (MQ=37)
                      tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                               <  2:2016621/100‑1 (MQ=255)
                       gTATTAGCAGACGTCAAACTTTCATCGCTGCCAGTACTCACCTCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTTATCTCGCTGCCCCTGGTTCttt                                                                              <  2:2937649/100‑1 (MQ=255)
                        tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                                             <  1:880662/100‑1 (MQ=255)
                         aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                            <  2:3570472/100‑1 (MQ=255)
                         aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                            >  2:3790289/1‑100 (MQ=255)
                             gCAGACGTCAAACTTTAATCGATAAAAGAAAACACAGCAAAAGACATTATCGATTAATG‑ACCAAGGGAAGAGATTTAGATGCTCATGATTATTAATCt                                                                          >  2:922547/1‑98 (MQ=255)
                              cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAATACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                        <  1:2681504/99‑1 (MQ=255)
                                gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                     >  2:1012012/1‑100 (MQ=255)
                                 aCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTCTCTGAAAtt                                                                    >  2:1884039/1‑100 (MQ=255)
                                  cGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTg                                                                   <  1:1678540/100‑1 (MQ=255)
                                     cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTTTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                                                <  1:2365878/100‑1 (MQ=255)
                                     cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGGTCTTTATCTGAAATTGCa                                                                 <  2:1465653/99‑1 (MQ=255)
                                        aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                                             <  2:1813643/100‑1 (MQ=255)
                                         cTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCaa                                                            >  2:477945/1‑100 (MQ=255)
                                         cTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCACTGGTTCTTTATCTGAAATTGCATTCaa                                                            >  1:1938477/1‑100 (MQ=255)
                                          tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                           <  2:1603564/100‑1 (MQ=255)
                                          tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGATTAATG‑ACAAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                           >  2:2516126/1‑100 (MQ=255)
                                           ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACt                                                          >  1:2261477/1‑100 (MQ=255)
                                             cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                                        >  2:512864/1‑100 (MQ=255)
                                             cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                                        >  1:2487022/1‑100 (MQ=255)
                                             cATCGCTGACAGAACTCACCGCAAAAGACCTTCCCGGTTAATG‑ACCAGGGCCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGACATTACATTCAACtca                                                        >  1:2894341/1‑98 (MQ=255)
                                                    gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTa                                                 >  2:267640/1‑100 (MQ=255)
                                                      cAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                                               <  2:1684455/100‑1 (MQ=255)
                                                       aGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                              <  1:3429426/100‑1 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                           >  1:1744776/1‑99 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                          >  2:2871894/1‑100 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                          >  2:1769243/1‑100 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                          >  2:1101323/1‑100 (MQ=255)
                                                              aCCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCTTGGTTCTTTATCTGAAATTGCATTCAACTGAC‑GATTAATccgtcagtt                                      <  1:107519/100‑9 (MQ=255)
                                                                cGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                                     <  2:3678647/100‑1 (MQ=255)
                                                                  cAAAAGACCGTGTCGGTTCATG‑GCTAGGGGCAGTGATCTCGCTGCTCTTGTTTCTTTATCTGAAATTGCATTCAACTGATGGATTAATCGTCAAttt                                      <  2:3386603/97‑1 (MQ=255)
                                                                   aaaaGTCCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                  <  1:2772907/100‑1 (MQ=255)
                                                                   aaaaGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                  >  2:2718489/1‑100 (MQ=255)
                                                                    aaaGACCTTCTCTGTTAATT‑ACCAGGGGCAGTTATCTCGCTTCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                  <  1:2275324/99‑1 (MQ=255)
                                                                    aaaGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag                                 <  1:2321263/100‑1 (MQ=255)
                                                                         ccTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTTAAAATGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAag                             <  2:2609550/99‑1 (MQ=255)
                                                                         ccTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                            >  2:400819/1‑100 (MQ=255)
                                                                         ccTTCTCGGTTAATG‑ACCAGGGGAAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGAATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                            >  2:937969/1‑100 (MQ=255)
                                                                          cTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                           <  1:2987234/100‑1 (MQ=255)
                                                                          cTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                           <  1:2742291/100‑1 (MQ=255)
                                                                             ctcGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                        >  2:3488415/1‑100 (MQ=255)
                                                                                ggTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                     >  2:821692/1‑100 (MQ=255)
                                                                                ggTTAATG‑ACCAGGGACAGTGATCTCGCTGCCCATGGTTCTTTATCTGAAATTGCATTCAACTGACGGATAACTCGTCAATTTAAGAGAAAGAGTTaca                      >  2:2353649/1‑99 (MQ=255)
                                                                                 tttAATT‑ATCAGGTGCACTGATCTCTCTGCCCCTTGTTCTTTATCTGAAATTACATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                     <  1:937980/98‑1 (MQ=255)
                                                                                 gTTAATG‑ACCAGGGGCAGTGATCTCGCTGTCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAGCGTCAATTTAAGAGAAAGAGTTACAcc                    <  1:1354389/100‑1 (MQ=255)
                                                                                 gTTAATG‑ACCAGGGGCAGTGATCTCGCTGCTCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAcc                    <  1:1287080/100‑1 (MQ=255)
                                                                                 gTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAcc                    >  2:3059756/1‑100 (MQ=255)
                                                                                  ttAATG‑ACCAGTGTCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCg                   <  2:3298064/100‑1 (MQ=255)
                                                                                     aagaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                 >  2:2025710/5‑100 (MQ=37)
                                                                                     aTG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                >  2:3115011/1‑100 (MQ=255)
                                                                                      tG‑ACCAGGGGCAGTGATCTCGCTACCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac               >  2:1554809/1‑100 (MQ=255)
                                                                                       g‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacc              <  2:3746140/100‑1 (MQ=255)
                                                                                       g‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacc              <  2:3196818/100‑1 (MQ=255)
                                                                                         aCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca             >  1:3280093/1‑100 (MQ=38)
                                                                                         aCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTAAATTTAAGAGAAAGAGTTACACCGTcacca             >  1:2681320/1‑100 (MQ=37)
                                                                                          ccAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcaccac            >  2:3855939/1‑100 (MQ=255)
                                                                                          ccAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcaccac            >  2:1238407/1‑100 (MQ=255)
                                                                                          ccAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTAAATTTAAAAGAAAGAGTTACACCGTcaccac            >  1:334118/1‑100 (MQ=38)
                                                                                           cAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATTTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacca             <  2:966981/98‑1 (MQ=255)
                                                                                           cAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt           <  2:2492775/100‑1 (MQ=255)
                                                                                           cAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt           <  1:2773330/100‑1 (MQ=255)
                                                                                           cAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGTATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt           >  1:1812173/1‑100 (MQ=38)
                                                                                           cAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt           >  1:828104/1‑100 (MQ=38)
                                                                                           cAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt           >  1:2960072/1‑100 (MQ=38)
                                                                                           cAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt           >  1:2867553/1‑100 (MQ=38)
                                                                                           cAGGGGCAGCGAGATCACTGCCCCTGGATCTATATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACt           >  2:178787/1‑100 (MQ=38)
                                                                                             ggggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTc         <  1:2075586/100‑1 (MQ=255)
                                                                                               ggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCg       >  2:1869653/1‑100 (MQ=255)
                                                                                               ggCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCg       <  1:813499/100‑1 (MQ=255)
                                                                                                 cAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTg     <  1:3469098/100‑1 (MQ=255)
                                                                                                  aGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGc    <  2:2780207/100‑1 (MQ=255)
                                                                                                   gTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCa   <  2:1447097/100‑1 (MQ=255)
                                                                                                   gTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCa   >  2:1325082/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
AAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATG‑ACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAC  >  NC_000913/1809963‑1810161

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