Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,954,256 T→A 6.6% I272N (ATT→AAT)  cmoB → tRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,954,2560TA 6.6% 379.9 / 11.6 105I272N (ATT→AAT) cmoBtRNA (cmo5U34)‑carboxymethyltransferase, carboxy‑SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (50/48);  new base A (3/4);  total (53/52)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 105 total aligned reads displayed.

GAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATT  >  NC_000913/1954114‑1954398
                                                                                                                                              |                                                                                                                                              
gAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCgg                                                                                                                                           >  4:848316/1‑148 (MQ=255)
 aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                          >  4:78041/1‑148 (MQ=255)
 aaaCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                          <  3:78041/148‑1 (MQ=255)
   aCGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAta                                                                                                                                        <  4:145626/148‑2 (MQ=255)
     gCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgtg                                                                                                                                      <  3:848316/148‑1 (MQ=255)
         gTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAt                                                                                                                                         >  2:139909/1‑141 (MQ=255)
         gTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCGGAt                                                                                                                                         <  1:139909/141‑1 (MQ=255)
              tGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacc                                                                                                                             >  3:450855/1‑148 (MQ=255)
                aTGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                          <  3:175984/133‑1 (MQ=255)
                aTGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                          >  4:175984/1‑133 (MQ=255)
                  ggTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgtg                                                                                                                                      >  1:1106027/1‑135 (MQ=255)
                  ggTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGAtgtg                                                                                                                                      <  2:1106027/135‑1 (MQ=255)
                        gagaACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAACCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaag                                                                                                                   <  3:464761/148‑1 (MQ=255)
                            acacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGAAATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagca                                                                                                               <  4:450855/148‑1 (MQ=255)
                             cacGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGagcag                                                                                                              >  3:284431/1‑148 (MQ=255)
                                  tGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                          >  1:414531/1‑115 (MQ=255)
                                  tGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGa                                                                                                                                          <  2:414531/115‑1 (MQ=255)
                                             ggCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGGGCGTTACCACCACAGAAGAGCAGCGa                                                                                                           >  3:203675/1‑135 (MQ=255)
                                             ggCGATCGTTACGCGCAAATGCGTAATGTCTATTTCAGCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGGGCGTTACCACCACAGAAGAGCAGCGa                                                                                                           <  4:203675/135‑1 (MQ=255)
                                                   cGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTaccac                                                                                                                           >  4:775548/1‑113 (MQ=255)
                                                   cGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTaccac                                                                                                                           <  3:775548/113‑1 (MQ=255)
                                                   cGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgg                                                                                               <  3:955180/141‑1 (MQ=255)
                                                             aaaTGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCt                                                                              >  2:23015/1‑148 (MQ=255)
                                                                tGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCg                                                                                                            <  4:851339/115‑1 (MQ=255)
                                                                tGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCg                                                                                                            >  3:851339/1‑115 (MQ=255)
                                                                tGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGc                                                                           >  1:54843/1‑148 (MQ=255)
                                                                    tAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGAt                                                                       <  2:54843/148‑1 (MQ=255)
                                                                     aaTGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                  >  2:61637/1‑120 (MQ=255)
                                                                     aaTGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                  <  1:285190/120‑1 (MQ=255)
                                                                     aaTGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                  >  2:285190/1‑120 (MQ=255)
                                                                       tGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTa                                                                                                                               >  1:1114720/1‑89 (MQ=255)
                                                                       tGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTa                                                                                                                               <  2:1114720/89‑1 (MQ=255)
                                                                       tGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACaga                                                                                                                     >  3:645301/1‑99 (MQ=255)
                                                                       tGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACaga                                                                                                                     <  4:645301/99‑1 (MQ=255)
                                                                            aTTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGATCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGa                                                               >  2:345166/1‑148 (MQ=255)
                                                                            aTTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCg                                                                                                                                            <  4:1112134/71‑1 (MQ=255)
                                                                            aTTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCAATGCg                                                                                                                                            >  3:1112134/1‑71 (MQ=255)
                                                                              ttCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacc                                                                                                                             >  1:1006835/1‑84 (MQ=255)
                                                                              ttCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTacc                                                                                                                             <  2:1006835/84‑1 (MQ=255)
                                                                                   cccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTaccacc                                                                                                                          >  1:863523/1‑82 (MQ=255)
                                                                                   cccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTaccacc                                                                                                                          <  2:863523/82‑1 (MQ=255)
                                                                                   cccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGAttt                                                                     >  1:691048/1‑135 (MQ=255)
                                                                                   cccTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGAttt                                                                     <  2:691048/135‑1 (MQ=255)
                                                                                     cTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                  >  4:534711/1‑104 (MQ=255)
                                                                                     cTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                  >  2:64504/1‑104 (MQ=255)
                                                                                     cTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                  <  3:534711/104‑1 (MQ=255)
                                                                                     cTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGaa                                                                                                  <  1:64504/104‑1 (MQ=255)
                                                                                     cTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTc                                                                    >  3:52475/1‑134 (MQ=255)
                                                                                     cTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTc                                                                    <  4:52475/134‑1 (MQ=255)
                                                                                       tCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTg                                                                             >  1:1114714/1‑123 (MQ=255)
                                                                                       tCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTg                                                                             <  2:1114714/123‑1 (MQ=255)
                                                                                       tCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGc                                                                           >  4:68814/1‑125 (MQ=255)
                                                                                       tCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGc                                                                           <  3:68814/125‑1 (MQ=255)
                                                                                                   ctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGATCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACgg                                        <  1:345166/148‑1 (MQ=255)
                                                                                                   ctgAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACgg                                        <  1:23015/148‑1 (MQ=255)
                                                                                                     gAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaag                                                                                                                   >  1:1094553/1‑71 (MQ=255)
                                                                                                     gAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACagaag                                                                                                                   <  2:1094553/71‑1 (MQ=255)
                                                                                                        aaaCTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                   <  4:284431/148‑1 (MQ=255)
                                                                                                            tGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                               >  3:903593/1‑148 (MQ=255)
                                                                                                                 gaagaaGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTg                                                                             >  1:655826/1‑97 (MQ=255)
                                                                                                                 gaagaaGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTg                                                                             <  2:655826/97‑1 (MQ=255)
                                                                                                                   agaagTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                   >  1:416400/1‑137 (MQ=255)
                                                                                                                           ggTTTTTTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgga                                                                                              >  2:615930/1‑70 (MQ=255)
                                                                                                                           ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgga                                                                                              <  1:615930/70‑1 (MQ=255)
                                                                                                                           ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAAt                                                                                  >  4:466447/1‑82 (MQ=255)
                                                                                                                           ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                               <  4:80971/85‑1 (MQ=255)
                                                                                                                           ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCca                                                           >  2:854596/1‑105 (MQ=255)
                                                                                                                           ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCca                                                           <  1:854596/105‑1 (MQ=255)
                                                                                                                           ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                       >  3:537631/1‑125 (MQ=255)
                                                                                                                           ggTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGt                                       <  4:537631/125‑1 (MQ=255)
                                                                                                                           ggTTTTGTAGAAAACCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                               >  3:80971/1‑85 (MQ=255)
                                                                                                                              tttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                               >  2:170391/1‑82 (MQ=255)
                                                                                                                              tttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                               <  1:170391/82‑1 (MQ=255)
                                                                                                                              tttGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcgcg             >  2:825374/1‑148 (MQ=255)
                                                                                                                                  tAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                               <  1:145287/126‑1 (MQ=255)
                                                                                                                                  tAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                               >  2:145287/1‑126 (MQ=255)
                                                                                                                                   aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTc                                                                    <  4:953373/88‑1 (MQ=255)
                                                                                                                                   aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTc                                                                    >  3:953373/1‑88 (MQ=255)
                                                                                                                                   aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACgg                                        <  2:58834/116‑1 (MQ=255)
                                                                                                                                   aGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACgg                                        >  1:58834/1‑116 (MQ=255)
                                                                                                                                    cacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                               <  3:773273/73‑1 (MQ=255)
                                                                                                                                    cacatCCGCAATGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGc                                                                               >  4:773273/4‑76 (MQ=255)
                                                                                                                                      tatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                               >  2:489383/1‑122 (MQ=255)
                                                                                                                                      tatCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGtt                               <  1:489383/122‑1 (MQ=255)
                                                                                                                                       atCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAAcg               >  3:178629/1‑137 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                           >  1:216573/1‑107 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTaaaa                                           <  2:216573/107‑1 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTAt                             <  3:771820/121‑1 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTAt                             >  4:771820/1‑121 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGAtt  >  2:72857/1‑148 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCAACGCACCGAATGGATGGTTACCGAATCGCTgg                                                                            >  1:340818/1‑74 (MQ=255)
                                                                                                                                         ccGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCAACGCACCGAATGGATGGTTACCGAATCGCTgg                                                                            <  2:340818/74‑1 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                   <  4:836925/113‑1 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGaa                                   >  3:836925/1‑113 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGc                        <  1:108096/124‑1 (MQ=255)
                                                                                                                                           gCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGc                        >  2:108096/1‑124 (MQ=255)
                                                                                                                                            cATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgg                                                                                               >  4:878993/1‑52 (MQ=255)
                                                                                                                                            cATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAatgg                                                                                               <  3:878993/52‑1 (MQ=255)
                                                                                                                                            cATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACg                                         <  3:210112/106‑1 (MQ=255)
                                                                                                                                            cATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACg                                         >  4:210112/1‑106 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GAAACGCTGGTTATTGATGGTGATGAGAACACGGTGCTGGTACCAGGCGATCGTTACGCGCAAATGCGTAATGTCTATTTCATCCCTTCCGCGCTGGCGCTGAAAAACTGGCTGAAGAAGTGTGGTTTTGTAGATATCCGCATTGCGGATGTGAGCGTTACCACCACAGAAGAGCAGCGACGCACCGAATGGATGGTTACCGAATCGCTGGCCGATTTCCTCGACCCACATGATCCGGGTAAAACGGTGGAAGGTTATCCGGCACCTAAACGCGCGGTGCTGATT  >  NC_000913/1954114‑1954398

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: