Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,028,708 A→G V78A (GTG→GCG)  vsr ← DNA mismatch endonuclease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,028,7080AG99.6% 808.9 / ‑5.0 252V78A (GTG→GCG) vsrDNA mismatch endonuclease
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (122/129);  total (122/130)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 257 total aligned reads displayed.

CCGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGT  >  NC_000913/2028618‑2028797
                                                                                          |                                                                                         
ccGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCt                                                                                         <  2:6802536/93‑1 (MQ=255)
 cGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCtt                                                                                        >  1:7812780/1‑93 (MQ=255)
  gAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCttt                                                                                       >  2:9682147/1‑93 (MQ=255)
  gAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCttt                                                                                       <  2:4694831/93‑1 (MQ=255)
   aGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTa                                                                                      >  1:6552674/1‑93 (MQ=255)
    gTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAGCATTTTTCCCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTaa                                                                                     <  1:7065692/93‑1 (MQ=255)
     ttCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTaaa                                                                                    >  2:3802007/1‑93 (MQ=255)
     ttCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTaaa                                                                                    <  2:7917075/93‑1 (MQ=255)
       ccTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACa                                                                                  >  2:19126/1‑93 (MQ=255)
       ccTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGACAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACa                                                                                  >  1:5338620/1‑93 (MQ=255)
        cTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAg                                                                                 >  2:9574517/1‑93 (MQ=255)
         tGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGa                                                                                <  1:3147155/93‑1 (MQ=255)
           cAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATa                                                                              >  1:1662988/1‑93 (MQ=255)
             aGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGc                                                                            <  1:5900430/93‑1 (MQ=255)
                gACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTCCTCCAGCCAGACCTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGt                                                                         <  1:9667534/93‑1 (MQ=255)
                 aCTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtg                                                                        >  1:5599463/1‑93 (MQ=255)
                  cTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtga                                                                       <  1:2869747/93‑1 (MQ=255)
                  cTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtga                                                                       <  2:8010672/93‑1 (MQ=255)
                  cTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtga                                                                       <  1:1575046/93‑1 (MQ=255)
                      tatCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtgatgat                                                                   >  1:4685638/1‑93 (MQ=255)
                       atCGCGGCGCGCGCTCTCAACATTTTTACCTATCTTCGCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtgatgatg                                                                  <  2:2602251/93‑1 (MQ=255)
                        tCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtgatgatga                                                                 >  2:925986/1‑93 (MQ=255)
                        tCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGtgatgatga                                                                 <  1:9728806/93‑1 (MQ=255)
                          gcggcgATCGCGCTCACCATTTTTACCTATCTTCTCCGGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGatgatgatg                                                               <  2:5228864/93‑1 (MQ=255)
                          gcggcgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGatgatgatg                                                               <  2:4982199/93‑1 (MQ=255)
                          gcggcgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGatgatgatg                                                               >  1:6180415/1‑93 (MQ=255)
                           cggcgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGc                                                              <  2:9200938/93‑1 (MQ=255)
                           cggcgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGc                                                              >  1:54411/1‑93 (MQ=255)
                             gcgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCa                                                            >  1:2555666/1‑93 (MQ=255)
                              cgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAg                                                           >  2:7078692/1‑93 (MQ=255)
                              cgATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAg                                                           <  2:2463158/93‑1 (MQ=255)
                               gATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGa                                                          >  2:8620534/1‑93 (MQ=255)
                                gtCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGaa                                                         <  1:6306629/92‑1 (MQ=255)
                                aTCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGaa                                                         >  1:1625206/1‑93 (MQ=255)
                                aTCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGaa                                                         >  1:9870721/1‑93 (MQ=255)
                                aTCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGaa                                                         >  2:2501072/1‑93 (MQ=255)
                                 tCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAag                                                        <  2:6755777/93‑1 (MQ=255)
                                  cgcgCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAagc                                                       <  2:8956068/93‑1 (MQ=255)
                                   gcgcTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAagca                                                      <  2:755488/93‑1 (MQ=255)
                                    cgcTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAagcag                                                     >  2:8412933/1‑93 (MQ=255)
                                     gcTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAagcagc                                                    <  1:8479425/93‑1 (MQ=255)
                                      cTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGcc                                                   >  1:2707573/1‑93 (MQ=255)
                                       tCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCa                                                  <  2:7634709/93‑1 (MQ=255)
                                        cAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCAt                                                 >  2:3050393/1‑93 (MQ=255)
                                        cAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCAt                                                 >  1:9080161/1‑93 (MQ=255)
                                        cAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATCCCAGAAGCAGCCAt                                                 >  1:4062045/1‑93 (MQ=255)
                                         aaCATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATg                                                <  1:8641448/93‑1 (MQ=255)
                                          aCATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATgg                                               <  2:9364935/93‑1 (MQ=255)
                                           cATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGGGATGATGATGCCAGAAGCAGCCATggg                                              >  1:7251669/1‑93 (MQ=255)
                                             tttttACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTa                                            >  2:8314869/1‑93 (MQ=255)
                                                ttACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAat                                         <  1:1422920/93‑1 (MQ=255)
                                                ttACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAat                                         <  2:1782256/93‑1 (MQ=255)
                                                ttACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAat                                         <  2:2441231/93‑1 (MQ=255)
                                                ttACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAat                                         >  2:4337036/1‑93 (MQ=255)
                                                 tACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAata                                        <  1:7026268/93‑1 (MQ=255)
                                                 tACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAata                                        >  1:5044372/1‑93 (MQ=255)
                                                  tccTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAatat                                       <  1:2602982/92‑1 (MQ=255)
                                                  aCCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAatat                                       <  1:4154583/93‑1 (MQ=255)
                                                   ccTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATc                                      <  2:5750556/93‑1 (MQ=255)
                                                    cTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCa                                     <  2:8180217/93‑1 (MQ=255)
                                                    cTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCa                                     >  2:3401367/1‑93 (MQ=255)
                                                     tATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCAc                                    >  2:2632198/1‑93 (MQ=255)
                                                     tATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCAc                                    >  1:4826520/1‑93 (MQ=255)
                                                     tATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGACAGAAGCAGCCATGGGTAAATATCAc                                    >  1:1768066/1‑93 (MQ=255)
                                                      aTCTTCTCCAGCCAGAACTCGTTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACg                                   <  1:4572962/93‑1 (MQ=255)
                                                      aTCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACg                                   >  2:7352873/1‑93 (MQ=255)
                                                      aTCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACg                                   >  1:4482543/1‑93 (MQ=255)
                                                       tcttctCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGc                                  >  1:7150947/1‑93 (MQ=255)
                                                       tcttctCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGc                                  >  1:9127696/1‑93 (MQ=255)
                                                       tcttctCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGc                                  >  1:3935479/1‑93 (MQ=255)
                                                        cttctCCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCa                                 >  1:10008612/1‑93 (MQ=255)
                                                          tctcCAGCCAGAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGc                               <  1:7304786/93‑1 (MQ=255)
                                                               agccagAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGatat                          >  2:2984525/1‑93 (MQ=255)
                                                               agccagAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGatat                          >  1:3059915/1‑93 (MQ=255)
                                                                gccagAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATAtt                         <  1:3620090/93‑1 (MQ=255)
                                                                 ccagAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtc                        <  2:2062455/93‑1 (MQ=255)
                                                                 ccagAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtc                        <  1:925986/93‑1 (MQ=255)
                                                                  cagAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtca                       <  2:1306678/93‑1 (MQ=255)
                                                                  cagAACTCGGTTCGAGCCGCAGGCGCTTAAAACAGATAGCAGTGATGATGATGCAAGAAGCAGCAATGGGTAAATATCACGCAGCGATATtca                       <  2:5436284/93‑1 (MQ=255)
                                                                   agAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtcat                      >  2:314045/1‑93 (MQ=255)
                                                                   agAACTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtcat                      >  1:2076732/1‑93 (MQ=255)
                                                                     aaCTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGGTGATGCCAGAAGCAGACATGGGTAAATATCACGCAGCGATATtcatga                    >  2:9938884/1‑91 (MQ=255)
                                                                     aaCTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtcatca                    >  1:1120253/1‑93 (MQ=255)
                                                                     aaCTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtcatca                    >  2:2528113/1‑93 (MQ=255)
                                                                     aaCTCGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATtcatca                    >  1:774044/1‑93 (MQ=255)
                                                                         cGGTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGa                <  1:197720/93‑1 (MQ=255)
                                                                          ggTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGAc               >  1:2355461/1‑93 (MQ=255)
                                                                           gTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACa              <  2:8895821/93‑1 (MQ=255)
                                                                           gTTCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACa              <  2:4457358/93‑1 (MQ=255)
                                                                            ttCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACaa             <  1:1531378/93‑1 (MQ=255)
                                                                             tCGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACaaa            <  1:67464/93‑1 (MQ=255)
                                                                              cGAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACaaaa           <  2:6088827/93‑1 (MQ=255)
                                                                               gAGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAAt          >  1:1851609/1‑93 (MQ=255)
                                                                                aGTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATc         >  1:124977/1‑93 (MQ=255)
                                                                                 gTCGCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATcc        >  1:2421494/1‑93 (MQ=255)
                                                                                    gCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGa     <  1:424548/93‑1 (MQ=255)
                                                                                    gCAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGa     >  2:4363977/1‑93 (MQ=255)
                                                                                     cAGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGAc    <  2:5925460/93‑1 (MQ=255)
                                                                                      aGGCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACg   >  2:9630961/1‑93 (MQ=255)
                                                                                       ggCGCTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGTGATATTCATCAACGACAAAATCCGGACGt  >  1:5618176/1‑93 (MQ=255)
                                                                                          |                                                                                         
CCGAGTTCCTGCAAGCGACTGATATCGCGGCGATCGCGCTCAACATTTTTACCTATCTTCTCCAGCCAGAACTCGGTTCGAGTCGCAGGCACTTTAAACAGATAGCAGTGATGATGATGCCAGAAGCAGCCATGGGTAAATATCACGCAGCGATATTCATCAACGACAAAATCCGGACGT  >  NC_000913/2028618‑2028797

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