Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,915,980 G→A G502G (GGC→GGT atpA ← membrane‑bound ATP synthase, F1 sector, alpha‑subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,915,9800GA84.8% 836.1 / 104.0 336G502G (GGC→GGTatpAmembrane‑bound ATP synthase, F1 sector, alpha‑subunit
Reads supporting (aligned to +/- strand):  ref base G (29/22);  new base A (152/132);  total (181/155)
Fisher's exact test for biased strand distribution p-value = 7.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.35e-02

Only 100 of 336 total aligned reads displayed.

CCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGC  >  NC_000913/3915890‑3916069
                                                                                          |                                                                                         
ccATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATAcc                                                                                         <  1:4282218/93‑1 (MQ=255)
  aTGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCtt                                                                                       >  2:6874784/1‑93 (MQ=255)
  aTGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCtt                                                                                       <  1:7428571/93‑1 (MQ=255)
   tGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCttt                                                                                      >  2:129663/1‑93 (MQ=255)
    gAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTc                                                                                     >  1:2116649/1‑93 (MQ=255)
    gAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTc                                                                                     <  1:6922130/93‑1 (MQ=255)
    gAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTc                                                                                     <  2:6174892/93‑1 (MQ=255)
     aGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCa                                                                                    <  2:749097/93‑1 (MQ=255)
     aGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCa                                                                                    >  1:1945288/1‑93 (MQ=255)
       cTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGc                                                                                  <  2:7031864/93‑1 (MQ=255)
       cTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGc                                                                                  <  2:5438725/93‑1 (MQ=255)
        ttCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTTAGCt                                                                                 >  1:6287718/1‑93 (MQ=255)
        ttCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCt                                                                                 <  1:6713543/93‑1 (MQ=255)
         tctcCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCtt                                                                                >  2:4095271/1‑93 (MQ=255)
         tctcCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCtt                                                                                <  2:4396512/93‑1 (MQ=255)
            cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGcc                                                                             >  1:346610/1‑93 (MQ=255)
            cctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGcc                                                                             <  2:2905669/93‑1 (MQ=255)
             ctcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCt                                                                            >  1:3876723/1‑93 (MQ=255)
              tcAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCtt                                                                           <  2:3620406/93‑1 (MQ=255)
               cAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTc                                                                          <  1:708236/93‑1 (MQ=255)
                aaTGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCg                                                                         <  1:6200834/93‑1 (MQ=255)
                   gCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGAtt                                                                      <  1:5105324/93‑1 (MQ=255)
                    ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGAttt                                                                     <  2:4991003/93‑1 (MQ=255)
                    ccTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGAttt                                                                     >  2:3318630/1‑93 (MQ=255)
                     cTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTtc                                                                    <  2:2730690/93‑1 (MQ=255)
                     cTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTtc                                                                    <  1:376951/93‑1 (MQ=255)
                       tGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTtcgt                                                                  >  1:272227/1‑93 (MQ=255)
                         cGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTtcgtcg                                                                <  2:1916292/93‑1 (MQ=255)
                          ggCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTcgtcgt                                                               >  2:2140331/1‑93 (MQ=255)
                             cTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGt                                                            <  1:4625501/93‑1 (MQ=255)
                              tGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTa                                                           <  2:260822/93‑1 (MQ=255)
                                cccTAAGGAAAGCCCCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATCCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGc                                                         <  1:2485533/93‑1 (MQ=255)
                                 ccTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGcc                                                        >  1:4410563/1‑93 (MQ=255)
                                  cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCa                                                       >  1:4517464/1‑93 (MQ=255)
                                  cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCa                                                       >  1:6401576/1‑93 (MQ=255)
                                  cTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCa                                                       <  2:6306950/93‑1 (MQ=255)
                                   tAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAc                                                      >  1:1382335/1‑93 (MQ=255)
                                   tAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAc                                                      >  1:1460919/1‑93 (MQ=255)
                                   tAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAc                                                      >  2:4601460/1‑93 (MQ=255)
                                   tAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAc                                                      >  2:2946683/1‑93 (MQ=255)
                                    aaGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCAcc                                                     <  1:5129185/93‑1 (MQ=255)
                                      ggCAAGCCTCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCgg                                                   >  2:3527943/1‑93 (MQ=255)
                                       gCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGt                                                  <  2:6496795/93‑1 (MQ=255)
                                        cAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTc                                                 <  1:3160726/93‑1 (MQ=255)
                                         aaGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCt                                                <  1:5043872/93‑1 (MQ=255)
                                          aGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTg                                               >  1:752901/1‑93 (MQ=255)
                                          aGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTg                                               >  2:1777338/1‑93 (MQ=255)
                                               ccAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGa                                          >  1:2597761/1‑93 (MQ=255)
                                                cAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAt                                         <  2:3768271/93‑1 (MQ=255)
                                                cAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAt                                         <  1:2186821/93‑1 (MQ=255)
                                                 aGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtc                                        >  2:1499157/1‑93 (MQ=255)
                                                 aGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtc                                        <  2:3235727/93‑1 (MQ=255)
                                                   aCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtctc                                      >  1:3981446/1‑93 (MQ=255)
                                                   aCGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGAtctc                                      >  2:7512986/1‑93 (MQ=255)
                                                    cGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATctct                                     <  1:1150059/93‑1 (MQ=255)
                                                    cGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATctct                                     >  1:362238/1‑93 (MQ=255)
                                                    cGTTACCAGGATTGGCTTGCTTTGAAGGAATCGCGGATGCCTTTCAGCTTGCCTTCGATCTCGTCGGTCTAGCCCCCGGTCCGCGTGAtctca                                     >  1:705483/1‑92 (MQ=255)
                                                       tACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGc                                  <  1:959479/93‑1 (MQ=255)
                                                         ccAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCAt                                >  2:3008724/1‑93 (MQ=255)
                                                         ccAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCAt                                >  1:7222253/1‑93 (MQ=255)
                                                          cAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATc                               >  2:4027871/1‑93 (MQ=255)
                                                          cAGGATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATc                               <  1:1561523/93‑1 (MQ=255)
                                                            ggATTGGGTTGCTTTGAAGGAATTGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                             >  1:5707990/1‑93 (MQ=255)
                                                            ggATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                             <  2:3117267/93‑1 (MQ=255)
                                                            ggATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                             <  1:5650095/93‑1 (MQ=255)
                                                            ggATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                             >  1:5891644/1‑93 (MQ=255)
                                                            ggATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                             >  1:4909216/1‑93 (MQ=255)
                                                            ggATTGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCaa                             >  1:5410530/1‑93 (MQ=255)
                                                              aTTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACg                           >  2:1949670/1‑93 (MQ=255)
                                                               ttGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACgg                          >  1:6762448/1‑93 (MQ=255)
                                                               ttGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACgg                          <  1:2442209/93‑1 (MQ=255)
                                                               ttGGGTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACgg                          >  1:3856714/1‑93 (MQ=255)
                                                                   gTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCg                      <  2:7178144/93‑1 (MQ=255)
                                                                   gTTGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCg                      <  2:5213247/93‑1 (MQ=255)
                                                                    ttGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGt                     <  2:5264944/93‑1 (MQ=255)
                                                                    ttGCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGt                     <  1:1038241/93‑1 (MQ=255)
                                                                      gCTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGa                   >  2:2502946/1‑93 (MQ=255)
                                                                       cTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGAt                  <  1:3079638/93‑1 (MQ=255)
                                                                       cTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGAt                  <  1:652248/93‑1 (MQ=255)
                                                                       cTTTGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGAt                  >  2:5565857/1‑93 (MQ=255)
                                                                        tttGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                 <  2:381774/93‑1 (MQ=255)
                                                                        tttGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                 <  1:6067543/93‑1 (MQ=255)
                                                                        tttGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATc                 >  2:5107475/1‑93 (MQ=255)
                                                                         ttGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCa                >  2:6675465/1‑93 (MQ=255)
                                                                          tGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCAc               <  1:3552587/93‑1 (MQ=255)
                                                                          tGAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCAc               <  2:1627013/93‑1 (MQ=255)
                                                                           gAAGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACg              >  2:2231040/1‑93 (MQ=255)
                                                                            aaGGAATCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACgg             >  2:6616642/1‑93 (MQ=255)
                                                                                aaTCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGa         <  1:1219742/93‑1 (MQ=255)
                                                                                aaTCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGa         >  2:5643954/1‑93 (MQ=255)
                                                                                 aTCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGAc        <  2:5921273/93‑1 (MQ=255)
                                                                                  tCGAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACg       >  2:1432875/1‑93 (MQ=255)
                                                                                    gAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTa     <  1:3016301/93‑1 (MQ=255)
                                                                                    gAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTa     >  2:4803535/1‑93 (MQ=255)
                                                                                    gAGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTa     >  2:2274746/1‑93 (MQ=255)
                                                                                     aGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTaa    >  1:29017/1‑93 (MQ=255)
                                                                                     aGGATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTaa    <  2:1134929/93‑1 (MQ=255)
                                                                                      ggATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAg   >  1:5318288/1‑93 (MQ=255)
                                                                                       gATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGc  >  2:4490085/1‑93 (MQ=255)
                                                                                       gATACCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTCCATCAACGGAGCGTGATCACGGTCGACGTAAGc  <  2:3910221/93‑1 (MQ=255)
                                                                                          |                                                                                         
CCATGAGCTTCTCCTCAATGCCTTGCGGCCTGCCCTAAGGCAAGCCGCCAGACGTTACCAGGATTGGGTTGCTTTGAAGGAATCGAGGATGCCTTTCAGCTTGCCTTCGATTTCGTCGTTGTAGCCACCGGTCTGGTTGATCTCTTGCATCAACGGAGCGTGATCACGGTCGACGTAAGC  >  NC_000913/3915890‑3916069

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