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.

AAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAC  >  NC_000913/1809963‑1810161
                                                                                                    |                                                                                                   
    gcAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtc                                                                                                 >  2:880128/1‑100 (MQ=255)
     cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtct                                                                                                >  1:1151039/1‑100 (MQ=255)
     cAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGAtct                                                                                                >  1:94452/1‑100 (MQ=255)
      agaaTGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAAACATTCTCGATTAATGACCAGGGGCAGTGAtct                                                                                                >  1:895153/3‑99 (MQ=255)
          tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTg                                                                                           <  2:1544196/100‑1 (MQ=255)
          tGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTg                                                                                           <  1:3232309/100‑1 (MQ=255)
           gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGc                                                                                          >  1:3678686/1‑100 (MQ=255)
           gATCCACGTTTTGTATTAGCAGACGTCAAACTTTCACCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGc                                                                                          <  1:2749044/100‑1 (MQ=255)
             tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                        >  1:3513650/1‑100 (MQ=255)
             tCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGccc                                                                                        >  2:1807484/1‑100 (MQ=255)
              ccACGTTTTTTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                       <  2:2982852/100‑1 (MQ=255)
              ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGcccc                                                                                       <  1:1807499/100‑1 (MQ=255)
              ccACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAAGGGCAGTGATCTCGATGcccc                                                                                       >  1:2456162/1‑100 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGTCCAGTGGCAGTGATCTCGCTGcccc                                                                                       <  2:2445254/99‑1 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                      >  2:419207/1‑100 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCt                                                                                      >  1:3795131/1‑100 (MQ=255)
               cACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAAAGAACCAGGGGCAGCGAGATCACTGcccc                                                                                       >  2:2116549/1‑100 (MQ=25)
                aCGTTTTGTATTAGCAGACGTCCCACTCTCATCGCTGATAGAACTCCCCGCAAAAGACGTTCTCGGTTACTGACCAGGGGCAGCGCGATCACTGccccca                                                                                     <  2:2128178/100‑3 (MQ=255)
                aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTg                                                                                     >  2:2971694/1‑100 (MQ=255)
                aCGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                     <  1:1994608/100‑1 (MQ=38)
                 gatcttGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTgn                                                                                    <  1:462391/96‑2 (MQ=255)
                 cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCTCACGACCTTCTCGGTTAATGACCAGCGGCAGCGAGATCACTGCCCCTg                                                                                     <  2:784068/99‑1 (MQ=37)
                 cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTgg                                                                                    <  2:3679863/100‑1 (MQ=255)
                 cGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAAAGAACCAGGGGCAGCGAGATCACTGCCCCTg                                                                                     <  1:2116568/100‑1 (MQ=25)
                  tttttGTATTATCATACGTCAACCTGTTATCGTTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGTCCAGGGGCAGCGAGATCACTGCCCCt                                                                                      <  2:3198082/96‑1 (MQ=21)
                   ttttGTGTTAGAAGAAGTCAAACTTTCATCGCTGGCAGAACTCACCGCAAAAGTCCTTCTCGGTTATTGACCAGGGTTAGCGATATCACTGCCCCTgg                                                                                    <  2:385174/98‑1 (MQ=255)
                   ttttGTATTAGCAGACTTCAAATTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGCGAGATCACTGCCCCTGGtc                                                                                  <  1:694566/100‑2 (MQ=37)
                      tGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCtt                                                                               <  2:2016621/100‑1 (MQ=255)
                       gTATTAGCAGACGTCAAACTTTCATCGCTGCCAGTACTCACCTCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTTATCTCGCTGCCCCTGGTTCttt                                                                              <  2:2937649/100‑1 (MQ=255)
                        tATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTa                                                                             <  1:880662/100‑1 (MQ=255)
                         aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                            <  2:3570472/100‑1 (MQ=255)
                         aTTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTAt                                                                            >  2:3790289/1‑100 (MQ=255)
                             gCAGACGTCAAACTTTAATCGATAAAAGAAAACACAGCAAAAGACATTATCGATTAATGACCAAGGGAAGAGATTTAGATGCTCATGATTATTAATCt                                                                          >  2:922547/1‑98 (MQ=255)
                              cAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAATACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGa                                                                        <  1:2681504/99‑1 (MQ=255)
                                gACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAAt                                                                     >  2:1012012/1‑100 (MQ=255)
                                 aCGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTCTCTGAAAtt                                                                    >  2:1884039/1‑100 (MQ=255)
                                  cGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTg                                                                   <  1:1678540/100‑1 (MQ=255)
                                     cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTTTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCAt                                                                <  1:2365878/100‑1 (MQ=255)
                                     cAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGGTCTTTATCTGAAATTGCa                                                                 <  2:1465653/99‑1 (MQ=255)
                                        aCTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCa                                                             <  2:1813643/100‑1 (MQ=255)
                                         cTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCaa                                                            >  2:477945/1‑100 (MQ=255)
                                         cTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCACTGGTTCTTTATCTGAAATTGCATTCaa                                                            >  1:1938477/1‑100 (MQ=255)
                                          tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                           <  2:1603564/100‑1 (MQ=255)
                                          tttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGATTAATGACAAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAAc                                                           >  2:2516126/1‑100 (MQ=255)
                                           ttCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACt                                                          >  1:2261477/1‑100 (MQ=255)
                                             cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                                        >  2:512864/1‑100 (MQ=255)
                                             cATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGa                                                        >  1:2487022/1‑100 (MQ=255)
                                             cATCGCTGACAGAACTCACCGCAAAAGACCTTCCCGGTTAATGACCAGGGCCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGACATTACATTCAACtca                                                        >  1:2894341/1‑98 (MQ=255)
                                                    gACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTa                                                 >  2:267640/1‑100 (MQ=255)
                                                      cAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAt                                               <  2:1684455/100‑1 (MQ=255)
                                                       aGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATc                                              <  1:3429426/100‑1 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTc                                           >  1:1744776/1‑99 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                          >  2:2871894/1‑100 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                          >  2:1769243/1‑100 (MQ=255)
                                                           cTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCa                                          >  2:1101323/1‑100 (MQ=255)
                                                              aCCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCTTGGTTCTTTATCTGAAATTGCATTCAACTGACGATTAATccgtcagtt                                      <  1:107519/100‑9 (MQ=255)
                                                                cGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTa                                     <  2:3678647/100‑1 (MQ=255)
                                                                  cAAAAGACCGTGTCGGTTCATGGCTAGGGGCAGTGATCTCGCTGCTCTTGTTTCTTTATCTGAAATTGCATTCAACTGATGGATTAATCGTCAAttt                                      <  2:3386603/97‑1 (MQ=255)
                                                                   aaaaGTCCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                  <  1:2772907/100‑1 (MQ=255)
                                                                   aaaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                  >  2:2718489/1‑100 (MQ=255)
                                                                    aaaGACCTTCTCTGTTAATTACCAGGGGCAGTTATCTCGCTTCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAaga                                  <  1:2275324/99‑1 (MQ=255)
                                                                    aaaGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAagag                                 <  1:2321263/100‑1 (MQ=255)
                                                                         ccTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTTAAAATGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAag                             <  2:2609550/99‑1 (MQ=255)
                                                                         ccTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                            >  2:400819/1‑100 (MQ=255)
                                                                         ccTTCTCGGTTAATGACCAGGGGAAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGAATTCAACTGACGGATTAATCGTCAATTTAAGAGAAaga                            >  2:937969/1‑100 (MQ=255)
                                                                          cTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                           <  1:2987234/100‑1 (MQ=255)
                                                                          cTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAagag                           <  1:2742291/100‑1 (MQ=255)
                                                                             ctcGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTa                        >  2:3488415/1‑100 (MQ=255)
                                                                                ggTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                     >  2:821692/1‑100 (MQ=255)
                                                                                ggTTAATGACCAGGGACAGTGATCTCGCTGCCCATGGTTCTTTATCTGAAATTGCATTCAACTGACGGATAACTCGTCAATTTAAGAGAAAGAGTTaca                      >  2:2353649/1‑99 (MQ=255)
                                                                                 tttAATTATCAGGTGCACTGATCTCTCTGCCCCTTGTTCTTTATCTGAAATTACATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTacac                     <  1:937980/98‑1 (MQ=255)
                                                                                 gTTAATGACCAGGGGCAGTGATCTCGCTGTCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAAGCGTCAATTTAAGAGAAAGAGTTACAcc                    <  1:1354389/100‑1 (MQ=255)
                                                                                 gTTAATGACCAGGGGCAGTGATCTCGCTGCTCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAcc                    <  1:1287080/100‑1 (MQ=255)
                                                                                 gTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACAcc                    >  2:3059756/1‑100 (MQ=255)
                                                                                  ttAATGACCAGTGTCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCg                   <  2:3298064/100‑1 (MQ=255)
                                                                                     aagaaCCAGGGGCAGCGAGATCACTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTc                 >  2:2025710/5‑100 (MQ=37)
                                                                                     aTGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTca                >  2:3115011/1‑100 (MQ=255)
                                                                                      tGACCAGGGGCAGTGATCTCGCTACCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcac               >  2:1554809/1‑100 (MQ=255)
                                                                                       gACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacc              <  2:3746140/100‑1 (MQ=255)
                                                                                       gACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTcacc              <  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)
                                                                                                    |                                                                                                   
AAGCGCAAAATGATCCACGTTTTGTATTAGCAGACGTCAAACTTTCATCGCTGACAGAACTCACCGCAAAAGACCTTCTCGGTTAATGACCAGGGGCAGTGATCTCGCTGCCCCTGGTTCTTTATCTGAAATTGCATTCAACTGACGGATTAATCGTCAATTTAAGAGAAAGAGTTACACCGTCACCACTTCCGTGCAC  >  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: