Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,965,229 A→G I400V (ATT→GTT)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,965,2290AG100.0% 882.0 / NA 269I400V (ATT→GTT) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (153/116);  total (153/116)

Only 100 of 270 total aligned reads displayed.

CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965139‑3965318
                                                                                          |                                                                                         
ctgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGtt                                                                                         >  2:6262341/1‑93 (MQ=255)
 tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCAGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTta                                                                                        >  1:4553997/1‑93 (MQ=255)
    tCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaata                                                                                     <  2:1878483/93‑1 (MQ=255)
    tCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaata                                                                                     <  1:539737/93‑1 (MQ=255)
    tCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaata                                                                                     <  1:3355697/93‑1 (MQ=255)
       cGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAAc                                                                                  >  2:2403562/1‑93 (MQ=255)
        gACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACt                                                                                 >  1:4477907/1‑93 (MQ=255)
        gACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACt                                                                                 <  2:5552141/93‑1 (MQ=255)
          cTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTgg                                                                               >  1:3092566/1‑93 (MQ=255)
           tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGc                                                                              >  1:887286/1‑93 (MQ=255)
            cAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCa                                                                             >  2:4205916/1‑93 (MQ=255)
             aGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCaa                                                                            >  2:3871245/1‑93 (MQ=255)
             aGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCaa                                                                            >  1:2160247/1‑93 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGTTGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           >  2:5075542/1‑93 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           <  1:1013071/93‑1 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           >  2:3584806/1‑93 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           >  1:1945838/1‑93 (MQ=255)
                 agaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgac                                                                        >  1:2780118/1‑93 (MQ=255)
                 agaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgac                                                                        >  2:4385330/1‑93 (MQ=255)
                  gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgacc                                                                       <  2:2792910/93‑1 (MQ=255)
                  gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgacc                                                                       >  2:2927641/1‑93 (MQ=255)
                    aCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaa                                                                     <  1:2697800/93‑1 (MQ=255)
                      tGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaaga                                                                   >  2:1813577/1‑93 (MQ=255)
                          gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagaccaa                                                               >  1:1851973/1‑93 (MQ=255)
                           aaaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAAt                                                              <  2:5104428/93‑1 (MQ=255)
                           aaaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAAt                                                              >  2:1747422/1‑93 (MQ=255)
                           aaaaTGTGGATCCTGCGCAAAATAATTCACCCGATGGGCGAAATCGATGCAATGAAATTCCTCGTTAATAAACTGGCAATgacaaagaccaag                                                              >  1:4193156/1‑92 (MQ=255)
                            aaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATg                                                             >  2:1471041/1‑93 (MQ=255)
                            aaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATg                                                             <  1:3546246/93‑1 (MQ=255)
                            aaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATg                                                             <  2:1276892/93‑1 (MQ=255)
                             aaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGa                                                            <  1:3779779/93‑1 (MQ=255)
                             aaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGa                                                            >  1:2407531/1‑93 (MQ=255)
                              aTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGAc                                                           <  2:4356650/93‑1 (MQ=255)
                                gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         <  1:6295467/93‑1 (MQ=255)
                                gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         <  2:1646331/93‑1 (MQ=255)
                                 tgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAt                                                        >  1:6106601/1‑93 (MQ=255)
                                  ggATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATTGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAtt                                                       >  2:1830765/1‑93 (MQ=255)
                                  ggATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAtt                                                       <  2:351283/93‑1 (MQ=255)
                                      ccTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttctt                                                   >  2:249042/1‑93 (MQ=255)
                                      ccTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttctt                                                   >  2:1364169/1‑93 (MQ=255)
                                        tGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCg                                                 <  2:5826480/93‑1 (MQ=255)
                                         gcgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGa                                                >  1:5592106/1‑93 (MQ=255)
                                         gcgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGa                                                >  2:4563066/1‑93 (MQ=255)
                                          cgcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGaa                                               <  1:2335839/93‑1 (MQ=255)
                                           gcAAAATCATTCGCCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGaaa                                              <  1:4004421/93‑1 (MQ=255)
                                           gcAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGaca                                              <  2:6209138/93‑3 (MQ=255)
                                            cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             <  1:505848/93‑1 (MQ=255)
                                             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            <  2:4070594/93‑1 (MQ=255)
                                             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  1:3161757/1‑93 (MQ=255)
                                             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  2:3763841/1‑93 (MQ=255)
                                              aaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                                           <  2:2233818/93‑1 (MQ=255)
                                              aaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                                           >  1:4232306/1‑93 (MQ=255)
                                               aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCATTGACGATTTCTTCGAAatgat                                          <  2:4007121/93‑1 (MQ=255)
                                                  cATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaa                                       >  2:5711724/1‑93 (MQ=255)
                                                  cATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaa                                       <  1:2366807/93‑1 (MQ=255)
                                                  cATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaa                                       <  2:5157321/93‑1 (MQ=255)
                                                  cATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaa                                       >  1:5546499/1‑93 (MQ=255)
                                                   aTTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa                                      <  1:2241769/93‑1 (MQ=255)
                                                   aTTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa                                      <  1:5243968/93‑1 (MQ=255)
                                                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAAc                                     >  2:2859076/1‑93 (MQ=255)
                                                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAAc                                     <  1:5178735/93‑1 (MQ=255)
                                                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAAc                                     >  2:6163351/1‑93 (MQ=255)
                                                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAAc                                     <  1:1691713/93‑1 (MQ=255)
                                                       aCCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCt                                  >  2:4826258/1‑93 (MQ=255)
                                                       aCCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCt                                  >  1:2979912/1‑93 (MQ=255)
                                                        cccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTc                                 >  1:1331304/1‑93 (MQ=255)
                                                        cccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTc                                 >  2:3213383/1‑93 (MQ=255)
                                                         ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                >  2:2697200/1‑93 (MQ=255)
                                                          cGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt                               >  2:4442113/1‑93 (MQ=255)
                                                          cGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt                               >  2:1132923/1‑93 (MQ=255)
                                                          cGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCAt                               >  2:5620201/1‑93 (MQ=255)
                                                           gATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATa                              >  2:3130866/1‑93 (MQ=255)
                                                           gATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATa                              >  2:3255723/1‑93 (MQ=255)
                                                            aTGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaa                             >  1:1165419/1‑93 (MQ=255)
                                                             tGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa                            <  2:5394991/93‑1 (MQ=255)
                                                             tGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaaa                            >  1:983971/1‑93 (MQ=255)
                                                                gCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAttt                         <  1:1121903/93‑1 (MQ=255)
                                                                   aaaTCGATGCAATGGGATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTc                      >  2:5990770/1‑93 (MQ=255)
                                                                   aaaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTc                      <  1:718016/93‑1 (MQ=255)
                                                                    aaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCt                     >  2:4479355/1‑93 (MQ=255)
                                                                     aTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCtt                    <  2:693526/93‑1 (MQ=255)
                                                                       cGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTAt                  >  1:1415071/1‑93 (MQ=255)
                                                                        gATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATg                 <  1:4646512/93‑1 (MQ=255)
                                                                         aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                >  2:1952125/1‑93 (MQ=255)
                                                                         aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                >  1:482307/1‑93 (MQ=255)
                                                                          tGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGcc               >  1:3676666/1‑93 (MQ=255)
                                                                           gCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCa              <  1:5438722/93‑1 (MQ=255)
                                                                            cAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaa             >  1:2799498/1‑93 (MQ=255)
                                                                             aaTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaa            <  1:1369599/93‑1 (MQ=255)
                                                                              aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           >  2:3009757/1‑93 (MQ=255)
                                                                              aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           <  2:5854971/93‑1 (MQ=255)
                                                                               tGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaaa          >  2:4610322/1‑93 (MQ=255)
                                                                                ggAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAc         >  1:5779901/1‑93 (MQ=255)
                                                                                ggAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAc         >  1:3541468/1‑93 (MQ=255)
                                                                                 gAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACg        >  2:2009217/1‑93 (MQ=255)
                                                                                   aTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGcc      >  1:4713243/1‑93 (MQ=255)
                                                                                   aTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGcc      <  1:3260884/93‑1 (MQ=255)
                                                                                    ttCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCa     <  2:3038431/93‑1 (MQ=255)
                                                                                     tCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAc    <  2:3439541/93‑1 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  1:1519134/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965139‑3965318

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