Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,754,691 T→C V171A (GTT→GCT)  intA → prophage CP4‑57 integrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,754,6910TC99.7% 1193.0 / ‑5.1 371V171A (GTT→GCT) intAprophage CP4‑57 integrase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (197/172);  total (197/174)
Fisher's exact test for biased strand distribution p-value = 4.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

Only 100 of 380 total aligned reads displayed.

CAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAATATCGGAA  >  NC_000913/2754601‑2754780
                                                                                          |                                                                                         
cAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTc                                                                                         <  1:5949860/93‑1 (MQ=255)
  aTCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTcgc                                                                                       <  2:7733018/93‑1 (MQ=255)
  aTCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTcgc                                                                                       <  2:179719/93‑1 (MQ=255)
   tCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTcgcc                                                                                      >  1:5390329/1‑93 (MQ=255)
   tCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTcgcc                                                                                      <  1:4703358/93‑1 (MQ=255)
     ggTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTcggcgc                                                                                    >  1:2595786/1‑93 (MQ=255)
     ggTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTcgccgc                                                                                    >  2:4460225/1‑93 (MQ=255)
      gTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCgccgcc                                                                                   >  1:6782129/1‑93 (MQ=255)
      gTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCgccgcc                                                                                   >  2:6049934/1‑93 (MQ=255)
        gATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCtt                                                                                 <  1:1776761/93‑1 (MQ=255)
         atatCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCttt                                                                                <  1:7298284/93‑1 (MQ=255)
         atatCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCttt                                                                                <  2:4035423/93‑1 (MQ=255)
           atCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGt                                                                              <  2:6600442/93‑1 (MQ=255)
            tCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAGAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGATCGCCGCCTTTGTc                                                                             >  2:3465642/1‑93 (MQ=255)
            tCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTc                                                                             >  1:6139912/1‑93 (MQ=255)
            tCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTc                                                                             <  2:9219850/93‑1 (MQ=255)
            acAGTATCGCTGATATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTc                                                                             <  2:7498634/92‑1 (MQ=255)
             cAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCa                                                                            >  2:5514008/1‑93 (MQ=255)
              aGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAg                                                                           <  2:9580667/93‑1 (MQ=255)
               gTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGc                                                                          >  1:9274108/1‑93 (MQ=255)
               gTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGc                                                                          <  2:2068644/93‑1 (MQ=255)
                tATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCg                                                                         >  1:2884067/1‑93 (MQ=255)
                tATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGACAGAGGAGCATTAGAGACTGCTCGCCGCCTTTGTCAGCg                                                                         >  2:7216490/1‑93 (MQ=255)
                  tCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTa                                                                       >  2:1494117/1‑93 (MQ=255)
                     cTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTa                                                                    >  1:6436915/1‑93 (MQ=255)
                      tGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCTGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTaa                                                                   <  1:3712991/93‑1 (MQ=255)
                      tGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTaa                                                                   >  1:4930372/1‑93 (MQ=255)
                       gagaTTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAAc                                                                  <  1:5467003/93‑1 (MQ=255)
                       gagaTTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAAc                                                                  <  2:1020287/93‑1 (MQ=255)
                          aTTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGaa                                                               >  1:9807080/1‑93 (MQ=255)
                            tAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGt                                                             >  1:3275663/1‑93 (MQ=255)
                            tAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGt                                                             >  2:8388134/1‑93 (MQ=255)
                              aGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCa                                                           <  1:7602105/93‑1 (MQ=255)
                                  tCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAt                                                       <  1:9681343/93‑1 (MQ=255)
                                  tCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAt                                                       <  2:2821851/93‑1 (MQ=255)
                                  tCATACTCTGGTTAAAGCACTTCAGCCGGATCAGGCCACAGGCGCATTAGAGACTGCTCGCCGACTTTGTCAGCGTATTAATCAAGCCATgct                                                       >  2:4260525/1‑91 (MQ=255)
                                   cATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAtt                                                      >  1:5773049/1‑93 (MQ=255)
                                   cATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAtt                                                      <  2:5199227/93‑1 (MQ=255)
                                    aTACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAttt                                                     <  2:3862933/93‑1 (MQ=255)
                                    aTACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAttt                                                     >  2:3194170/1‑93 (MQ=255)
                                     tACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTa                                                    <  1:6251431/93‑1 (MQ=255)
                                     tACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTa                                                    >  1:5700774/1‑93 (MQ=255)
                                     tACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTa                                                    >  1:10009748/1‑93 (MQ=255)
                                     tACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTa                                                    <  2:9321893/93‑1 (MQ=255)
                                     tACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTa                                                    >  1:3766761/1‑93 (MQ=255)
                                     tACTCTGGTTAAAGCAGTTAAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTa                                                    <  1:9490348/93‑1 (MQ=255)
                                       ctctGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTCTg                                                  >  2:885799/1‑93 (MQ=255)
                                       ctctGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATg                                                  <  1:8247664/93‑1 (MQ=255)
                                        tctGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATgc                                                 <  1:8824530/93‑1 (MQ=255)
                                           ggTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCa                                              <  2:4611979/93‑1 (MQ=255)
                                            gTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAg                                             >  2:1102369/1‑93 (MQ=255)
                                              tAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGaa                                           <  1:3663503/93‑1 (MQ=255)
                                               aaaGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAac                                          >  1:8494370/1‑93 (MQ=255)
                                               aaaGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGAAGAac                                          >  1:2244288/1‑93 (MQ=255)
                                                aaGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAaca                                         <  2:2477452/93‑1 (MQ=255)
                                                 aGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAacac                                        <  1:5023101/93‑1 (MQ=255)
                                                 aGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAacac                                        >  1:5305437/1‑93 (MQ=255)
                                                  gCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAAcaca                                       <  1:5214870/93‑1 (MQ=255)
                                                     gTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGc                                    <  1:6906220/93‑1 (MQ=255)
                                                      ttCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGcc                                   <  1:8967793/93‑1 (MQ=255)
                                                       tCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCt                                  <  1:5627884/93‑1 (MQ=255)
                                                       tCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCt                                  <  2:2223308/93‑1 (MQ=255)
                                                       tCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCt                                  <  2:3348116/93‑1 (MQ=255)
                                                        cAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtg                                 >  1:3544806/1‑93 (MQ=255)
                                                          gCCGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgat                               <  1:7743826/93‑1 (MQ=255)
                                                           ccGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgatt                              <  1:4749055/93‑1 (MQ=255)
                                                           ccGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgatt                              >  1:9936147/1‑93 (MQ=255)
                                                            cGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAATGAAGTCATGATTTATGCGCAGAACACAGGCCtgattg                             <  2:8586977/93‑1 (MQ=255)
                                                            cGGTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattg                             >  2:5918133/1‑93 (MQ=255)
                                                            cGGTTCAGGCAAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattg                             <  1:2916375/93‑1 (MQ=255)
                                                             ggTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattga                            <  1:1173804/93‑1 (MQ=255)
                                                             ggTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattga                            <  2:3660880/93‑1 (MQ=255)
                                                             ggTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattga                            <  1:7915211/93‑1 (MQ=255)
                                                             ggTTCAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCtgattga                            <  2:1267819/93‑1 (MQ=255)
                                                                ttacgCCAGAGGCGCATTAAACACTGCTCGCCGCCTTTGTCAGCGTATTAACGACGTCATTATTTAGGCGCAGAACACAGGCCTGATTGATGc                         >  2:7009400/5‑93 (MQ=255)
                                                                 cAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCt                        <  1:2429469/93‑1 (MQ=255)
                                                                 cAGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCt                        <  1:907894/93‑1 (MQ=255)
                                                                  aGGCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTg                       >  2:6919622/1‑93 (MQ=255)
                                                                   ggCCAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGt                      <  1:3970353/93‑1 (MQ=255)
                                                                     ccAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTc                    >  1:2008157/1‑93 (MQ=255)
                                                                      cAGAGGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTcc                   >  2:2964235/1‑93 (MQ=255)
                                                                        gagGTGCATTAGAGACTGCTCGCCGCCTTTGTCGGCGTATTAACGAAGCCATGATTTATGCGCAGAAAACAGGCCTGATTGATGCTCTGCCtg                 >  2:560845/1‑92 (MQ=255)
                                                                        gagGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTa                 >  2:9762101/1‑93 (MQ=255)
                                                                        gagGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTa                 >  1:7121897/1‑93 (MQ=255)
                                                                         agGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAg                <  2:1639980/93‑1 (MQ=255)
                                                                         agGTGCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAg                >  1:315470/1‑93 (MQ=255)
                                                                             gCATTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTa            >  1:4253620/1‑93 (MQ=255)
                                                                             gCAATAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGAGGTATGCGCAGAACACAGGCCTGAGTGATGCTGGTCCGAGTGTa            >  2:8969359/1‑93 (MQ=255)
                                                                               aTTAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTaaa          <  2:9850646/93‑1 (MQ=255)
                                                                                ttAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAat         >  2:8216557/1‑93 (MQ=255)
                                                                                ttAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAat         >  2:4888039/1‑93 (MQ=255)
                                                                                ttAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAat         >  1:1559483/1‑93 (MQ=255)
                                                                                 tAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAata        <  2:5360168/93‑1 (MQ=255)
                                                                                 tAGAGACTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAata        <  2:5762326/93‑1 (MQ=255)
                                                                                   gagaCTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAATATc      <  2:194609/93‑1 (MQ=255)
                                                                                    agaCTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAATATCg     >  1:3025442/1‑93 (MQ=255)
                                                                                     gaCTGCTCGCAGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGAGCAGAACACAGGCCTGATTGATGCTGTTACTAGTGTAAATATCgg    >  2:8507749/1‑93 (MQ=255)
                                                                                      aCTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAATATCGGa   >  1:1665581/1‑93 (MQ=255)
                                                                                       cTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAATATCGGaa  >  2:648501/1‑93 (MQ=255)
                                                                                       cTGCTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAATATCGGaa  <  2:5829676/93‑1 (MQ=255)
                                                                                          |                                                                                         
CAATCGGTGATATCAGTATCACTGAGATTAAGGCTCATACTCTGGTTAAAGCAGTTCAGCCGGTTCAGGCCAGAGGTGCATTAGAGACTGTTCGCCGCCTTTGTCAGCGTATTAACGAAGTCATGATTTATGCGCAGAACACAGGCCTGATTGATGCTGTTCCTAGTGTAAATATCGGAA  >  NC_000913/2754601‑2754780

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