Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 561,644 A→G Y27C (TAT→TGT)  sfmH → involved in fimbrial asembly

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913561,6440AG92.4% 612.8 / 32.2 197Y27C (TAT→TGT) sfmHinvolved in fimbrial asembly
Reads supporting (aligned to +/- strand):  ref base A (7/8);  new base G (76/106);  total (83/114)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 202 total aligned reads displayed.

ATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAG  >  NC_000913/561552‑561733
                                                                                            |                                                                                         
  tattattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACtgtg                                                                                         >  1:4763526/1‑93 (MQ=255)
   ntnattGNCGATGGCATGTTTNTNTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGa                                                                                        <  2:5331066/92‑1 (MQ=255)
   attattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGa                                                                                        <  2:5966372/93‑1 (MQ=255)
   attattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGa                                                                                        <  2:13393/93‑1 (MQ=255)
   attattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGa                                                                                        >  2:1427585/1‑93 (MQ=255)
    ttattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAt                                                                                       >  2:474864/1‑93 (MQ=255)
    ttattGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAt                                                                                       <  2:1763368/93‑1 (MQ=255)
     tattGACGATGACATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAtc                                                                                      <  1:5839335/93‑1 (MQ=255)
       ttGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAtctc                                                                                    <  2:1961138/93‑1 (MQ=255)
        tGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAtctct                                                                                   <  1:5730646/93‑1 (MQ=255)
         ggCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAtctctc                                                                                  <  1:3802161/93‑1 (MQ=255)
         ggCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAtctctc                                                                                  >  1:1156688/1‑93 (MQ=255)
         ggCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAtctctc                                                                                  <  1:1122220/93‑1 (MQ=255)
         ggCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGAtctctc                                                                                  >  2:78269/1‑93 (MQ=255)
            gATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCaaa                                                                               <  1:254798/93‑1 (MQ=255)
             aTGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAat                                                                              >  2:2270229/1‑93 (MQ=255)
             aTGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAat                                                                              <  2:4398930/93‑1 (MQ=255)
             aTGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAat                                                                              >  2:5810862/1‑93 (MQ=255)
             aTGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAat                                                                              <  1:3950816/93‑1 (MQ=255)
             aTGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAat                                                                              >  2:1068366/1‑93 (MQ=255)
              tGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAata                                                                             >  1:1715198/1‑93 (MQ=255)
              tGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAata                                                                             <  1:2705858/93‑1 (MQ=255)
              tGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAata                                                                             >  2:4096449/1‑93 (MQ=255)
                gCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATc                                                                           >  1:455086/1‑93 (MQ=255)
                gCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATc                                                                           <  2:3666324/93‑1 (MQ=255)
                 cATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCt                                                                          >  2:5005457/1‑93 (MQ=255)
                  aTGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCtt                                                                         <  1:1886930/93‑1 (MQ=255)
                   tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCttt                                                                        >  1:242009/1‑93 (MQ=255)
                   tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCttt                                                                        <  1:4557925/93‑1 (MQ=255)
                   tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCttt                                                                        <  2:4192294/93‑1 (MQ=255)
                   tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCttt                                                                        >  1:543655/1‑93 (MQ=255)
                   tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCttt                                                                        <  2:1551083/93‑1 (MQ=255)
                   tGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCttt                                                                        >  1:1963190/1‑93 (MQ=255)
                     tttGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTaa                                                                      <  1:1380765/93‑1 (MQ=255)
                     tttGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTaa                                                                      >  2:2040443/1‑93 (MQ=255)
                     tttGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTaa                                                                      >  2:6038409/1‑93 (MQ=255)
                       tgtgTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATa                                                                    <  1:918649/93‑1 (MQ=255)
                         tgtCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGc                                                                  <  1:426145/93‑1 (MQ=255)
                         tgtCTGGCAAACATATCCGGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGc                                                                  <  1:3038898/93‑1 (MQ=255)
                               gccaaCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAAt                                                            <  1:291738/90‑1 (MQ=255)
                                gaaaCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATa                                                           <  2:2881502/92‑1 (MQ=255)
                                cAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATa                                                           >  2:2724718/1‑93 (MQ=255)
                                  aaCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAc                                                         <  1:4901856/93‑1 (MQ=255)
                                  aaCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAc                                                         >  1:2129960/1‑93 (MQ=255)
                                   aCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAcc                                                        <  1:4690687/93‑1 (MQ=255)
                                   aCATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAcc                                                        >  1:1621647/1‑93 (MQ=255)
                                    cATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAcca                                                       <  1:2871538/93‑1 (MQ=255)
                                     atatCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAccag                                                      <  1:1481042/93‑1 (MQ=255)
                                       atCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAccagcc                                                    <  2:167793/93‑1 (MQ=255)
                                       atCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAccagcc                                                    <  2:3779566/93‑1 (MQ=255)
                                       atCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagcc                                                    <  1:5809053/93‑1 (MQ=255)
                                        tCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAAccagcca                                                   <  1:5134626/93‑1 (MQ=255)
                                         ccTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAccagccag                                                  <  1:2320448/93‑1 (MQ=255)
                                         ccTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAAccagccag                                                  <  1:3263890/93‑1 (MQ=255)
                                          cTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAgg                                                 <  1:2587889/93‑1 (MQ=255)
                                           tGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAggg                                                <  1:376505/93‑1 (MQ=255)
                                            gggCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGc                                               <  2:4552740/93‑1 (MQ=255)
                                            gggCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGc                                               >  1:3635093/1‑93 (MQ=255)
                                             ggCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCa                                              <  1:2819807/93‑1 (MQ=255)
                                              gCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAg                                             <  2:5358051/93‑1 (MQ=255)
                                                tACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAt                                           <  1:4140994/93‑1 (MQ=255)
                                                tACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAt                                           >  2:2748068/1‑93 (MQ=255)
                                                tACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAt                                           >  2:3298706/1‑93 (MQ=255)
                                                  cTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttg                                         >  2:4255975/1‑93 (MQ=255)
                                                   tGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgt                                        <  2:5440148/93‑1 (MQ=255)
                                                   tGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgt                                        >  1:917596/1‑93 (MQ=255)
                                                    gTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgtt                                       <  1:1281420/93‑1 (MQ=255)
                                                    gTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGAttgtt                                       >  2:4871554/1‑93 (MQ=255)
                                                      ttGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATtgttgt                                     <  2:1636842/93‑1 (MQ=255)
                                                       tgtgCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTgttgtt                                    <  1:4834948/93‑1 (MQ=255)
                                                         tgCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGtttt                                  <  2:1455772/93‑1 (MQ=255)
                                                         tgCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGtttt                                  <  2:31305/93‑1 (MQ=255)
                                                          gCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTa                                 <  2:2579076/93‑1 (MQ=255)
                                                          gCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTa                                 <  2:6128268/93‑1 (MQ=255)
                                                          gCAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTa                                 >  1:1000409/1‑93 (MQ=255)
                                                           cAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTAc                                <  2:511927/93‑1 (MQ=255)
                                                           cAAATAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTAc                                <  2:5859494/93‑1 (MQ=255)
                                                             aaTAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCg                              <  1:5597589/93‑1 (MQ=255)
                                                              aTAGTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCgg                             <  1:3374598/93‑1 (MQ=255)
                                                                 gTACTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGaaa                          >  1:3721916/1‑93 (MQ=255)
                                                                   aCTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGaaaaa                        >  1:3176304/1‑93 (MQ=255)
                                                                    cTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAAt                       >  2:32341/1‑93 (MQ=255)
                                                                    cTGGCGTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAAt                       >  2:5027873/1‑93 (MQ=255)
                                                                         gTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAGTCCGGc                  <  1:5107601/93‑1 (MQ=255)
                                                                         gTAGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGc                  <  1:6117684/93‑1 (MQ=255)
                                                                           aGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTg                >  1:1670448/1‑93 (MQ=255)
                                                                           aGCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTg                <  2:3066281/93‑1 (MQ=255)
                                                                            gCAGAAGATGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTgg               >  1:2270237/1‑93 (MQ=255)
                                                                              agaagaTGAACACTGTGATCTCTCAAGTATCTTTAATATCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGt             <  2:746049/93‑1 (MQ=255)
                                                                              agaagaTGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTAGGt             <  2:196542/93‑1 (MQ=255)
                                                                              agaagaTGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGt             >  2:2931367/1‑93 (MQ=255)
                                                                               gaagaTGAACACTGTGATCTCTCATATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTa            >  2:1846857/1‑93 (MQ=255)
                                                                               gaagaTGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTa            >  2:548764/1‑93 (MQ=255)
                                                                               gaagaTGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCAGGCTGGGTa            >  2:516431/1‑93 (MQ=255)
                                                                                aagaTGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAg           <  1:1510434/93‑1 (MQ=255)
                                                                                    tGAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGgtgt       >  2:5241532/1‑93 (MQ=255)
                                                                                    tGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGgtgt       >  2:1170880/1‑93 (MQ=255)
                                                                                     gAACACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGtc      <  2:5593917/93‑1 (MQ=255)
                                                                                     gAACACTGTGATCTCTCAAATATCTTTAATAGAACCAATAACCAGCAAGGGCAGTTGGTTGTTTTACCGGAAAAATCCGGCTGGGAGGgtgcc      >  1:1063966/1‑91 (MQ=255)
                                                                                      aaCACTGTGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGtct     >  2:2260132/1‑93 (MQ=255)
                                                                                            |                                                                                         
ATTATTATTGGCGATGGCATGTTTGTGTCTGGCAAACATATCCTGGGCTACTGTTTGTGCAAATAGTACTGGCGTAGCAGAAGATGAACACTATGATCTCTCAAATATCTTTAATAGCACCAATAACCAGCCAGGGCAGATTGTTGTTTTACCGGAAAAATCCGGCTGGGTAGGTGTCTCAG  >  NC_000913/561552‑561733

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