Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 562,546 A→G intergenic (+4/‑7) sfmH → / → sfmF involved in fimbrial asembly/putative fimbrial‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913562,5460AG84.5% 576.1 / 97.7 258intergenic (+4/‑7)sfmH/sfmFinvolved in fimbrial asembly/putative fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base A (19/21);  new base G (99/119);  total (118/140)
Fisher's exact test for biased strand distribution p-value = 8.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

Only 100 of 258 total aligned reads displayed.

TCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATA  >  NC_000913/562455‑562637
                                                                                            |                                                                                           
tCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTAc                                                                                            <  2:5281518/93‑1 (MQ=255)
 cgcgCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCg                                                                                           >  2:2553727/1‑93 (MQ=255)
 cgcgCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACg                                                                                           <  2:753403/93‑1 (MQ=255)
    gcCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGgag                                                                                        >  2:6059311/1‑93 (MQ=255)
     ccTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGgaga                                                                                       >  1:1129117/1‑93 (MQ=255)
      cTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGaa                                                                                      <  1:1254200/93‑1 (MQ=255)
      cTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCACTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGaa                                                                                      <  1:19820/93‑1 (MQ=255)
       tGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAAt                                                                                     <  2:168168/93‑1 (MQ=255)
          ccTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATgag                                                                                  <  1:1571764/93‑1 (MQ=255)
            tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaa                                                                                <  1:2456086/93‑1 (MQ=255)
             aTCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaag                                                                               <  2:3492992/93‑1 (MQ=255)
             aTCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGagaag                                                                               >  2:3258213/1‑93 (MQ=255)
              tCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaaga                                                                              >  1:5783349/1‑93 (MQ=255)
              tCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaaga                                                                              <  1:4341403/93‑1 (MQ=255)
               taGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAagag                                                                             >  1:2371048/2‑93 (MQ=255)
               cAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAagag                                                                             >  2:7151787/1‑93 (MQ=255)
               cAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAagag                                                                             <  2:931389/93‑1 (MQ=255)
                aGTATTACCGGTCAAAAACCGACCGAAGGGCGGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGt                                                                            >  1:4488165/1‑93 (MQ=255)
                aGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGt                                                                            <  1:7133012/93‑1 (MQ=255)
                 gTATTACCGGTCAAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGt                                                                            <  1:2356679/93‑1 (MQ=255)
                 gTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTa                                                                           >  1:2951261/1‑93 (MQ=255)
                  tATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTAc                                                                          >  2:417855/1‑93 (MQ=255)
                    ttACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActc                                                                        >  2:6178623/1‑93 (MQ=255)
                    ttACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActc                                                                        >  2:2167671/1‑93 (MQ=255)
                    ttACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActc                                                                        >  2:3550826/1‑93 (MQ=255)
                    ttACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActc                                                                        >  2:788803/1‑93 (MQ=255)
                     tACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActct                                                                       <  1:3295045/93‑1 (MQ=255)
                     tACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActct                                                                       <  2:6944594/93‑1 (MQ=255)
                     tACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActct                                                                       <  1:6633269/93‑1 (MQ=255)
                          gTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGc                                                                  <  1:3963944/93‑1 (MQ=255)
                          gTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGc                                                                  >  2:4404303/1‑93 (MQ=255)
                            cAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTg                                                                <  1:6789636/93‑1 (MQ=255)
                              aaaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGt                                                               <  2:5154714/93‑1 (MQ=255)
                              aaaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGt                                                               >  2:3781737/1‑93 (MQ=255)
                              aaaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGt                                                               >  1:415262/1‑93 (MQ=255)
                               aaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGtt                                                              <  1:535967/93‑1 (MQ=255)
                                aaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGttt                                                             >  1:6528010/1‑93 (MQ=255)
                                 aaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTc                                                            <  2:6212325/93‑1 (MQ=255)
                                 aaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTc                                                            <  1:1002085/93‑1 (MQ=255)
                                 aaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTc                                                            >  1:971798/1‑93 (MQ=255)
                                  accgaccgAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCt                                                           >  2:6696683/1‑93 (MQ=255)
                                     gaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCg                                                        >  1:6267174/1‑93 (MQ=255)
                                     gaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCg                                                        >  1:1875121/1‑93 (MQ=255)
                                     gaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCg                                                        <  2:1350946/93‑1 (MQ=255)
                                     gaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCg                                                        >  2:4692852/1‑93 (MQ=255)
                                       ccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGGGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCggg                                                      >  1:3387036/1‑93 (MQ=255)
                                       ccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCggg                                                      >  1:3813706/1‑93 (MQ=255)
                                        cgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGc                                                     <  1:5659859/93‑1 (MQ=255)
                                         gaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCt                                                    >  1:1437603/1‑93 (MQ=255)
                                           aGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTAc                                                  >  1:6983686/1‑93 (MQ=255)
                                            gggCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACt                                                 >  1:1416673/1‑93 (MQ=255)
                                             ggCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtg                                                <  1:885055/93‑1 (MQ=255)
                                             ggCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtg                                                >  1:5136076/1‑93 (MQ=255)
                                             ggCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtg                                                <  1:5586986/93‑1 (MQ=255)
                                              gCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtgt                                               <  1:1824484/93‑1 (MQ=255)
                                              gCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtgt                                               <  2:1578400/93‑1 (MQ=255)
                                               ccGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtgtg                                              <  1:3120357/93‑1 (MQ=255)
                                               ccGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtgtg                                              <  2:831341/93‑1 (MQ=255)
                                               ccGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtgtg                                              >  1:3590098/1‑93 (MQ=255)
                                                 gTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGa                                            >  2:2404488/1‑93 (MQ=255)
                                                 gTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGa                                            >  2:5547999/1‑93 (MQ=255)
                                                 gTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGa                                            >  1:6010511/1‑93 (MQ=255)
                                                        gcgcTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAg                                     <  1:4639272/93‑1 (MQ=255)
                                                         cgcTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGt                                    <  1:3086261/93‑1 (MQ=255)
                                                           cTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTgg                                  >  1:5051485/1‑93 (MQ=255)
                                                            tGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGa                                 <  2:6678327/93‑1 (MQ=255)
                                                            tGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGa                                 <  2:6370636/93‑1 (MQ=255)
                                                            tGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGa                                 <  2:5290669/93‑1 (MQ=255)
                                                            tGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGa                                 <  1:179706/93‑1 (MQ=255)
                                                             ggggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGAt                                >  2:4357034/1‑93 (MQ=255)
                                                              gggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATg                               >  1:6353687/1‑93 (MQ=255)
                                                               ggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATgg                              <  2:499251/93‑1 (MQ=255)
                                                                 tATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGc                            <  2:1158493/93‑1 (MQ=255)
                                                                  aTTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCa                           <  2:2105212/93‑1 (MQ=255)
                                                                  aTTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCa                           >  2:198990/1‑93 (MQ=255)
                                                                  aTTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCa                           <  1:772218/93‑1 (MQ=255)
                                                                        gcgTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAc                     >  2:3138300/1‑93 (MQ=255)
                                                                        gcgTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAc                     <  2:4493537/93‑1 (MQ=255)
                                                                         cgTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAcc                    <  2:3231425/93‑1 (MQ=255)
                                                                         cgTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAcc                    >  1:7178963/1‑93 (MQ=255)
                                                                          gTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAccc                   >  1:1924752/1‑93 (MQ=255)
                                                                          gTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAccc                   <  2:2131466/93‑1 (MQ=255)
                                                                           tCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCt                  <  2:4092077/93‑1 (MQ=255)
                                                                           tCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCt                  <  1:1311028/93‑1 (MQ=255)
                                                                              aTTCTCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTa               <  2:4886898/93‑1 (MQ=255)
                                                                              aTTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTa               <  1:6598803/93‑1 (MQ=255)
                                                                              aTTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTa               <  1:7459342/93‑1 (MQ=255)
                                                                              aTTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTa               <  2:6544258/93‑1 (MQ=255)
                                                                               ttATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAg              >  2:1841229/1‑93 (MQ=255)
                                                                               ttATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAg              <  2:5622412/93‑1 (MQ=255)
                                                                               ttATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAg              <  1:7027425/93‑1 (MQ=255)
                                                                                 aTCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGa            <  1:4724439/93‑1 (MQ=255)
                                                                                  tCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGaa           >  1:5846279/1‑93 (MQ=255)
                                                                                    aaTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACg         >  1:655672/1‑93 (MQ=255)
                                                                                     aTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGa        <  1:3678182/93‑1 (MQ=255)
                                                                                     aTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGa        <  2:3024591/93‑1 (MQ=255)
                                                                                     aTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGa        >  2:4604785/1‑93 (MQ=255)
                                                                                      tGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGAt       >  2:2732365/1‑93 (MQ=255)
                                                                                        aGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTa     <  1:2058086/93‑1 (MQ=255)
                                                                                         ggTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACATTTaa    <  2:2870970/93‑1 (MQ=255)
                                                                                            |                                                                                           
TCGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAATA  >  NC_000913/562455‑562637

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