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,5460AG99.7% 964.0 / ‑4.7 284intergenic (+4/‑7)sfmH/sfmFinvolved in fimbrial asembly/putative fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (109/174);  total (110/174)
Fisher's exact test for biased strand distribution p-value = 3.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.30e-01

Only 100 of 290 total aligned reads displayed.

CGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAA  >  NC_000913/562456‑562635
                                                                                          |                                                                                         
  cgcCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGga                                                                                       >  2:4479265/1‑93 (MQ=255)
    ccTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGGATTTACGCGTCGATTATCAATGAGGTGCGgaga                                                                                     >  1:2240857/1‑93 (MQ=255)
      tgcccTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTTGGGTATTTACGCTTCGATTATCAATGAGGTGCGGAGAAt                                                                                   <  2:7147766/90‑1 (MQ=255)
      tGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAAt                                                                                   <  1:8596353/93‑1 (MQ=255)
      tGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAAt                                                                                   <  2:1409740/93‑1 (MQ=255)
       ggCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATg                                                                                  <  1:4686019/93‑1 (MQ=255)
        gCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATga                                                                                 >  1:1008666/1‑93 (MQ=255)
        gCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATga                                                                                 <  1:917546/93‑1 (MQ=255)
         ccTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATgag                                                                                <  1:55029/93‑1 (MQ=255)
          cTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATgaga                                                                               <  1:7066410/93‑1 (MQ=255)
           tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATgggag                                                                              >  2:3926327/1‑92 (MQ=255)
           tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaa                                                                              >  1:7927462/1‑93 (MQ=255)
           tATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaa                                                                              >  1:7423741/1‑93 (MQ=255)
            aTCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaag                                                                             <  2:8370106/93‑1 (MQ=255)
            aTCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaag                                                                             >  2:7462166/1‑93 (MQ=255)
             tCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaaga                                                                            <  2:6791543/93‑1 (MQ=255)
             tCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGagaaga                                                                            >  1:3200620/1‑93 (MQ=255)
               aGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGAGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGt                                                                          <  1:2958825/93‑1 (MQ=255)
                gTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTa                                                                         <  2:4833861/93‑1 (MQ=255)
                gTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTa                                                                         <  2:2084681/93‑1 (MQ=255)
                gTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGGATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTa                                                                         >  2:5181677/1‑93 (MQ=255)
                 tATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTAc                                                                        <  1:2475917/93‑1 (MQ=255)
                 tATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTAc                                                                        <  2:7252986/93‑1 (MQ=255)
                 tATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTAc                                                                        >  2:4894927/1‑93 (MQ=255)
                  aTTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTAct                                                                       <  1:982597/93‑1 (MQ=255)
                  aTTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTAct                                                                       >  1:2320761/1‑93 (MQ=255)
                   ttACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActc                                                                      <  1:4755480/93‑1 (MQ=255)
                    tACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActct                                                                     <  2:4582211/93‑1 (MQ=255)
                    tACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTActct                                                                     >  1:7777232/1‑93 (MQ=255)
                      ccGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCttt                                                                   <  2:3308652/93‑1 (MQ=255)
                      ccGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCttt                                                                   >  2:7420149/1‑93 (MQ=255)
                      ccGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGAACTCttt                                                                   >  1:1966868/1‑93 (MQ=255)
                       cGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTa                                                                  <  1:8546698/93‑1 (MQ=255)
                         gTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGc                                                                >  1:2718500/1‑93 (MQ=255)
                            aaaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGt                                                             <  1:2080132/93‑1 (MQ=255)
                             aaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGtt                                                            >  2:3709583/1‑93 (MQ=255)
                             aaaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGtt                                                            >  1:5315757/1‑93 (MQ=255)
                              aaaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGGATGAGAAGAGTACTCTTTAGCTGttt                                                           >  1:6217455/1‑93 (MQ=255)
                               aaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTc                                                          <  2:1562381/93‑1 (MQ=255)
                               aaCCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTc                                                          >  2:7775872/1‑93 (MQ=255)
                                 ccgaccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTg                                                        >  2:2693219/1‑93 (MQ=255)
                                     ccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCggg                                                    >  2:8491587/1‑93 (MQ=255)
                                     ccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCggg                                                    >  1:4084155/1‑93 (MQ=255)
                                     ccgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCggg                                                    >  2:6475101/1‑93 (MQ=255)
                                      cgaAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGc                                                   <  1:8113135/93‑1 (MQ=255)
                                         aGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTAc                                                <  2:2845697/93‑1 (MQ=255)
                                         aGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTAc                                                <  2:3137356/93‑1 (MQ=255)
                                           ggCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtg                                              <  2:8619053/93‑1 (MQ=255)
                                             ccGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtgtg                                            <  2:2893236/93‑1 (MQ=255)
                                             ccGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACtgtg                                            <  2:8647121/93‑1 (MQ=255)
                                              cGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTgg                                           >  1:6049873/1‑93 (MQ=255)
                                              cGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTgg                                           <  2:6982749/93‑1 (MQ=255)
                                               gTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGATTGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGa                                          >  1:8174601/1‑93 (MQ=255)
                                                tttAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAg                                         <  1:616338/93‑1 (MQ=255)
                                                tttAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAg                                         >  1:2531402/1‑93 (MQ=255)
                                                  tAGCGCGGTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGtt                                       <  1:5160020/93‑1 (MQ=255)
                                                  tAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGtt                                       <  1:2918290/93‑1 (MQ=255)
                                                     cgcgcTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCa                                    <  1:14118/93‑1 (MQ=255)
                                                      gcgcTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAg                                   <  1:8190428/93‑1 (MQ=255)
                                                      gcgcTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAg                                   <  2:4208642/93‑1 (MQ=255)
                                                       cgcTGGGGTAATTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGt                                  <  1:7511518/93‑1 (MQ=255)
                                                        gcTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTg                                 <  1:7008632/93‑1 (MQ=255)
                                                         cTGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTgg                                >  2:4295443/1‑93 (MQ=255)
                                                          tGGGGTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACCGTGGATTTCCAGTGGa                               <  2:8084554/93‑1 (MQ=255)
                                                           ggggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGAt                              <  1:54155/93‑1 (MQ=255)
                                                           ggggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGAt                              <  1:4182120/93‑1 (MQ=255)
                                                           ggggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGAt                              >  1:3709095/1‑93 (MQ=255)
                                                           ggggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGAt                              <  2:3426220/93‑1 (MQ=255)
                                                            gggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATg                             >  1:338714/1‑93 (MQ=255)
                                                             ggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATgg                            <  1:188662/93‑1 (MQ=255)
                                                             ggTATTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATgg                            <  2:589079/93‑1 (MQ=255)
                                                                aTTTACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCa                         >  2:6156585/1‑93 (MQ=255)
                                                                 tttACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTGCGGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAg                        <  2:8456710/93‑1 (MQ=255)
                                                                  ttACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGt                       <  2:6395343/93‑1 (MQ=255)
                                                                  ttACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGt                       <  2:2247365/93‑1 (MQ=255)
                                                                   tACGCGTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGtt                      <  1:1669443/93‑1 (MQ=255)
                                                                      gcgTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAc                   >  1:4907772/1‑93 (MQ=255)
                                                                       cgTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAcc                  >  2:1285101/1‑93 (MQ=255)
                                                                       cgTCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGAcc                  <  1:2166306/93‑1 (MQ=255)
                                                                         tCGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCt                <  1:8456507/93‑1 (MQ=255)
                                                                          cGATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCtt               <  2:1511349/93‑1 (MQ=255)
                                                                           gATTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGGTGACCCttt              >  1:1220741/1‑93 (MQ=255)
                                                                            aTTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTa             <  2:5840610/93‑1 (MQ=255)
                                                                            aTTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTa             <  2:4047950/93‑1 (MQ=255)
                                                                            aTTATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAACTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTTACCCTTTa             >  2:5486478/1‑93 (MQ=255)
                                                                             ttATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAg            <  1:1132165/93‑1 (MQ=255)
                                                                             ttATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAg            <  2:1749192/93‑1 (MQ=255)
                                                                             ttATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAg            <  2:750145/93‑1 (MQ=255)
                                                                             ttATCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAg            <  1:6908057/93‑1 (MQ=255)
                                                                                tCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGaa         <  2:2642807/93‑1 (MQ=255)
                                                                                tCAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGaa         <  1:5099424/93‑1 (MQ=255)
                                                                                 cAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAAc        <  1:6009559/93‑1 (MQ=255)
                                                                                 cAATGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAAc        <  2:4159308/93‑1 (MQ=255)
                                                                                  aaTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACg       >  2:3010388/1‑93 (MQ=255)
                                                                                  aaTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACg       >  1:5733221/1‑93 (MQ=255)
                                                                                  aaTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACg       >  2:538604/1‑93 (MQ=255)
                                                                                   aTGAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGa      <  1:1229831/93‑1 (MQ=255)
                                                                                     gAGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGAtt    <  1:1608740/93‑1 (MQ=255)
                                                                                      aGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTa   <  2:7387952/93‑1 (MQ=255)
                                                                                      aGGTGCGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTa   <  1:6577497/93‑1 (MQ=255)
                                                                                          |                                                                                         
CGCGCCTGGCCTATCAGTATTACCGGTCAAAAACCGACCGAAGGGCCGTTTAGCGCGCTGGGGTATTTACGCGTCGATTATCAATGAGGTACGGAGAATGAGAAGAGTACTCTTTAGCTGTTTCTGCGGGCTACTGTGGAGTTCCAGTGGATGGGCAGTTGACCCTTTAGGAACGATTAA  >  NC_000913/562456‑562635

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