Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,566 G→A 17.5% pseudogene (128/129 nt) ykgQ → pseudogene, putative dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913317,5660GA17.5% 307.6 / 54.5 138pseudogene (128/129 nt)ykgQpseudogene, putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (51/62);  new base A (15/9);  total (66/72)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

Only 100 of 154 total aligned reads displayed.

CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  NC_000913/317467‑317665
                                                                                                   |                                                                                                     
cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattg                                                                                                       >  2:1194517/1‑100 (MQ=255)
 acaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattga                                                                                                      <  2:1268164/100‑1 (MQ=255)
 acaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattg                                                                                                       >  2:2676194/1‑99 (MQ=255)
   aaaTTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGata                                                                                                    <  2:3021718/100‑1 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                              >  1:1196258/1‑100 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                              <  1:2193025/100‑1 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                              >  2:2568965/1‑100 (MQ=255)
            cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                           <  2:1524690/100‑1 (MQ=255)
            cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                           <  2:3330283/100‑1 (MQ=255)
              cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                         >  1:307116/1‑100 (MQ=255)
              cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTa                                                                                         >  1:368876/1‑100 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑t                                                                                       >  2:1305789/1‑100 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTa                                                                                         <  1:1559417/99‑1 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑tt                                                                                      >  2:3032383/1‑99 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑ttt                                                                                     <  1:2034654/99‑1 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTa                                                                                   >  1:2865956/1‑100 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTa                                                                                   >  1:1125556/1‑100 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTTTTCGTTGTAAAAAAAAAAAAAAAATACATGATTGATATTAATCATGTA‑TTTTAc                                                                                  <  2:1628/100‑1 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAc                                                                                  >  1:236124/1‑100 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATACTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAc                                                                                  <  2:1454643/100‑1 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAc‑c                                                                                >  1:484591/1‑100 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAc‑c                                                                                >  2:2002038/1‑100 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAc‑c                                                                                <  1:511287/100‑1 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAc‑c                                                                                <  2:1629704/100‑1 (MQ=255)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAGATACATGCTTGATATTAATCATGTA‑TTTTAc‑c                                                                                <  1:3458145/99‑1 (MQ=255)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTA‑tttttaca                                                                                <  1:1137808/100‑5 (MQ=38)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTa                                                                              <  2:1438549/100‑1 (MQ=255)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTa                                                                              <  1:1855105/100‑1 (MQ=255)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTA‑tttttacag                                                                               <  1:1370906/100‑6 (MQ=37)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAAGTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATCAATCATGGA‑tttttaca                                                                                <  2:1972507/99‑5 (MQ=255)
                        gTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAAGACATGATTAATATCAATCATGTA‑tttttacagc                                                                              <  2:1554393/100‑7 (MQ=37)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAg                                                                            <  1:435486/100‑1 (MQ=255)
                           aTGTTGGGGAGTCGGTAAAACATTTTAGTCCTGAATTTGTTAGTGGTTGTCTTGCTGTAACAATAAATGATTGATATTATTCATGTG‑TTTTAA‑CTaa                                                                             >  2:1195559/1‑97 (MQ=255)
                           aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAgac                                                                          <  2:953793/100‑1 (MQ=255)
                            tGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCCGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAgacg                                                                         <  1:2101687/100‑1 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAgacga                                                                        >  2:1988165/1‑100 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAgacga                                                                        >  2:316503/1‑100 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTCTAATCCTGAATTTGTTTGTGGGTGCCTTGATGTAAAAACACACGACTGATAATCATCAGGTA‑TTTTGC‑CTAAgacga                                                                        >  2:1918129/1‑100 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATTTA‑TTTTAC‑CTAAGACGACa                                                                      >  2:697112/1‑100 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACa                                                                      >  1:1978538/1‑100 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACaa                                                                     <  2:2431511/100‑1 (MQ=255)
                                       cGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCtt                                                              <  1:1011010/100‑1 (MQ=255)
                                        ggcaaaaCATTTTAATCCTGACTTTGTTCGTCGTTGTCTTGCTGTAAAAATACATGATTGATTTTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAAcc                                                                >  2:823123/4‑97 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCtttt                                                            <  1:160949/100‑1 (MQ=255)
                                          tctaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTa                                                           <  2:1245369/97‑1 (MQ=255)
                                            aaaCATTTAAATCCTGAATTTTTTCGTGTTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAAc                                                         <  2:3269830/100‑1 (MQ=255)
                                             aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACa                                                        >  1:3232480/1‑100 (MQ=255)
                                              aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAg                                                       <  2:2075055/100‑1 (MQ=255)
                                                aTTTTAAACCTAAAATTGTTCGTGGTTGTCTTGCTATAAAAAAAAATGATTGATATTAATCATGTA‑TTTTAA‑AAAAAAAAAAAAAAACCTTTTAACAggg                                                     <  2:900279/100‑1 (MQ=255)
                                                 ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAg                                                       >  2:1513848/1‑97 (MQ=255)
                                                 ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGa                                                    <  1:3615160/100‑1 (MQ=255)
                                                  tttAATCCTGAATTTGTTCGTGGTTGGCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGGCGACAAAAACATTTTAACAGGGAg                                                   <  2:284183/100‑1 (MQ=255)
                                                     aaTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGa                                                <  2:697203/100‑1 (MQ=255)
                                                       tCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAATTACATGATTGGTATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAGACCTTTTAACCGGGCGTgc                                                >  2:1718029/1‑97 (MQ=255)
                                                        ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                             >  2:2137996/1‑100 (MQ=255)
                                                        ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                             >  2:359710/1‑100 (MQ=255)
                                                         cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                             >  2:3307828/1‑99 (MQ=255)
                                                         cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                            <  1:2878034/100‑1 (MQ=255)
                                                                tGTTCGTGGTTGTCTTGCTGTAAAAATACATGTTTGATATTACTCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTc                                     <  1:1279143/100‑1 (MQ=255)
                                                                 gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCa                                    >  2:2315582/1‑100 (MQ=255)
                                                                   tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                  >  1:1395869/1‑100 (MQ=255)
                                                                       ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                              <  2:1734560/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  1:1275583/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  1:1008134/100‑1 (MQ=255)
                                                                         cttccgatcTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  2:436575/92‑1 (MQ=255)
                                                                          tGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑ATAAAAAAAAAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                           <  2:1308657/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          <  1:748500/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          >  2:920263/1‑100 (MQ=255)
                                                                               tGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtctc                      >  2:2038423/1‑100 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                     <  1:2517272/100‑1 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                     >  2:1826724/1‑100 (MQ=255)
                                                                                 taggtaaaaTACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                    >  1:1398563/6‑100 (MQ=37)
                                                                                 taggtaaaaTACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                    >  1:2306129/6‑100 (MQ=37)
                                                                                 cTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                    <  1:328842/100‑1 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   <  2:3438409/100‑1 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   <  2:2745689/100‑1 (MQ=255)
                                                                                   ggtaaaaTACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  >  1:411598/4‑100 (MQ=38)
                                                                                   ggtaaaaTACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  >  1:2758502/4‑100 (MQ=38)
                                                                                   ggtaaaaTACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  >  1:603339/4‑100 (MQ=38)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCaagatcggaag                  >  1:436577/1‑90 (MQ=255)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  <  1:3386872/100‑1 (MQ=255)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAAGAAGTGATATTGATCTTCACTCTGTCATATAACCGGt                  >  2:3476835/1‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:704964/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:1356980/3‑100 (MQ=255)
                                                                                     aaaaaTACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTaa                >  1:2764106/1‑100 (MQ=255)
                                                                                       aaaTACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              <  1:1609834/100‑1 (MQ=255)
                                                                                       aaaTACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              <  1:2429826/100‑1 (MQ=255)
                                                                                          tACATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgn           <  1:1239494/100‑2 (MQ=255)
                                                                                           aCATGATTGATACTAATCATGTA‑TTTTAC‑CTAAGACTACAAAACCCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGc          <  1:2970844/100‑1 (MQ=255)
                                                                                            cATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg         >  1:3564135/1‑100 (MQ=255)
                                                                                            cATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCACCAGTaccagcggg         >  2:803532/1‑92 (MQ=255)
                                                                                            cATGATTAATATCAATCATGTATTTTTACAGCAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg           >  1:1870696/1‑100 (MQ=21)
                                                                                            cATGATTAATATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg         >  1:3203307/1‑100 (MQ=255)
                                                                                              tgattgatATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc       <  1:1808904/100‑1 (MQ=255)
                                                                                              tgattgatATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc       <  1:1816846/100‑1 (MQ=255)
                                                                                              tgattgatATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc       <  1:2270790/100‑1 (MQ=255)
                                                                                              tgattaatATCAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg         >  2:498955/1‑98 (MQ=255)
                                                                                                  tgatATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGc   >  1:3506227/1‑100 (MQ=255)
                                                                                                   gatATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCt  >  2:2088937/1‑100 (MQ=255)
                                                                                                   gatATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCt  <  1:2924208/100‑1 (MQ=255)
                                                                                                   |                                                                                                     
CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAC‑CTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  NC_000913/317467‑317665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑