Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,060,280 (T)9→8 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.95.1% 412.7 / 14.0 102intergenic (+42/‑175)typA/yihLGTP‑binding protein/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base . (57/40);  total (60/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 134 total aligned reads displayed.

ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCA  >  NC_000913/4060173‑4060366
                                                                                                    |                                                                                              
aCGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGt                                                                                                <  1:2351395/100‑1 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtt                                                                                               <  1:1006009/100‑1 (MQ=255)
 cGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtt                                                                                               <  2:2986443/100‑1 (MQ=255)
   tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                             <  1:1497828/100‑1 (MQ=255)
   tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                             <  1:2898559/100‑1 (MQ=255)
   tCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttt                                                                                             <  1:1033429/100‑1 (MQ=255)
    caacTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAATCCTGCCATCGCTGGCAGGtttt                                                                                             <  2:1791282/96‑1 (MQ=39)
    cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                            >  2:1477670/1‑100 (MQ=255)
    cACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATGAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttt                                                                                            <  2:315987/100‑1 (MQ=255)
     aCCTGACGGAAAACGATCGTCGCCGCGTTATCCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAGACCTGCCATCGCTGGCAGGtttttt                                                                                           <  2:2449584/100‑1 (MQ=38)
     aCCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGtttttt                                                                                           <  2:2552036/100‑1 (MQ=255)
       cTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGtttttttt                                                                                         <  1:620805/100‑1 (MQ=255)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCATCGCTGGCAGGttttttatg                                                                                        <  2:1468501/100‑4 (MQ=38)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                        >  1:3010622/1‑99 (MQ=255)
        tGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                        <  1:2459968/100‑2 (MQ=255)
         gACGGAAAAGGATCGTCGCCGCGCCTACCGCGGACCGAAAGCCGATTAATTTCGTTTTTTAGTCATAGAAAACCTGCCAGCGATGGCAGGttttttttg                                                                                        <  1:996625/99‑2 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAACCCTGCCAGCGATGGCAGGTTTTTTTTTGc                                                                                      <  1:240515/100‑1 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                      <  2:2878831/100‑4 (MQ=255)
          aCGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcc                                                                                      <  1:2225741/100‑4 (MQ=255)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                                     >  2:2580466/1‑96 (MQ=255)
           cGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGttttttttgcct                                                                                     >  1:875251/1‑96 (MQ=255)
            ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                                   >  2:1748902/1‑100 (MQ=255)
            ggAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTg                                                                                   >  1:528895/1‑100 (MQ=255)
             gAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGAAGGttttttttgcct                                                                                     >  1:2884176/1‑94 (MQ=255)
                 aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTg                                                                              <  2:191167/100‑1 (MQ=255)
                 aCGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTg                                                                              <  2:480308/100‑1 (MQ=255)
                   gATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAg                                                                            >  1:392194/1‑100 (MQ=255)
                    aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                           >  1:1891760/1‑100 (MQ=255)
                    aTCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGGCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAgc                                                                           <  1:131704/100‑1 (MQ=255)
                       gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGca                                                                        <  2:1707877/100‑1 (MQ=255)
                       gtcgCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGca                                                                        >  1:45053/1‑100 (MQ=255)
                          gccgcGCCAACCGCGCACCGAAAGACGATTCATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTAGTGAGCGCacac                                                                     >  2:1335135/1‑100 (MQ=255)
                          gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                                     >  2:1330616/1‑100 (MQ=255)
                          gccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCacac                                                                     <  1:960796/100‑1 (MQ=255)
                           ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                    <  2:1469031/100‑1 (MQ=255)
                           ccgcGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACAcc                                                                    >  1:2095429/1‑100 (MQ=255)
                             gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                                  >  1:1972700/1‑100 (MQ=255)
                             gcgcCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCaa                                                                  >  1:735902/1‑100 (MQ=255)
                                     cgcgCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTa                                                          <  1:2746547/100‑1 (MQ=255)
                                      gcgcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAt                                                         >  1:1139297/1‑100 (MQ=255)
                                        gcACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtct                                                       >  2:1175631/1‑100 (MQ=255)
                                         cACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTAtctc                                                      >  1:116135/1‑100 (MQ=255)
                                          aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                     <  1:595216/100‑1 (MQ=255)
                                          aCCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTcc                                                     >  2:640668/1‑100 (MQ=255)
                                             gAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCtt                                                  >  1:1756867/1‑100 (MQ=255)
                                              aaaGACGATTAATTTCGTTTTTTAGTTATAAAATACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTta                                                 >  2:2036013/1‑100 (MQ=255)
                                                aGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTtata                                               <  2:457631/100‑1 (MQ=255)
                                                 gACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATaa                                              >  2:638461/1‑100 (MQ=255)
                                                  aCGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAg                                             <  1:1785428/100‑1 (MQ=255)
                                                      ttAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtg                                         <  2:68463/100‑1 (MQ=255)
                                                        aaTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                       <  1:1518052/100‑1 (MQ=255)
                                                         aTTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGAGAGTATCTCCCTTATAAGTCTGTgg                                      >  2:1357360/1‑100 (MQ=255)
                                                          tttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                       <  1:248533/100‑1 (MQ=255)
                                                          tttCGTTTTTTAATCATAAAAAACCTCCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCtgtg                                       <  2:1257039/98‑1 (MQ=255)
                                                           ttCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGtt                                    >  2:2282702/1‑100 (MQ=255)
                                                            tCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTa                                   >  1:2713619/1‑100 (MQ=255)
                                                              gTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                 <  1:2580489/100‑1 (MQ=255)
                                                               ttttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACg                                 <  2:2095405/99‑1 (MQ=255)
                                                                tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                               <  1:2993495/100‑1 (MQ=255)
                                                                tttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTc                               <  2:1939352/100‑1 (MQ=255)
                                                                 ttttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCa                              <  1:2725689/100‑1 (MQ=255)
                                                                   ttAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                            >  2:2362148/1‑100 (MQ=255)
                                                                   acaggcaaaaaaaaCCTGCCATTGCTGGCAGGTTTTTTATGACTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAAt                            >  2:2942979/9‑100 (MQ=21)
                                                                     aGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATag                          <  2:735898/100‑1 (MQ=255)
                                                                      ggcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATaga                         >  2:1281704/6‑100 (MQ=37)
                                                                       tCATAAAAAACCTGCCAGCGATCGCAGGTTTCTTTTGTCTGTCGTGAACGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                        <  2:2084116/100‑1 (MQ=255)
                                                                       gcaaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagag                        >  1:2908060/5‑100 (MQ=37)
                                                                        naTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       >  2:1640785/2‑100 (MQ=255)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       >  1:1370948/1‑100 (MQ=255)
                                                                        cATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagaga                       <  1:1073803/100‑1 (MQ=255)
                                                                         aaaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  2:2913138/3‑100 (MQ=38)
                                                                         aTAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATagagag                      >  1:2223231/1‑100 (MQ=255)
                                                                          tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                     >  1:3055035/1‑100 (MQ=255)
                                                                          tAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCACTTATAAGTATGTGGTTACGTCCATAGAGAGc                     >  2:1671397/1‑100 (MQ=255)
                                                                          aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                     <  2:2667292/99‑1 (MQ=38)
                                                                          aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                     >  1:2006426/2‑100 (MQ=38)
                                                                          aaaaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGc                     >  2:2138226/2‑100 (MQ=38)
                                                                            aaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                   >  1:1482859/1‑100 (MQ=38)
                                                                            aaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                   >  1:914921/1‑100 (MQ=38)
                                                                            aaaaaCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                   >  1:1436319/1‑100 (MQ=38)
                                                                            aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                   <  2:1285815/100‑1 (MQ=255)
                                                                            aaaaaCCTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCtt                   <  2:2764338/100‑1 (MQ=255)
                                                                             aaaaCCTGCCAGCGATGGCAGGTTTTTTTTGTTTTTTTTTAACGCACACCAATGCGAGTATCTTTCTTATTAGTCTGTGGTTCCGTCATTTGAGAGCttt                  >  2:1199639/1‑100 (MQ=255)
                                                                                aCCTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc               >  1:2641424/1‑100 (MQ=255)
                                                                                 ccTGCCATCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTc               >  1:1216455/1‑99 (MQ=255)
                                                                                 ccTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTca              <  2:1587323/100‑1 (MQ=255)
                                                                                 catgacAGCGCTGGCAGGTTTTTTTTGCCTGTCGTGAGCGAACACCAATGAGAGTATCTCCAGTATAAGTCTGTGTTTACGTCAATAGAGACCTTGGTca              <  2:2884147/95‑1 (MQ=255)
                                                                                  cTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac             <  2:45053/100‑1 (MQ=255)
                                                                                  cTGCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTcac             <  1:2282725/100‑1 (MQ=255)
                                                                                    gCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa           <  2:2425322/100‑1 (MQ=255)
                                                                                    gCCAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACaa           >  1:2551138/1‑100 (MQ=255)
                                                                                     ccAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt          >  1:2990684/1‑100 (MQ=255)
                                                                                     ccAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt          <  1:607359/100‑1 (MQ=255)
                                                                                     ccAGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAAt          >  1:2374581/1‑100 (MQ=255)
                                                                                       aGCGATGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTa        >  1:731419/1‑100 (MQ=255)
                                                                                           aTGGCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTg    <  1:2415507/100‑1 (MQ=255)
                                                                                             ggCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa  >  2:2784566/1‑100 (MQ=255)
                                                                                             ggCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa  >  1:696926/1‑100 (MQ=255)
                                                                                             ggCAGGTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCa  >  1:698869/1‑100 (MQ=255)
                                                                                              gCAGGTTTTTTTTTTTTTTTTTGAGCGCACACCAATGCGAGTATCTTCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCt     >  2:437057/1‑96 (MQ=255)
                                                                                                    |                                                                                              
ACGTCACCTGACGGAAAACGATCGTCGCCGCGCCAACCGCGCACCGAAAGACGATTAATTTCGTTTTTTAGTCATAAAAAACCTGCCAGCGATGGCAGGTTTTTTTTTGCCTGTCGTGAGCGCACACCAATGCGAGTATCTCCCTTATAAGTCTGTGGTTACGTCAATAGAGAGCTTTGTCACAATTATCTGCA  >  NC_000913/4060173‑4060366

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