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,2290AG99.8% 1286.6 / ‑5.2 403I400V (ATT→GTT) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (213/189);  total (214/189)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.34e-01

Only 100 of 406 total aligned reads displayed.

CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965139‑3965318
                                                                                          |                                                                                         
 tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTta                                                                                        >  1:4181587/1‑93 (MQ=255)
 tgctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTta                                                                                        >  2:1449925/1‑93 (MQ=255)
   ctCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaat                                                                                      <  2:2875329/93‑1 (MQ=255)
    tCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaata                                                                                     <  1:3416011/93‑1 (MQ=255)
     cACGACTCAGGAAGAACTGCAGAGAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaataa                                                                                    <  2:5651630/93‑1 (MQ=255)
     cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaataa                                                                                    <  2:1923753/93‑1 (MQ=255)
     cACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTtaataa                                                                                    <  2:7432144/93‑1 (MQ=255)
      aCGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATaaa                                                                                   <  2:3488388/93‑1 (MQ=255)
          cTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTgg                                                                               >  2:1546435/1‑93 (MQ=255)
           tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGc                                                                              >  1:1762753/1‑93 (MQ=255)
           tCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGc                                                                              >  2:2537649/1‑93 (MQ=255)
             aGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCaa                                                                            <  2:1352189/93‑1 (MQ=255)
             aGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCaa                                                                            >  1:6453395/1‑93 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           <  1:5035046/93‑1 (MQ=255)
              ggAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAAt                                                                           >  1:7437230/1‑93 (MQ=255)
               gaagaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATg                                                                          >  2:6879239/1‑93 (MQ=255)
                  gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgacc                                                                       >  1:3832662/1‑93 (MQ=255)
                  gaaCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgacc                                                                       >  1:4371033/1‑93 (MQ=255)
                    aCTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaa                                                                     <  1:2551332/93‑1 (MQ=255)
                     cTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaag                                                                    >  1:5915260/1‑93 (MQ=255)
                      tGCAGAAAATGTGGATCCTGCGCAAAATCATTCACGCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaaga                                                                   <  1:3121930/93‑1 (MQ=255)
                      tGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaaga                                                                   <  1:2493760/93‑1 (MQ=255)
                         aGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagacca                                                                <  2:3960041/93‑1 (MQ=255)
                          gAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATgaccaagaccaa                                                               >  1:7578734/1‑93 (MQ=255)
                            aaaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATCCCTCGTTAATAAACTGGCAATGACCAAGACCAATg                                                             <  2:4874426/93‑1 (MQ=255)
                             aaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGa                                                            >  2:4331673/1‑93 (MQ=255)
                             aaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGa                                                            <  1:3670448/93‑1 (MQ=255)
                             aaTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGa                                                            <  2:7304642/93‑1 (MQ=255)
                              aTGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGAc                                                           >  2:7501760/1‑93 (MQ=255)
                               tgtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACg                                                          >  2:7538957/1‑93 (MQ=255)
                                gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         <  2:7182246/93‑1 (MQ=255)
                                gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         <  2:3982410/93‑1 (MQ=255)
                                gtgGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGa                                                         <  1:683704/93‑1 (MQ=255)
                                   gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAttt                                                      <  1:4812093/93‑1 (MQ=255)
                                   gATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAttt                                                      >  2:6567234/1‑93 (MQ=255)
                                    aTCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGAtttt                                                     >  1:6402163/1‑92 (MQ=255)
                                      ccTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttctt                                                   >  1:5984747/1‑93 (MQ=255)
                                       cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                                                  <  2:3800623/93‑1 (MQ=255)
                                       cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                                                  <  2:3756808/93‑1 (MQ=255)
                                       cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                                                  >  1:1088889/1‑93 (MQ=255)
                                       cTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATttcttc                                                  <  1:1358854/93‑1 (MQ=255)
                                            cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             >  2:218587/1‑93 (MQ=255)
                                            cAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAat                                             >  2:6229544/1‑93 (MQ=255)
                                             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  1:2169823/1‑93 (MQ=255)
                                             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  2:1038838/1‑93 (MQ=255)
                                             aaaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatg                                            >  2:6354742/1‑93 (MQ=255)
                                              aaaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatga                                           <  2:2343096/93‑1 (MQ=255)
                                               aaTCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAatgat                                          >  1:6583556/1‑93 (MQ=255)
                                                 tCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAAtgatga                                        >  1:7893692/1‑93 (MQ=255)
                                                   aTTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGaaa                                      <  1:6967984/93‑1 (MQ=255)
                                                    ttCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAAc                                     <  1:4870706/93‑1 (MQ=255)
                                                     tCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACg                                    >  1:4546042/1‑93 (MQ=255)
                                                      cACCCGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGc                                   <  2:6733496/93‑1 (MQ=255)
                                                        cccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTc                                 <  1:2234763/93‑1 (MQ=255)
                                                         ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                >  2:3415037/1‑93 (MQ=255)
                                                         ccGATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCa                                <  2:2743808/93‑1 (MQ=255)
                                                           gATGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATa                              <  1:6327392/93‑1 (MQ=255)
                                                            aTGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaa                             >  1:1732056/1‑93 (MQ=255)
                                                            aTGGGCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATaa                             <  1:1000986/93‑1 (MQ=255)
                                                              gggCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAt                           <  1:7217794/93‑1 (MQ=255)
                                                               ggCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAtt                          >  2:5980988/1‑93 (MQ=255)
                                                               ggCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAtt                          >  2:7788843/1‑93 (MQ=255)
                                                               ggCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAtt                          <  1:950932/93‑1 (MQ=255)
                                                                gCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAttt                         >  2:3168707/1‑93 (MQ=255)
                                                                gCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAttt                         >  1:5775390/1‑93 (MQ=255)
                                                                gCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAttt                         <  2:5290341/93‑1 (MQ=255)
                                                                gCGAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAAttt                         <  2:4167551/93‑1 (MQ=255)
                                                                  gAAATCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGt                       <  1:1559047/93‑1 (MQ=255)
                                                                   aaaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTc                      <  1:3921925/93‑1 (MQ=255)
                                                                   aaaTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTc                      <  2:3339659/93‑1 (MQ=255)
                                                                     aTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCtt                    >  1:1484532/1‑93 (MQ=255)
                                                                     aTCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCtt                    >  1:4685735/1‑93 (MQ=255)
                                                                      tCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTa                   >  2:1807924/1‑93 (MQ=255)
                                                                      tCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTa                   >  1:2658971/1‑93 (MQ=255)
                                                                      tCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTa                   >  2:6614980/1‑93 (MQ=255)
                                                                      tCGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTa                   >  2:4976354/1‑93 (MQ=255)
                                                                       cGATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTAt                  <  1:7653737/93‑1 (MQ=255)
                                                                        gATGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATg                 >  1:6772187/1‑93 (MQ=255)
                                                                         aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                <  1:5264072/93‑1 (MQ=255)
                                                                         aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                <  1:1895106/93‑1 (MQ=255)
                                                                         aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                <  2:6941096/93‑1 (MQ=255)
                                                                         aTGCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGc                <  2:4548014/93‑1 (MQ=255)
                                                                           gCAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCa              <  2:3613438/93‑1 (MQ=255)
                                                                            cAATGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaa             <  1:2081833/93‑1 (MQ=255)
                                                                              aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           >  2:5383131/1‑93 (MQ=255)
                                                                              aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           >  2:5938170/1‑93 (MQ=255)
                                                                              aTGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaa           >  2:7692086/1‑93 (MQ=255)
                                                                               tGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaaa          >  1:5522992/1‑93 (MQ=255)
                                                                               tGGAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCaaaaa          >  1:3539635/1‑93 (MQ=255)
                                                                                ggAATTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAAc         >  2:1146817/1‑93 (MQ=255)
                                                                                   aTTCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGcc      >  1:4136141/1‑93 (MQ=255)
                                                                                    tgccTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCa     <  1:7828446/91‑1 (MQ=255)
                                                                                     tCCTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCAc    >  1:7740661/1‑93 (MQ=255)
                                                                                      ccTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAACCGCCACg   >  1:6117538/1‑93 (MQ=255)
                                                                                      ccTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACg   <  2:1281191/93‑1 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  <  2:2455860/93‑1 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  1:7285035/1‑93 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  1:2980629/1‑93 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  <  2:4776690/93‑1 (MQ=255)
                                                                                       cTCGTTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACgt  >  2:1674415/1‑93 (MQ=255)
                                                                                          |                                                                                         
CTGCTCACGACTCAGGAAGAACTGCAGAAAATGTGGATCCTGCGCAAAATCATTCACCCGATGGGCGAAATCGATGCAATGGAATTCCTCATTAATAAACTGGCAATGACCAAGACCAATGACGATTTCTTCGAAATGATGAAACGCTCATAAATTTGTCTTATGCCAAAAACGCCACGT  >  NC_000913/3965139‑3965318

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