Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,545,134 A→G T13T (ACA→ACG fimD → type I fimbriae usher protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,545,1340AG98.9% 315.4 / ‑4.5 88T13T (ACA→ACGfimDtype I fimbriae usher protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (43/44);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.10e-01

TTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGGATGATC  >  NC_000913/4544998‑4545267
                                                                                                                                        |                                                                                                                                     
ttACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCaa                                                                                                                                    <  1:680743/140‑1 (MQ=255)
   cccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGc                                                                                                                                 <  1:776153/140‑1 (MQ=255)
   cccccAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGc                                                                                                                                 <  1:891227/140‑1 (MQ=255)
        aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAGGACTTTACCAGCGAAACACGCAATGCTTGCa                                                                                                                            <  2:701052/140‑1 (MQ=255)
        aaaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGGGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCa                                                                                                                            >  2:338503/1‑140 (MQ=255)
         aaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCa                                                                                                                            >  1:187476/1‑139 (MQ=255)
         aaaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCa                                                                                                                            <  2:187476/139‑1 (MQ=255)
          aaTGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTGCCAGCGAATCACGCAATGCTTGCata                                                                                                                          <  2:278096/140‑1 (MQ=255)
              aCGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCg                                                                                                                      >  1:1111648/1‑140 (MQ=255)
                  gCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAg                                                                                                                  >  1:1065777/1‑140 (MQ=255)
                   cGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGc                                                                                                                 >  1:1273248/1‑140 (MQ=255)
                    gTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCa                                                                                                                >  1:1004335/1‑140 (MQ=255)
                     tAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCAt                                                                                                               >  2:635845/1‑140 (MQ=255)
                            aTAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTgg                                                                                                        <  1:1263080/140‑1 (MQ=255)
                             tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGc                                                                                                       <  2:1299855/140‑1 (MQ=255)
                             tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGc                                                                                                       <  2:1065777/140‑1 (MQ=255)
                             tAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGc                                                                                                       <  1:50009/140‑1 (MQ=255)
                               aCGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTg                                                                                                     <  1:736192/140‑1 (MQ=255)
                                cGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTgg                                                                                                    >  1:87248/1‑140 (MQ=255)
                                cGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTgg                                                                                                    >  2:1208627/1‑140 (MQ=255)
                                  cAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTAGtt                                                                                                  >  2:267141/1‑140 (MQ=255)
                                     ggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGttttt                                                                                               >  2:735847/1‑140 (MQ=255)
                                     ggggAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGttttt                                                                                               >  1:49933/1‑140 (MQ=255)
                                          aTTTTTCGCCTGAATTAAAAGAATTTGCTGCCGGGTGATTTTAAGCCGGAGGAATAAATTCATTTTCTAATTTTAGGCTTTTCCCACCAAACCCCCAATGCCTGGATTTTTGGAAGCATCGGTTGGCTGGTTTTTTTGTc                                                                                          <  2:348956/140‑1 (MQ=255)
                                                cGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTc                                                                                    <  2:1004335/140‑1 (MQ=255)
                                                     gAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGt                                                                               <  2:1351328/140‑1 (MQ=255)
                                                      aaTAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTc                                                                              >  2:867680/1‑140 (MQ=255)
                                                                aTTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtt                                                                    >  2:827470/1‑140 (MQ=255)
                                                                aTTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtt                                                                    >  1:763593/1‑140 (MQ=255)
                                                                aTTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtt                                                                    >  2:358820/1‑140 (MQ=255)
                                                                aTTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtt                                                                    >  2:395526/1‑140 (MQ=255)
                                                                aTTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtt                                                                    >  1:334334/1‑140 (MQ=255)
                                                                aTTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtt                                                                    >  1:605894/1‑140 (MQ=255)
                                                                aTTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtt                                                                    >  2:296016/1‑140 (MQ=255)
                                                                 ttGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCttt                                                                   >  2:770769/1‑140 (MQ=255)
                                                                  tGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCtttt                                                                  <  1:827470/140‑1 (MQ=255)
                                                                   gACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTg                                                                 <  2:1337984/140‑1 (MQ=255)
                                                                    aCTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGc                                                                >  2:567974/1‑140 (MQ=255)
                                                                      tGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCg                                                              >  1:191506/1‑140 (MQ=255)
                                                                         cGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGcac                                                           >  2:287967/1‑140 (MQ=255)
                                                                         cGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGcac                                                           <  1:635845/140‑1 (MQ=255)
                                                                          gggTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGcaca                                                          <  1:735847/140‑1 (MQ=255)
                                                                          gggTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGcaca                                                          <  2:1373347/140‑1 (MQ=255)
                                                                           ggTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAg                                                         <  1:296016/140‑1 (MQ=255)
                                                                            gTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAgg                                                        <  2:763593/140‑1 (MQ=255)
                                                                              gATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCa                                                      <  1:567974/140‑1 (MQ=255)
                                                                              gATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCa                                                      <  1:337438/140‑1 (MQ=255)
                                                                                        cGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCAt                                            <  2:24703/140‑1 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGcac                                                           <  2:1382072/123‑1 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGcac                                                           >  1:1382072/1‑123 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGca                                                            >  1:313076/1‑122 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGca                                                            <  2:313076/122‑1 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCt                                                   <  1:307284/131‑1 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCt                                                   >  2:307284/1‑131 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCt                                          >  1:1360933/1‑140 (MQ=255)
                                                                                          gaggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCt                                          >  1:560953/1‑140 (MQ=255)
                                                                                           aggaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTg                                         >  2:1277356/1‑140 (MQ=255)
                                                                                             gaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGcc                                       <  1:770769/140‑1 (MQ=255)
                                                                                             gaATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGcc                                       <  1:395526/140‑1 (MQ=255)
                                                                                              aataatGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCg                                      <  2:891440/140‑1 (MQ=255)
                                                                                                taatGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGAc                                    <  2:334334/140‑1 (MQ=255)
                                                                                                taatGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGAc                                    <  2:1111648/140‑1 (MQ=255)
                                                                                                  atGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACct                                  <  1:924833/140‑1 (MQ=255)
                                                                                                  atGTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACct                                  <  2:802220/140‑1 (MQ=255)
                                                                                                    gTCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACctct                                >  1:913120/1‑140 (MQ=255)
                                                                                                     tCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGAc                                    <  1:698165/135‑1 (MQ=255)
                                                                                                     tCATATCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGAc                                    >  2:698165/1‑135 (MQ=255)
                                                                                                       atatCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTAtt                             >  1:77773/1‑140 (MQ=255)
                                                                                                       atatCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTAtt                             >  1:139248/1‑140 (MQ=255)
                                                                                                        tatCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCg                                      >  1:1376175/1‑130 (MQ=255)
                                                                                                        tatCTGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCg                                      <  2:1376175/130‑1 (MQ=255)
                                                                                                            tGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAAt                        <  2:412278/140‑1 (MQ=255)
                                                                                                            tGAATTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAAt                        >  2:591124/1‑140 (MQ=255)
                                                                                                              aaTTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAAt                        <  2:79532/138‑1 (MQ=255)
                                                                                                              aaTTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAAt                        >  1:79532/1‑138 (MQ=255)
                                                                                                               aTTTAAGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCcg                     <  2:77773/140‑1 (MQ=255)
                                                                                                                    aGACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCtt                >  1:1142253/1‑140 (MQ=255)
                                                                                                                     gACTTTACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCttt               >  2:693080/1‑140 (MQ=255)
                                                                                                                         ttACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAg           >  1:540094/1‑140 (MQ=255)
                                                                                                                          tACCAGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGc          >  2:365902/1‑140 (MQ=255)
                                                                                                                          tACCAGCGAAACACGCAATGCATGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGc          <  2:1273248/140‑1 (MQ=255)
                                                                                                                              aGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGTCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGgat      <  2:265388/140‑1 (MQ=255)
                                                                                                                              aGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGgat      <  2:605894/140‑1 (MQ=255)
                                                                                                                              aGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGgat      <  2:540094/140‑1 (MQ=255)
                                                                                                                              aGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGgat      <  2:49933/140‑1 (MQ=255)
                                                                                                                              aGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGgat      <  2:191506/140‑1 (MQ=255)
                                                                                                                              aGCGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGgat      <  1:693080/140‑1 (MQ=255)
                                                                                                                                cGAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGgatga    >  2:1015842/1‑140 (MQ=255)
                                                                                                                                 gAAACACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCCCCTTTGTCATCTGCCGAGCTCTATTTTAATCCGCGCTTTTTAGCGgatgat   <  2:1142253/140‑1 (MQ=255)
                                                                                                                                  aaaCACGCAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGGATGATc  >  2:121340/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                     
TTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAGAATTGACTGCCGGGTGATTTTAAGCCGGAGGAATAATGTCATATCTGAATTTAAGACTTTACCAGCGAAACACACAATGCTTGCATATTCGTAAGCATCGTTTGGCTGGTTTTTTTGTCCGACTCGTTGTCGCCTGTGCTTTTGCCGCACAGGCACCTTTGTCATCTGCCGACCTCTATTTTAATCCGCGCTTTTTAGCGGATGATC  >  NC_000913/4544998‑4545267

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