Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,289,828 G→A V283I (GTA→ATA)  yraK → putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,289,8280GA98.1% 737.3 / ‑1.2 212V283I (GTA→ATA) yraKputative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (95/113);  total (99/113)
Fisher's exact test for biased strand distribution p-value = 4.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 230 total aligned reads displayed.

ACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGA  >  NC_000913/3289738‑3289917
                                                                                          |                                                                                         
tcTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGata                                                                                         <  2:3538749/92‑1 (MQ=255)
 cTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGatag                                                                                        >  2:851800/1‑93 (MQ=255)
    cAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGAtagtag                                                                                     >  2:3438662/1‑93 (MQ=255)
    cAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGAtagtag                                                                                     >  1:2618456/1‑93 (MQ=255)
      gTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGAc                                                                                   <  2:3803370/93‑1 (MQ=255)
        aGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCa                                                                                 <  2:145232/93‑1 (MQ=255)
        aGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCa                                                                                 <  1:2064243/93‑1 (MQ=255)
        aGACTTTGATATCACTTATGACTGTGGTGATACTCCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCa                                                                                 <  2:424214/93‑1 (MQ=255)
          aCTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGt                                                                               >  1:2114883/1‑93 (MQ=255)
             ttGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACa                                                                            <  2:1958149/93‑1 (MQ=255)
               gATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAAc                                                                          <  2:947410/93‑1 (MQ=255)
               gATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAAc                                                                          <  1:735088/93‑1 (MQ=255)
                atatCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAAcc                                                                         >  1:930530/1‑93 (MQ=255)
                 tatCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCt                                                                        >  1:3567307/1‑93 (MQ=255)
                    cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGt                                                                     <  1:3714895/93‑1 (MQ=255)
                    cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGt                                                                     <  1:4682667/93‑1 (MQ=255)
                    cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGt                                                                     <  2:2682297/93‑1 (MQ=255)
                    cACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGt                                                                     <  1:610345/93‑1 (MQ=255)
                     aCTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTc                                                                    >  2:636327/1‑93 (MQ=255)
                      cTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCg                                                                   >  2:3268767/1‑93 (MQ=255)
                      cTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCg                                                                   >  2:3005501/1‑93 (MQ=255)
                       ttATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGc                                                                  >  2:3639705/1‑93 (MQ=255)
                        tATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGcc                                                                 <  1:2828013/93‑1 (MQ=255)
                            aCTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGc                                                             >  2:3390052/1‑93 (MQ=255)
                            aCTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGc                                                             <  1:437201/93‑1 (MQ=255)
                            aCTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGc                                                             <  1:741624/93‑1 (MQ=255)
                            aCTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGc                                                             <  1:1573621/93‑1 (MQ=255)
                            aCTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGc                                                             <  1:2428428/93‑1 (MQ=255)
                                tggtgATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCgaag                                                         <  2:2564396/93‑1 (MQ=255)
                                tggtgATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCgaag                                                         <  1:2921144/93‑1 (MQ=255)
                                    gATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGt                                                     <  1:2113042/93‑1 (MQ=255)
                                    gATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACATGGTCGCCAGGCGAAGAAGt                                                     >  2:1883213/1‑93 (MQ=255)
                                    gATACATCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGt                                                     >  2:4709058/1‑93 (MQ=255)
                                     aTACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGtt                                                    <  2:3106821/93‑1 (MQ=255)
                                      tACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTc                                                   <  1:3925142/93‑1 (MQ=255)
                                      tACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTc                                                   <  1:1728143/93‑1 (MQ=255)
                                       aCTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGACGAAGTTCa                                                  >  2:173310/1‑93 (MQ=255)
                                       aCTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCa                                                  >  2:304733/1‑93 (MQ=255)
                                         ttCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGa                                                >  1:2992345/1‑93 (MQ=255)
                                         ttCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGa                                                >  2:4599886/1‑93 (MQ=255)
                                          tCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGAc                                               >  2:2355750/1‑93 (MQ=255)
                                            aaaGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACaa                                             <  1:3649249/93‑1 (MQ=255)
                                            aaaGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACaa                                             <  1:3354877/93‑1 (MQ=255)
                                             aaGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAt                                            >  2:3584919/1‑93 (MQ=255)
                                              aGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtg                                           <  2:1571833/93‑1 (MQ=255)
                                              aGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtg                                           <  2:2497177/93‑1 (MQ=255)
                                              aGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtg                                           <  2:230909/93‑1 (MQ=255)
                                               gATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtgt                                          >  2:2279660/1‑93 (MQ=255)
                                               aaTTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtgt                                          <  1:2589422/92‑1 (MQ=255)
                                                aTTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAAtgtg                                         >  1:1167139/1‑93 (MQ=255)
                                                 ttAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGc                                        <  1:1119692/93‑1 (MQ=255)
                                                   aaaaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGccc                                      <  1:2900749/93‑1 (MQ=255)
                                                   aaaaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGccc                                      <  2:3622626/93‑1 (MQ=255)
                                                    aaaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCg                                     >  2:1417370/1‑93 (MQ=255)
                                                    aaaaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCg                                     >  1:455575/1‑93 (MQ=255)
                                                      aaCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGAt                                   >  2:1203099/1‑93 (MQ=255)
                                                       aCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATg                                  >  2:4464778/1‑93 (MQ=255)
                                                       aCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATg                                  <  1:2279660/93‑1 (MQ=255)
                                                       aCTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATg                                  >  1:3235904/1‑93 (MQ=255)
                                                        cTCGTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGt                                 <  1:3804857/93‑1 (MQ=255)
                                                           gTTGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCgg                              >  1:197465/1‑93 (MQ=255)
                                                            ttGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGt                             <  2:3328995/93‑1 (MQ=255)
                                                             tGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTa                            <  1:2420758/93‑1 (MQ=255)
                                                             tGCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTa                            <  1:1062827/93‑1 (MQ=255)
                                                              gCAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTAt                           >  1:325271/1‑93 (MQ=255)
                                                               cAAATGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTAtt                          <  1:3741911/93‑1 (MQ=255)
                                                                aaaTGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTc                         >  2:4656725/1‑93 (MQ=255)
                                                                 aaTGCGCATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCg                        >  2:4235692/1‑93 (MQ=255)
                                                                 aaTGCGCATCGACGGTACAACTGGGATAGTAGACCAGGACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTct                        >  2:3399969/1‑92 (MQ=255)
                                                                      gcATCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTg                   <  1:860061/93‑1 (MQ=255)
                                                                        aTCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCTGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGGATTCGTATTGaa                 >  2:4209706/1‑93 (MQ=255)
                                                                         tCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGaaa                >  2:1366804/1‑93 (MQ=255)
                                                                         tCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGaaa                >  2:2993972/1‑93 (MQ=255)
                                                                         tCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGaaa                <  1:2042285/93‑1 (MQ=255)
                                                                         tCGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGaaa                >  1:4267051/1‑93 (MQ=255)
                                                                          cGACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGaaaa               <  1:303851/93‑1 (MQ=255)
                                                                           gACGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAt              <  2:1637441/93‑1 (MQ=255)
                                                                            aCGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTTCCCGATGTCGGTATTCGTATTGAAAAtc             <  1:970438/93‑1 (MQ=255)
                                                                             cGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAtct            <  2:3981563/93‑1 (MQ=255)
                                                                             cGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAtct            <  1:1103870/93‑1 (MQ=255)
                                                                             cGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAtct            >  2:1969282/1‑93 (MQ=255)
                                                                             cGGTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAtct            <  2:1480024/93‑1 (MQ=255)
                                                                              ggTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAAtctc           <  1:1650738/93‑1 (MQ=255)
                                                                               gTACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTGGGATTTGTATTGAAAATCTcgg         >  1:3381876/1‑93 (MQ=255)
                                                                               acaCAAATGGTGAAGTAGAGCAAGACAAACTGGTCGCCAGGTCGAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcg          <  1:3589847/91‑1 (MQ=255)
                                                                                tACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcgg         <  2:2076025/93‑1 (MQ=255)
                                                                                tACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcgg         <  1:3094814/93‑1 (MQ=255)
                                                                                tACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcgg         >  2:4386160/1‑93 (MQ=255)
                                                                                tACAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATATcgg         >  1:4190900/1‑93 (MQ=255)
                                                                                 aCAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggc        <  1:2430751/93‑1 (MQ=255)
                                                                                 aCAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggc        <  2:763497/93‑1 (MQ=255)
                                                                                 aCAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggc        >  1:2706804/1‑93 (MQ=255)
                                                                                 aCAACTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggc        <  1:408541/93‑1 (MQ=255)
                                                                                    aCTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggcggc     <  2:3211342/93‑1 (MQ=255)
                                                                                     cTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggcggcg    >  2:4021123/1‑93 (MQ=255)
                                                                                     cTGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggcggcg    >  1:1954628/1‑93 (MQ=255)
                                                                                      tGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggcggcgg   >  1:3734999/1‑93 (MQ=255)
                                                                                      tGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggcggcgg   <  1:3118107/93‑1 (MQ=255)
                                                                                      tGGGATAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTcggcggcgg   >  1:2486234/1‑93 (MQ=255)
                                                                                       gggaTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGa  <  2:4404586/93‑1 (MQ=255)
                                                                                          |                                                                                         
ACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAAAAACTCGTTGCAAATGCGCATCGACGGTACAACTGGGGTAGTAGACCAGTACAACCTGGTCGCCAGGCGAAGAAGTTCAGACAATGTGCCCGATGTCGGTATTCGTATTGAAAATCTCGGCGGCGGA  >  NC_000913/3289738‑3289917

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