Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,060,280 (T)9→8 100% intergenic (+50/‑167) 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,2720T.92.5% 512.4 / 34.3 133intergenic (+42/‑175)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base . (68/55);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 7.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 179 total aligned reads displayed.

ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCA  >  NC_000913/4060173‑4060366
                                                                                                    |                                                                                              
aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTTAT‑TTATTACCTGCAATCGCTGGCAGGt                                                                                                <  2:674957/100‑1 (MQ=255)
aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGt                                                                                                <  2:1953897/100‑1 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  1:441895/1‑100 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  1:1162137/1‑100 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  2:341018/1‑100 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtt                                                                                               >  2:1249188/1‑100 (MQ=255)
  gTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGtt                                                                                               <  1:1852899/99‑1 (MQ=255)
    cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGttttt                                                                                            <  1:1994523/100‑1 (MQ=255)
     aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtttttt                                                                                           <  2:2818476/100‑1 (MQ=255)
     aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtttt                                                                                             >  2:56559/1‑98 (MQ=255)
      ccTGACGGAAAACGTTCGTTGTCGCGCTTACCGTGCAACGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                           <  1:171345/99‑1 (MQ=25)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAATACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                          <  2:1549760/100‑1 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                          >  2:2904830/1‑100 (MQ=255)
      ccTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGttttttt                                                                                          <  2:2207932/100‑1 (MQ=255)
       cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGTACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGtttttta                                                                                          <  1:344061/99‑2 (MQ=38)
       cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                         <  1:3695125/100‑1 (MQ=255)
       cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                         >  1:2827237/1‑100 (MQ=255)
        tGACGGAAAACTTTCTTCGCTGCGCTAATCTCGCACCGAAAGACGATTAATTTCGTTTTTTATTCAT‑AAAAAACTTGTTATCGCTGGCAGGtttttta                                                                                          <  2:271046/98‑2 (MQ=255)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCTCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                        <  1:180461/100‑4 (MQ=38)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                        <  1:208462/100‑4 (MQ=255)
        tGACGGAAAACGATCGTCGCCGCGCCAACCACGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                        <  2:2646438/100‑2 (MQ=255)
         gACTGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                        <  2:655983/100‑4 (MQ=37)
         gACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGttttttttgc                                                                                       >  1:2212335/1‑98 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGttttttttgac                                                                                      <  2:546263/100‑4 (MQ=255)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCATCGCTGGCAGGttttttatgact                                                                                     <  2:3768296/100‑7 (MQ=37)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGttttttttgcat                                                                                     >  1:1638558/1‑96 (MQ=255)
            ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTg                                                                                   >  1:3379183/1‑100 (MQ=255)
            ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTg                                                                                   >  1:902431/1‑100 (MQ=255)
              aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACTATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTc                                                                                 <  1:2402528/100‑1 (MQ=255)
              aaaaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTACTCAT‑AAAAAACCTGCCAGCCACGGCAGTctttttttcc                                                                                       >  2:3757962/1‑93 (MQ=255)
                aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGt                                                                               <  1:3528242/100‑1 (MQ=255)
                aaCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTACCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGt                                                                               >  1:431629/1‑100 (MQ=255)
                   gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAg                                                                            <  2:1699289/100‑1 (MQ=255)
                    cgatctgccgcGCCAACCGCGCACCGAAAGTCGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgc                                                                           <  1:2995110/94‑1 (MQ=255)
                     tcgtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCGGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAg                                                                            <  1:2753996/99‑1 (MQ=255)
                      cgtcgCCGCGCCAACCGCGCACTGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgcgc                                                                         <  1:1485521/100‑1 (MQ=255)
                       gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTAAAAGCAATGGCAGG‑AATTTTTTGCCTGTCGTGAGCGca                                                                        >  1:2196104/1‑100 (MQ=255)
                        tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcac                                                                       <  1:2257241/100‑1 (MQ=255)
                        tcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGGTTTTTAGTCAT‑AAAAGACCTGCCAGCGATCGCACG‑TTTTTTTTGCCTGTCGTGAGCGcac                                                                       <  2:3461371/100‑1 (MQ=255)
                         cgccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGcaca                                                                      <  1:1103681/100‑1 (MQ=255)
                          gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAgcagatcg                                                                     >  2:2995083/1‑94 (MQ=255)
                          gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCacac                                                                     >  2:829087/1‑100 (MQ=255)
                             gcgcCAACCGCGCAACGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                                  >  2:927589/1‑100 (MQ=255)
                                ccAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCATGTAGTGAGCGCACACCAATGc                                                               >  1:3114709/1‑100 (MQ=255)
                                 cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                              >  2:58965/1‑100 (MQ=255)
                                 cAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCg                                                              >  1:2096032/1‑100 (MQ=255)
                                  aaCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGa                                                             >  2:1824498/1‑100 (MQ=255)
                                  aaCCGCGCACCGAAAGACGATTAAAATCATTTTTTAGTGAT‑AAAAAACATGACAGCGATGACAGG‑TTATTTTAACTAGTAGTGAACGCACACCAATgaga                                                             >  2:1016195/1‑100 (MQ=255)
                                   aCCGCGTACCGAAAGACGATTAATTTCGTTGTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCA‑GTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                            <  2:755609/100‑1 (MQ=255)
                                   aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                            <  1:3450391/100‑1 (MQ=255)
                                   aCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAg                                                            >  2:42099/1‑100 (MQ=255)
                                     cgtgtaCCTAAAGACGATTAATTTCGTTTTTTGGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGGCTGTCGTTAGCGCACACCAATGCGAGTa                                                          <  2:1641151/95‑1 (MQ=255)
                                       cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                        <  2:2028937/100‑1 (MQ=255)
                                       cgcACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtc                                                        <  2:73504/100‑1 (MQ=255)
                                         cACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                      >  2:3100297/1‑100 (MQ=255)
                                         cACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                      >  2:2445727/1‑100 (MQ=255)
                                          aCCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                     >  1:599722/1‑100 (MQ=255)
                                           ccGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAATACCTGCCAGCGATGGCAGG‑TTTTTTGTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTccc                                                    <  1:3609338/100‑1 (MQ=255)
                                            cGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCt                                                   <  2:3205365/100‑1 (MQ=255)
                                             ggaaGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGTTTTTTTTT‑CCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                  >  1:101655/3‑100 (MQ=255)
                                             gAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                  >  2:3807252/1‑100 (MQ=255)
                                                      ttAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                         <  1:1008336/100‑1 (MQ=255)
                                                        aaTTTCGTTTTTTAGTCAT‑AAAAAACCTGTCAGCGATGGCAGGTTTTTTTT‑‑CCTGTCGTGAGCGCACACTAATGCGAGTATTTCCCTTATAAGTCTGTgg                                      <  2:1673726/100‑1 (MQ=255)
                                                         aTTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                      >  1:882177/1‑100 (MQ=255)
                                                         aTTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                      <  1:3445799/100‑1 (MQ=255)
                                                         aTTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTgg                                      <  1:2598641/100‑1 (MQ=255)
                                                          tttGGTTTTTTAGTCAT‑AAAACACCTGGCAGCGAGGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACCCCAATGCGAGTATCTCCCTTATAAGTCTGTGGt                                     <  2:3706005/100‑1 (MQ=255)
                                                             cGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTTAGCGCACACCAATGCGAGTATCTCCCTCATAAGTCTGTGATTAc                                  >  2:3000141/1‑100 (MQ=255)
                                                              gTTTTTTAGTTAT‑AAAAATCCTGCCAGTGATGGTAGGTTTTTTTTT‑CCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTAc                                  <  1:766726/99‑1 (MQ=255)
                                                              gTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                 >  1:2800691/1‑100 (MQ=255)
                                                              gTTTTTTAGTCAT‑AAAAAACCTGCCAGATATGTCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                 <  1:2287015/100‑1 (MQ=255)
                                                                      gTCAT‑AAAAAACCTGCCAGCGAGGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAGCTCCCTTATAAGTCTGTGGTTACGTCAATag                          >  2:3898826/1‑99 (MQ=255)
                                                                        gcaaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                        >  2:2774968/5‑100 (MQ=37)
                                                                        cAT‑AAAAAACCTGTCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAATTCTGTGGTTAAGTCAATata                         >  2:549861/1‑96 (MQ=255)
                                                                        cAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       <  2:2479935/100‑1 (MQ=255)
                                                                        cAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGACAATagaaa                       >  1:3680476/1‑98 (MQ=255)
                                                                        cAT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAATATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       >  1:596486/1‑100 (MQ=255)
                                                                         caaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCACTTATAAGTCTGTGGTTACGTCAATagaga                       >  1:1532326/4‑100 (MQ=37)
                                                                         aT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      <  1:804374/100‑1 (MQ=255)
                                                                         aT‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  1:1888897/1‑100 (MQ=255)
                                                                          t‑AAAAAACCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       >  1:2837469/1‑98 (MQ=255)
                                                                          aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  1:3344987/3‑100 (MQ=38)
                                                                          aaaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  1:2957881/3‑100 (MQ=38)
                                                                           aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                     >  1:2725461/2‑100 (MQ=38)
                                                                           aaaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTAACGTCAATAGAGAGc                     >  1:3875849/2‑100 (MQ=38)
                                                                            aaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTT‑CCTGTCATGAGCGCACACCCATGCAAGTATCTCTGTTATCAGTCTGTGGTTACGTAAATAGACAGCt                    >  2:3237417/1‑100 (MQ=21)
                                                                            aaaaaaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                    >  1:2324585/1‑100 (MQ=38)
                                                                            aaaaaaCCTGCCAGCGATGGCAGGTTTTTTTTT‑TTTTTCTTGAGCGCACAACAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                    >  2:191447/1‑100 (MQ=255)
                                                                            aaaaaaCCTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGCGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCt                    >  2:191522/1‑100 (MQ=255)
                                                                              aaaaCCTGCCATCGCTGGCAGG‑TTTTTTATGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCttt                  >  1:1625109/1‑100 (MQ=38)
                                                                                aaCCTGCCATCGCTGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAATCTGTGGTTACGTCAATAGAGAGCTTTGt                >  1:532916/1‑100 (MQ=38)
                                                                                  ccTGCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca              >  1:1744835/1‑100 (MQ=255)
                                                                                   cTGCCAGCGATGGCAGG‑TTTTTTTTGCTTGTCGTGAGCGCACACCAATGCGAGTATCTCTCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac             >  2:2556234/1‑100 (MQ=255)
                                                                                    tGCCAGTGATGGCAGG‑TTTTTTTTGCCTGTCTTGAGCGCATACCAATGCGAGTATCTCCCTAATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcaca            <  1:583367/100‑1 (MQ=255)
                                                                                     gCCAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcataa           >  2:2163903/1‑97 (MQ=255)
                                                                                      ccAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt          >  2:3286128/1‑100 (MQ=255)
                                                                                       cAGCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAtt         >  2:1757853/1‑100 (MQ=255)
                                                                                         gCGATGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTAt       >  1:3748646/1‑100 (MQ=255)
                                                                                             tGGCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGc   <  2:2202317/100‑1 (MQ=255)
                                                                                               gCAGG‑TTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa  <  1:659932/99‑1 (MQ=255)
                                                                                                    |                                                                                              
ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCAT‑AAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCA  >  NC_000913/4060173‑4060366

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