Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,548,150 A→G T137A (ACG→GCG)  fimF → minor component of type 1 fimbriae

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,548,1500AG99.1% 393.2 / ‑4.8 110T137A (ACG→GCG) fimFminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (60/49);  total (60/50)
Fisher's exact test for biased strand distribution p-value = 4.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 110 total aligned reads displayed.

TACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGG  >  NC_000913/4548056‑4548246
                                                                                              |                                                                                                
tACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGAcc                                                                                             <  1:2356473/100‑1 (MQ=255)
tACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGAcc                                                                                             >  2:2836760/1‑100 (MQ=255)
 aCGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTCATGCTCCATCGTCCGCGCTTTCGTGGGCGAccc                                                                                            <  1:1685084/100‑1 (MQ=255)
 aCGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGAccc                                                                                            <  1:2468994/100‑1 (MQ=255)
       tCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGc                                                                                      >  2:1900717/1‑100 (MQ=255)
        cAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGcc                                                                                     >  2:2811546/1‑100 (MQ=255)
        cAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGcc                                                                                     <  2:1293070/100‑1 (MQ=255)
        cAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGcc                                                                                     >  2:1606565/1‑100 (MQ=255)
         aGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCg                                                                                    <  1:331908/100‑1 (MQ=255)
         aGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCg                                                                                    >  2:1071962/1‑100 (MQ=255)
          gCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCgg                                                                                   >  2:1614037/1‑100 (MQ=255)
          gCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCgg                                                                                   <  2:287906/100‑1 (MQ=255)
           cGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCggg                                                                                  <  1:609378/100‑1 (MQ=255)
           cGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCggg                                                                                  <  2:2777283/100‑1 (MQ=255)
              cTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTaa                                                                               >  2:214034/1‑100 (MQ=255)
              cTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTaa                                                                               <  2:2372513/100‑1 (MQ=255)
              cTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTaa                                                                               <  2:2499720/100‑1 (MQ=255)
                 cGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAAcc                                                                            >  2:763116/1‑100 (MQ=255)
                   ggACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCaa                                                                          >  1:1589892/1‑100 (MQ=255)
                        gggAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACg                                                                     <  2:1228098/100‑1 (MQ=255)
                        gggAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACg                                                                     >  1:523484/1‑100 (MQ=255)
                         ggAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGc                                                                    >  1:1101947/1‑100 (MQ=255)
                         ggAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGc                                                                    <  2:2018943/100‑1 (MQ=255)
                          gAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCt                                                                   >  1:2153529/1‑100 (MQ=255)
                           aaTACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGc                                                                    >  2:1046163/1‑98 (MQ=255)
                             tACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGaa                                                                >  1:1815588/1‑100 (MQ=255)
                             tACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGaa                                                                <  1:31199/100‑1 (MQ=255)
                                aGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAAttt                                                             >  2:625915/1‑100 (MQ=255)
                                aGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAAttt                                                             >  2:625932/1‑100 (MQ=255)
                                  cTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAAttttt                                                           <  1:419018/100‑1 (MQ=255)
                                   ttCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTa                                                          <  1:717192/100‑1 (MQ=255)
                                   ttCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCATGACGCCGGGTAAAGCCAATACGCTGAATTTTTa                                                          >  1:2619774/1‑100 (MQ=255)
                                      tGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACg                                                        <  1:2002122/99‑1 (MQ=255)
                                       gAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGcc                                                      >  2:195702/1‑100 (MQ=255)
                                       gAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGcc                                                      >  2:1987479/1‑100 (MQ=255)
                                       gAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGcc                                                      >  2:2184617/1‑100 (MQ=255)
                                       gAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGcc                                                      >  1:661494/1‑100 (MQ=255)
                                       gAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGcc                                                      >  2:1315455/1‑100 (MQ=255)
                                         aTGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCg                                                    >  2:2403998/1‑100 (MQ=255)
                                          tGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCgg                                                   >  2:850722/1‑100 (MQ=255)
                                           gAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGc                                                  >  1:1878763/1‑100 (MQ=255)
                                            agcagcAAAATCAAATACCCCTTAATGCTCCCTCGTCTGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCt                                                 <  1:1954325/100‑1 (MQ=255)
                                            agcagcAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCt                                                 <  1:2776104/100‑1 (MQ=255)
                                             gcagcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTa                                                >  2:609309/1‑100 (MQ=255)
                                             gcagcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTa                                                >  2:2530743/1‑100 (MQ=255)
                                              cagcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTaa                                               >  2:1975206/1‑100 (MQ=255)
                                              cagcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTaa                                               >  1:168854/1‑100 (MQ=255)
                                               agcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAAt                                              >  1:1056543/1‑100 (MQ=255)
                                                gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATg                                             <  1:957591/100‑1 (MQ=255)
                                                gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATg                                             >  1:1361761/1‑100 (MQ=255)
                                                gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATg                                             >  1:612051/1‑100 (MQ=255)
                                                gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATg                                             >  2:2121995/1‑100 (MQ=255)
                                                     aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                        >  2:2266247/1‑100 (MQ=255)
                                                     aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                        >  2:2265863/1‑100 (MQ=255)
                                                     aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                        <  2:436937/100‑1 (MQ=255)
                                                     aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                        <  1:301239/100‑1 (MQ=255)
                                                     aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                        >  1:2407054/1‑100 (MQ=255)
                                                     aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTAGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                        >  2:134351/1‑100 (MQ=255)
                                                       cAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGacac                                      <  1:1140285/100‑1 (MQ=255)
                                                       cAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGacac                                      <  1:1139838/100‑1 (MQ=255)
                                                        aaaTACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGAcaca                                     <  2:1422588/100‑1 (MQ=255)
                                                             ccccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGc                                >  2:2768518/1‑100 (MQ=255)
                                                             ccccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGc                                >  1:1139775/1‑100 (MQ=255)
                                                              cccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAGCCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGcc                               <  1:284502/100‑1 (MQ=255)
                                                              cccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGcc                               <  1:1394735/100‑1 (MQ=255)
                                                              cccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGcc                               <  2:1580757/100‑1 (MQ=255)
                                                              cccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGcc                               <  1:2836722/100‑1 (MQ=255)
                                                              cccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGcc                               <  2:2081504/100‑1 (MQ=255)
                                                               ccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCt                              >  2:1848662/1‑100 (MQ=255)
                                                               cattAATGCTCCATCGTCCGCGCTTTCGTGGGAGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCt                              >  2:2756698/3‑100 (MQ=255)
                                                                cTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTg                             <  2:2194488/100‑1 (MQ=255)
                                                                  tAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTc                           <  1:1487695/100‑1 (MQ=255)
                                                                   aaTGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCa                          >  2:1403509/1‑100 (MQ=255)
                                                                        tCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCg                     >  1:304794/1‑100 (MQ=255)
                                                                        tCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCg                     >  1:304971/1‑100 (MQ=255)
                                                                          cATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCATATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                   <  2:661500/100‑1 (MQ=255)
                                                                          cATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                   <  2:2865330/100‑1 (MQ=255)
                                                                          cATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                   <  2:2863791/100‑1 (MQ=255)
                                                                          cATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                   <  2:565999/100‑1 (MQ=255)
                                                                          cATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                   >  1:859358/1‑100 (MQ=255)
                                                                          cATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                   <  1:1892798/100‑1 (MQ=255)
                                                                          cATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                   <  2:1945271/100‑1 (MQ=255)
                                                                                ccGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCatat             >  1:1651737/1‑100 (MQ=255)
                                                                                  gcgcTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCa           >  1:67660/1‑100 (MQ=255)
                                                                                  gcgcTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCa           >  2:2775601/1‑100 (MQ=255)
                                                                                    gcTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAAt         <  2:136881/100‑1 (MQ=255)
                                                                                     cTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATg        <  1:2215547/100‑1 (MQ=255)
                                                                                     cTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATg        >  2:600431/1‑100 (MQ=255)
                                                                                     cTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATg        >  1:2836948/1‑100 (MQ=255)
                                                                                      tttCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGc       <  2:2577473/100‑1 (MQ=255)
                                                                                       ttCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGcc      >  2:2588693/1‑100 (MQ=255)
                                                                                       ttCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGcc      >  2:1587867/1‑100 (MQ=255)
                                                                                       ttCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGcc      >  1:382705/1‑100 (MQ=255)
                                                                                         cGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCAc    <  1:807442/100‑1 (MQ=255)
                                                                                         cGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCAc    <  1:807441/100‑1 (MQ=255)
                                                                                         cGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCAc    <  2:804150/100‑1 (MQ=255)
                                                                                          gTGGGCGACCTTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACACGTGCCTGTCACTGCGGGGCATATCAATGCCACg   <  1:1900691/100‑1 (MQ=255)
                                                                                          gTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATTCCAc    <  2:884353/99‑1 (MQ=255)
                                                                                          gTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACg   >  1:1885915/1‑100 (MQ=255)
                                                                                           tGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACAAAGGTGCCTGTCACTGCGGGGCATATCAATGCGACgg  >  1:1286068/1‑100 (MQ=255)
                                                                                              |                                                                                                
TACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGG  >  NC_000913/4548056‑4548246

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