Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,548,150 A→G 25.8% 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,1500AG25.8% 220.8 / 95.9 144T137A (ACG→GCG) fimFminor component of type 1 fimbriae
Reads supporting (aligned to +/- strand):  ref base A (68/38);  new base G (19/18);  total (87/57)
Fisher's exact test for biased strand distribution p-value = 1.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.17e-01

Only 100 of 148 total aligned reads displayed.

AAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGC  >  NC_000913/4548052‑4548247
                                                                                                  |                                                                                                 
aaaaTACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGgac                                                                                                  <  2:2985423/100‑1 (MQ=255)
aaaaTACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGgac                                                                                                  <  2:1861402/100‑1 (MQ=255)
 aaaTACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGgacg                                                                                                 <  1:3745856/100‑1 (MQ=255)
  aaTACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGggcga                                                                                                >  2:2153659/1‑100 (MQ=255)
    tACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGAcc                                                                                              >  2:1606665/1‑100 (MQ=255)
     aCGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGAccc                                                                                             <  1:2984970/100‑1 (MQ=255)
     aCGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGAccc                                                                                             <  1:3370822/100‑1 (MQ=255)
     aCGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGAccc                                                                                             <  2:1360510/100‑1 (MQ=255)
     aCGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGAccc                                                                                             <  1:3727956/100‑1 (MQ=255)
      cGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCt                                                                                            <  2:1323948/100‑1 (MQ=255)
      cGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCt                                                                                            >  2:369316/1‑100 (MQ=255)
        gtgtCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGa                                                                                          >  2:686473/1‑100 (MQ=255)
         tgtCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGAc                                                                                         <  1:3370982/100‑1 (MQ=255)
          gtCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACg                                                                                        >  1:724556/1‑100 (MQ=255)
          gtCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACg                                                                                        <  1:1731148/100‑1 (MQ=255)
            cAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGcc                                                                                      >  1:1813935/1‑100 (MQ=255)
            cAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGcc                                                                                      >  1:965367/1‑100 (MQ=255)
            cAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGcc                                                                                      >  2:481358/1‑100 (MQ=255)
              gCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCgg                                                                                    >  1:1533749/1‑100 (MQ=255)
                ggCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGt                                                                                  >  1:3503471/1‑100 (MQ=255)
                ggCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACACTGACGCCGGGt                                                                                  >  1:1268885/1‑100 (MQ=255)
                  cTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTaa                                                                                >  1:354990/1‑100 (MQ=255)
                  cTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTaa                                                                                >  2:902879/1‑100 (MQ=255)
                     cggaattGGGAATACGGCCTCTCTATGAGCAGCAAAATCAATTACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAAc                                                                              <  2:2396462/93‑1 (MQ=255)
                      gggACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCa                                                                            >  1:2498867/1‑100 (MQ=255)
                      gggACTGGGAATACAGCGTCTGAATGAGCAGCAAAATCAAATACGCCTTAATGCTCCATCGTCCGCGCGTTCGTGGACGACCATGACGCCAGGTAAACCa                                                                            >  2:1980645/1‑100 (MQ=255)
                        gACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAAcc                                                                             <  2:432021/97‑1 (MQ=255)
                         aCTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAAt                                                                         <  2:1015917/100‑1 (MQ=255)
                         aCTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAAt                                                                         <  2:1551002/100‑1 (MQ=255)
                          cTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATa                                                                        >  2:3817431/1‑100 (MQ=255)
                          cTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATa                                                                        >  2:1762912/1‑100 (MQ=255)
                          cTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCAGGTAAACCAAATa                                                                        >  1:354192/1‑100 (MQ=255)
                            gggAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACg                                                                      >  1:3554134/1‑100 (MQ=255)
                            gggAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACg                                                                      >  1:3013982/1‑100 (MQ=255)
                             ggAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGc                                                                     >  2:549263/1‑100 (MQ=255)
                                 tACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGaa                                                                 >  1:1139439/1‑100 (MQ=255)
                                 tACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGaa                                                                 >  1:2851889/1‑100 (MQ=255)
                                 tACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGaa                                                                 >  1:80141/1‑100 (MQ=255)
                                   cAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAAtt                                                               <  2:3391280/100‑1 (MQ=255)
                                   cAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAAtt                                                               >  1:715499/1‑100 (MQ=255)
                                   cAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACACTGAAtt                                                               >  1:1267173/1‑100 (MQ=255)
                                    aGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAAttt                                                              <  2:3481248/100‑1 (MQ=255)
                                     gCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAAtttt                                                             >  2:2018408/1‑100 (MQ=255)
                                      cTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAAttttt                                                            >  2:1449187/1‑100 (MQ=255)
                                       ttCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTa                                                           <  1:1876158/100‑1 (MQ=255)
                                         cTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACg                                                         >  2:1231043/1‑100 (MQ=255)
                                          tGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGc                                                        >  1:1003717/1‑100 (MQ=255)
                                          tGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGc                                                        >  1:2238528/1‑100 (MQ=255)
                                           gAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGcc                                                       >  1:800015/1‑100 (MQ=255)
                                             aTGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCg                                                     >  2:1065586/1‑100 (MQ=255)
                                              tGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCgg                                                    >  1:2901570/1‑100 (MQ=255)
                                               gAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGc                                                   <  1:1411777/100‑1 (MQ=255)
                                               gAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGc                                                   <  1:323254/100‑1 (MQ=255)
                                                 gcagcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTa                                                 >  1:2783707/1‑100 (MQ=255)
                                                  cagcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAACTTTTACGCCCGGCTaa                                                >  2:3723089/1‑100 (MQ=255)
                                                   agcaAAATTAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAAt                                               >  2:2867165/1‑100 (MQ=255)
                                                   agcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAAt                                               >  2:3195128/1‑100 (MQ=255)
                                                   agcaAAATCAAATACCCCTTAATGCTACATCCTCCGCGCTCTCGTGGACGACCTTGAAGCGGGGTAAACGAAATAGGTTGAATTTTTACGCCCGGATAAt                                               >  2:1990719/1‑100 (MQ=255)
                                                    gcaAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATg                                              >  2:424407/1‑100 (MQ=255)
                                                     caAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATg                                              >  2:346422/1‑99 (MQ=255)
                                                        aaTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGa                                          <  2:3288197/100‑1 (MQ=255)
                                                         aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                         >  1:1954476/1‑100 (MQ=255)
                                                         aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                         <  1:2695409/100‑1 (MQ=255)
                                                         aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGCCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                         <  2:1951332/100‑1 (MQ=255)
                                                         aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                         <  2:136550/100‑1 (MQ=255)
                                                         aTCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGac                                         >  1:2467278/1‑100 (MQ=255)
                                                          tctagtACCCCATAATCCCCCATCTTTCCCGCTTTCGTGGACGACCCTGACGCCGCGTAAACCAAATACGCTTAATATTAACGCACGCCTAATGGc                                            <  2:3809164/91‑1 (MQ=255)
                                                              aTACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAgg                                    <  1:1014242/100‑1 (MQ=255)
                                                              aTACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAgg                                    >  1:1777060/1‑100 (MQ=255)
                                                                aCCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTg                                  >  2:249834/1‑100 (MQ=255)
                                                                aCCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTg                                  <  1:1718737/100‑1 (MQ=255)
                                                                   ccTTAATGCTCCATCGTCCGCGCTTTCGTGGGCGACCCTGACGCCGGGTCAACCAACTACGCTGAATTTTTACGCCCGGCTAATGGCGACACAAGTGCCt                               >  2:3291945/1‑100 (MQ=255)
                                                                   ccTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCt                               >  2:3050288/1‑100 (MQ=255)
                                                                   ccTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCt                               >  1:2263847/1‑100 (MQ=255)
                                                                    cTTAATGCTCCATCGTCCGCGCTTTCGTTGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTg                              >  2:3819311/1‑100 (MQ=255)
                                                                      tAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGt                             <  1:1112985/99‑1 (MQ=255)
                                                                          gCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTg                        >  1:2801159/1‑100 (MQ=255)
                                                                          gCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTg                        >  1:1024954/1‑100 (MQ=255)
                                                                           cTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGc                       <  1:1189664/100‑1 (MQ=255)
                                                                            tCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCg                      >  2:2764496/1‑100 (MQ=255)
                                                                              cATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                    <  1:3208116/100‑1 (MQ=255)
                                                                              cATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                    <  1:2769197/100‑1 (MQ=255)
                                                                              cATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCggg                    <  1:1558618/100‑1 (MQ=255)
                                                                                tCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGc                  >  1:1218247/1‑100 (MQ=255)
                                                                                  gTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAACCCAAATCCGCTGAATTTTTACACCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCa                 >  2:2173676/1‑99 (MQ=255)
                                                                                  gTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCat                >  2:924826/1‑100 (MQ=255)
                                                                                   tCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCata               >  2:3370653/1‑100 (MQ=255)
                                                                                   tCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCata               >  1:1850307/1‑100 (MQ=255)
                                                                                     cgcgCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATc             <  1:2153673/100‑1 (MQ=255)
                                                                                     cgcgCTTTCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATc             <  2:2717360/100‑1 (MQ=255)
                                                                                         cTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCACTg         >  2:1744957/1‑100 (MQ=255)
                                                                                         cTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATg         <  2:1102855/100‑1 (MQ=255)
                                                                                          tttCGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGc        >  1:202367/1‑100 (MQ=255)
                                                                                           ttCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACCCAGGTGCCTGTCACTGCGGGGCATATCAATGcc       <  2:1075846/100‑1 (MQ=255)
                                                                                            tCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCa      <  2:666320/100‑1 (MQ=255)
                                                                                             cGTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCAc     <  1:196216/100‑1 (MQ=255)
                                                                                              gTGGGCGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACg    <  2:3791607/100‑1 (MQ=255)
                                                                                               tGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACgg   <  1:1345506/100‑1 (MQ=255)
                                                                                                ggACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGc  <  1:2531276/100‑1 (MQ=255)
                                                                                                ggACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGc  >  2:2332224/1‑100 (MQ=255)
                                                                                                  |                                                                                                 
AAAATACGGTGTCAGCGGCTTCGGGACTGGGAATACAGCTTCTGAATGAGCAGCAAAATCAAATACCCCTTAATGCTCCATCGTCCGCGCTTTCGTGGACGACCCTGACGCCGGGTAAACCAAATACGCTGAATTTTTACGCCCGGCTAATGGCGACACAGGTGCCTGTCACTGCGGGGCATATCAATGCCACGGC  >  NC_000913/4548052‑4548247

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