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.

GTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATT  >  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)
          aCGTCACCTGACGGAAATCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTAAAAAAAAAAAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                         <  2:2760008/100‑1 (MQ=25)
          aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCTATCGCTGGCAgtt                                                                                          <  1:2569307/100‑3 (MQ=255)
          aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑t                                                                                         <  1:2822258/100‑1 (MQ=255)
          aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑t                                                                                         <  2:1584690/100‑1 (MQ=255)
           cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tt                                                                                        <  1:345862/100‑1 (MQ=255)
           cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑tt                                                                                        >  1:2128866/1‑100 (MQ=255)
            gTCACCTGACGGAAAACGATCGTCGCCGCGCTAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑ttt                                                                                       <  2:486112/100‑1 (MQ=255)
            gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttt                                                                                       <  1:197877/100‑1 (MQ=255)
             tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttt                                                                                      <  2:2773564/100‑1 (MQ=255)
                ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑tttttta                                                                                   <  1:3437051/100‑2 (MQ=255)
                    tcGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGc                                                                               <  1:1409532/99‑1 (MQ=255)
                    aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatgac                                                                               <  1:3659308/100‑6 (MQ=38)
                    aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGG‑ttttttatgac                                                                               <  1:1548713/100‑6 (MQ=38)
                     cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                              >  2:2244178/1‑100 (MQ=255)
                     cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGcc                                                                              >  2:737418/1‑100 (MQ=255)
                     cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGca                                                                              >  2:1409522/1‑99 (MQ=255)
                      ggAAAACGATTGTCGCCGCGCCAACCGCGCAGCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTg                                                                                <  2:1682913/97‑1 (MQ=255)
                        aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGt                                                                           >  1:38925/1‑100 (MQ=255)
                         aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                          >  1:2453127/1‑100 (MQ=255)
                         aaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTc                                                                          >  1:3134739/1‑100 (MQ=255)
                            cGATCGTCGCCGCGCCACCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCg                                                                         <  2:437023/98‑1 (MQ=255)
                            cGATCGTCGCCACGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCCGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTg                                                                       >  2:2002298/1‑100 (MQ=255)
                             gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAACTTCGTCTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑CTTTTTTTTGCCTGTCCTGa                                                                      >  2:868946/1‑100 (MQ=255)
                               tcgtcgCCGCGCCAACCGTGCACCGAAAGTCGATTAATTTCTTTTGTTAGTCATAAAAAGCCTGCCAGCGATGGCAGG‑TTTTGCTTTGCCTGTCGTGAg                                                                     <  2:711944/99‑1 (MQ=255)
                                 gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAgcgc                                                                  <  1:607821/100‑1 (MQ=255)
                                    gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGcac                                                                >  1:1757371/1‑100 (MQ=255)
                                    gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGcaca                                                               >  1:2797679/1‑100 (MQ=255)
                                      cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                             <  1:917241/100‑1 (MQ=255)
                                      cgcgCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                             <  2:2943463/100‑1 (MQ=255)
                                          ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCTTGTGTTGGGCGCccacca                                                            >  1:3725681/1‑98 (MQ=255)
                                           cAACCGCGCACCGAAGGCCGATTAATGTCGTTTTTTAGTCATAAAATACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                             <  2:3125217/95‑1 (MQ=255)
                                           cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGc                                                        >  1:1908412/1‑100 (MQ=255)
                                              ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAGCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     <  2:1863153/100‑1 (MQ=255)
                                              ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     <  2:557929/100‑1 (MQ=255)
                                              ccGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     >  1:2631400/1‑100 (MQ=255)
                                               cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                     <  2:2902084/99‑1 (MQ=255)
                                                gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                   <  2:2186783/100‑1 (MQ=255)
                                                gcgcACCGAAAGACAATTAATTTCGTTTTTTTGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCCGTCGTGAGCGCACACCAATGCGAGTa                                                   <  2:352381/100‑1 (MQ=255)
                                                 cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                   <  1:1066086/100‑1 (MQ=255)
                                                   cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                >  1:362014/1‑100 (MQ=255)
                                                      cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                             <  1:661379/100‑1 (MQ=255)
                                                      cGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                             >  1:3233913/1‑100 (MQ=255)
                                                       gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                             >  1:3489786/1‑99 (MQ=255)
                                                         atgACGTTTATTTCCGTTTTTTAGTCATAAAAAACTTCCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                          <  2:2320182/98‑1 (MQ=255)
                                                         aaGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                          <  1:2076650/100‑1 (MQ=255)
                                                          aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                         <  2:3489759/100‑1 (MQ=255)
                                                          aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtat                                         >  2:2856310/1‑100 (MQ=255)
                                                           gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                        <  2:1536403/100‑1 (MQ=255)
                                                            aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                       <  1:73847/100‑1 (MQ=255)
                                                             cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                      <  2:589923/100‑1 (MQ=255)
                                                             cGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                      <  2:3241198/100‑1 (MQ=255)
                                                              gATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGAACACCAATGCGAGTATCTCCCTTATAAg                                      >  2:41741/1‑99 (MQ=255)
                                                               aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGt                                     <  2:2355118/100‑1 (MQ=255)
                                                               aTTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTc                                    >  1:1250399/1‑100 (MQ=255)
                                                                 tAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCt                                   <  2:1951423/100‑1 (MQ=255)
                                                                  aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgt                                 <  2:1191702/100‑1 (MQ=255)
                                                                    tttCGTTTTTTAGTCATAAAAGACCTGCCGGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGGGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                               <  2:1351216/100‑1 (MQ=255)
                                                                    tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                               <  2:1328613/100‑1 (MQ=255)
                                                                     ttCGTTTTTTAGTCATAAAACACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                              <  2:1652460/100‑1 (MQ=255)
                                                                     ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTACCTGTCGTGAGCGCACACCAATGCGTGTATCTCCCTTATATGTCTGTGGt                              >  2:767689/1‑100 (MQ=255)
                                                                        gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                           >  2:1755803/1‑100 (MQ=255)
                                                                         ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                          <  2:2629996/100‑1 (MQ=255)
                                                                               aGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATa                    <  2:1599363/100‑1 (MQ=255)
                                                                                gcaaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                   >  1:1829791/6‑100 (MQ=38)
                                                                                gcaaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                   >  1:2389119/6‑100 (MQ=38)
                                                                                 caaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                  >  1:2414999/5‑100 (MQ=255)
                                                                                  cATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                 >  1:597492/1‑100 (MQ=255)
                                                                                  aaaaaaaaaCCTGCCATCACTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                  >  1:2326262/4‑99 (MQ=38)
                                                                                   aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                >  1:3073154/3‑100 (MQ=255)
                                                                                   aTAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                <  2:1523476/100‑1 (MQ=255)
                                                                                    aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag               >  1:1210445/2‑100 (MQ=255)
                                                                                     aaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc              >  1:2207004/1‑100 (MQ=255)
                                                                                     aaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc              >  1:2252024/1‑100 (MQ=255)
                                                                                       aaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTTTTTTTTCTTTATCGCACACCAATGCGCGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt            >  2:3447288/1‑100 (MQ=25)
                                                                                       aaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt            >  1:3424982/1‑100 (MQ=255)
                                                                                        aaaCCTGCCAGCGATGGCAGG‑TTTTTTTTTTTTTGTTTTTATCGCACACCAATGCGCGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt           >  2:2546819/1‑100 (MQ=255)
                                                                                         aaCCTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTg          >  2:3815271/1‑100 (MQ=255)
                                                                                         aaCCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt         <  1:103045/100‑1 (MQ=255)
                                                                                          aCCTGCCAGCGATGGCAGG‑TTTTTTTTTTCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt         >  1:1662847/1‑100 (MQ=255)
                                                                                          aCCTGCCAGCGATGGCAGG‑‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGt         <  1:878991/99‑1 (MQ=255)
                                                                                           ccTGCCATCGCTGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc        >  1:3195262/1‑100 (MQ=255)
                                                                                           ccTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc        <  2:1217528/100‑1 (MQ=255)
                                                                                            cTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca       <  2:3603964/100‑1 (MQ=255)
                                                                                            cTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca       <  1:3468588/100‑1 (MQ=255)
                                                                                             tGCCAGCGTTGGAAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac      <  1:2680846/100‑1 (MQ=255)
                                                                                              gCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac      <  2:2054219/99‑1 (MQ=255)
                                                                                               ccAGCGATGGCAGG‑TTTTTTTTTTTTTTTTGTTTGCGCATACGAATATGAGTATCTTCCTTATAATTTTGTGGTTACGTCAATAGAGTGCTTTTTcaca     >  2:1137803/1‑99 (MQ=255)
                                                                                               ccAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa    <  2:667004/100‑1 (MQ=255)
                                                                                               ccAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa    >  1:321147/1‑100 (MQ=255)
                                                                                                cAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAGGTCACTAGAGGGCTTTGTTACAAt   >  2:1442006/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
GTATTCGTAAACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATT  >  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: ‑