Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,261 G→T 17.1% intergenic (+31/‑186) typA → / → yihL GTP‑binding protein/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,060,2610GT17.1% 574.8 / 100.0 222intergenic (+31/‑186)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (78/106);  new base T (21/17);  total (99/123)
Fisher's exact test for biased strand distribution p-value = 1.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.30e-01

Only 100 of 224 total aligned reads displayed.

GTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATT  >  NC_000913/4060163‑4060359
                                                                                                  |                                                                                                   
  aTTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGa                                                                                                  <  1:3305029/100‑1 (MQ=255)
     cGTAAACGTCACCTGATGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                               >  2:1639788/1‑100 (MQ=255)
     cGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACAGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATgg                                                                                               >  2:548995/1‑100 (MQ=255)
      gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                              >  2:1010892/1‑100 (MQ=255)
      gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                              <  1:3603348/100‑1 (MQ=255)
      gTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGc                                                                                              >  1:3754659/1‑100 (MQ=255)
        aaaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAg                                                                                            >  1:2334996/1‑100 (MQ=255)
         aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAgg                                                                                           <  1:51291/100‑1 (MQ=255)
         aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                           <  2:1470378/100‑1 (MQ=255)
         aaCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAgg                                                                                           <  1:697563/100‑1 (MQ=255)
          aCGTCACCTGACGGAAATCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTAAAAAAAAAAAAAAAACCTGCCATCGCTGGCAGGt                                                                                         <  2:2760008/100‑1 (MQ=25)
          aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCTATCGCTGGCAgtt                                                                                          <  1:2569307/100‑3 (MQ=255)
          aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                         <  1:2822258/100‑1 (MQ=255)
          aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGt                                                                                         <  2:1584690/100‑1 (MQ=255)
           cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                        <  1:345862/100‑1 (MQ=255)
           cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                        >  1:2128866/1‑100 (MQ=255)
            gTCACCTGACGGAAAACGATCGTCGCCGCGCTAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttt                                                                                       <  2:486112/100‑1 (MQ=255)
            gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttt                                                                                       <  1:197877/100‑1 (MQ=255)
             tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                      <  2:2773564/100‑1 (MQ=255)
                ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                   <  1:3437051/100‑2 (MQ=255)
                    tcGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                               <  1:1409532/99‑1 (MQ=255)
                    aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgac                                                                               <  1:3659308/100‑6 (MQ=38)
                    aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatgac                                                                               <  1:1548713/100‑6 (MQ=38)
                     cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                              >  2:2244178/1‑100 (MQ=255)
                     cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGcc                                                                              >  2:737418/1‑100 (MQ=255)
                     cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGca                                                                              >  2:1409522/1‑99 (MQ=255)
                      ggAAAACGATTGTCGCCGCGCCAACCGCGCAGCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTg                                                                                <  2:1682913/97‑1 (MQ=255)
                        aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGt                                                                           >  1:38925/1‑100 (MQ=255)
                         aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                                                                          >  1:2453127/1‑100 (MQ=255)
                         aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTc                                                                          >  1:3134739/1‑100 (MQ=255)
                            cGATCGTCGCCGCGCCACCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCg                                                                         <  2:437023/98‑1 (MQ=255)
                            cGATCGTCGCCACGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCCGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTg                                                                       >  2:2002298/1‑100 (MQ=255)
                             gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAACTTCGTCTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGCTTTTTTTTGCCTGTCCTGa                                                                      >  2:868946/1‑100 (MQ=255)
                               tcgtcgCCGCGCCAACCGTGCACCGAAAGTCGATTAATTTCTTTTGTTAGTCATAAAAAGCCTGCCAGCGATGGCAGGTTTTGCTTTGCCTGTCGTGAg                                                                     <  2:711944/99‑1 (MQ=255)
                                 gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAgcgc                                                                  <  1:607821/100‑1 (MQ=255)
                                    gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGcac                                                                >  1:1757371/1‑100 (MQ=255)
                                    gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGcaca                                                               >  1:2797679/1‑100 (MQ=255)
                                      cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                             <  1:917241/100‑1 (MQ=255)
                                      cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                             <  2:2943463/100‑1 (MQ=255)
                                          ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCTTGTGTTGGGCGCccacca                                                            >  1:3725681/1‑98 (MQ=255)
                                           cAACCGCGCACCGAAGGCCGATTAATGTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                             <  2:3125217/95‑1 (MQ=255)
                                           cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                        >  1:1908412/1‑100 (MQ=255)
                                              ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAGCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     <  2:1863153/100‑1 (MQ=255)
                                              ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     <  2:557929/100‑1 (MQ=255)
                                              ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     >  1:2631400/1‑100 (MQ=255)
                                               cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     <  2:2902084/99‑1 (MQ=255)
                                                gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                   <  2:2186783/100‑1 (MQ=255)
                                                gcgcACCGAAAGACAATTAATTTCGTTTTTTTGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCCGTCGTGAGCGCACACCAATGCGAGTa                                                   <  2:352381/100‑1 (MQ=255)
                                                 cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                   <  1:1066086/100‑1 (MQ=255)
                                                   cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                >  1:362014/1‑100 (MQ=255)
                                                      cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                             <  1:661379/100‑1 (MQ=255)
                                                      cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                             >  1:3233913/1‑100 (MQ=255)
                                                       gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                             >  1:3489786/1‑99 (MQ=255)
                                                         atgACGTTTATTTCCGTTTTTTAGTCATAAAAAACTTCCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                          <  2:2320182/98‑1 (MQ=255)
                                                         aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                          <  1:2076650/100‑1 (MQ=255)
                                                          aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                         <  2:3489759/100‑1 (MQ=255)
                                                          aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                         >  2:2856310/1‑100 (MQ=255)
                                                           gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                        <  2:1536403/100‑1 (MQ=255)
                                                            aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                       <  1:73847/100‑1 (MQ=255)
                                                             cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                      <  2:589923/100‑1 (MQ=255)
                                                             cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                      <  2:3241198/100‑1 (MQ=255)
                                                              gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGAACACCAATGCGAGTATCTCCCTTATAAg                                      >  2:41741/1‑99 (MQ=255)
                                                               aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                     <  2:2355118/100‑1 (MQ=255)
                                                               aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                    >  1:1250399/1‑100 (MQ=255)
                                                                 tAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                   <  2:1951423/100‑1 (MQ=255)
                                                                  aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                 <  2:1191702/100‑1 (MQ=255)
                                                                    tttCGTTTTTTAGTCATAAAAGACCTGCCGGCGATGGCAGGTTTTTTTTTGCCTGTCGTGGGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                               <  2:1351216/100‑1 (MQ=255)
                                                                    tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                               <  2:1328613/100‑1 (MQ=255)
                                                                     ttCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                              <  2:1652460/100‑1 (MQ=255)
                                                                     ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTACCTGTCGTGAGCGCACACCAATGCGTGTATCTCCCTTATATGTCTGTGGt                              >  2:767689/1‑100 (MQ=255)
                                                                        gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                           >  2:1755803/1‑100 (MQ=255)
                                                                         ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                          <  2:2629996/100‑1 (MQ=255)
                                                                               aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                    <  2:1599363/100‑1 (MQ=255)
                                                                                gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                   >  1:1829791/6‑100 (MQ=38)
                                                                                gcaaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                   >  1:2389119/6‑100 (MQ=38)
                                                                                 caaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                  >  1:2414999/5‑100 (MQ=255)
                                                                                  cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                 >  1:597492/1‑100 (MQ=255)
                                                                                  aaaaaaaaaCCTGCCATCACTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                  >  1:2326262/4‑99 (MQ=38)
                                                                                   aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                >  1:3073154/3‑100 (MQ=255)
                                                                                   aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                <  2:1523476/100‑1 (MQ=255)
                                                                                    aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag               >  1:1210445/2‑100 (MQ=255)
                                                                                     aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc              >  1:2207004/1‑100 (MQ=255)
                                                                                     aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc              >  1:2252024/1‑100 (MQ=255)
                                                                                       aaaaCCTGCCATCGCTGGCAGGTTTTTTTTTTTTTTTCTTTATCGCACACCAATGCGCGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt            >  2:3447288/1‑100 (MQ=25)
                                                                                       aaaaCCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt            >  1:3424982/1‑100 (MQ=255)
                                                                                        aaaCCTGCCAGCGATGGCAGGTTTTTTTTTTTTTGTTTTTATCGCACACCAATGCGCGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt           >  2:2546819/1‑100 (MQ=255)
                                                                                         aaCCTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg          >  2:3815271/1‑100 (MQ=255)
                                                                                         aaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt         <  1:103045/100‑1 (MQ=255)
                                                                                          aCCTGCCAGCGATGGCAGGTTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt         >  1:1662847/1‑100 (MQ=255)
                                                                                          aCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt         <  1:878991/99‑1 (MQ=255)
                                                                                           ccTGCCATCGCTGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc        >  1:3195262/1‑100 (MQ=255)
                                                                                           ccTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc        <  2:1217528/100‑1 (MQ=255)
                                                                                            cTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca       <  2:3603964/100‑1 (MQ=255)
                                                                                            cTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca       <  1:3468588/100‑1 (MQ=255)
                                                                                             tGCCAGCGTTGGAAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac      <  1:2680846/100‑1 (MQ=255)
                                                                                              gCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac      <  2:2054219/99‑1 (MQ=255)
                                                                                               ccAGCGATGGCAGGTTTTTTTTTTTTTTTTGTTTGCGCATACGAATATGAGTATCTTCCTTATAATTTTGTGGTTACGTCAATAGAGTGCTTTTTcaca     >  2:1137803/1‑99 (MQ=255)
                                                                                               ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa    <  2:667004/100‑1 (MQ=255)
                                                                                               ccAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa    >  1:321147/1‑100 (MQ=255)
                                                                                                cAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAGGTCACTAGAGGGCTTTGTTACAAt   >  2:1442006/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATT  >  NC_000913/4060163‑4060359

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: