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.

CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  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)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAt                                                                                       >  2:1305789/1‑100 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTa                                                                                         <  1:1559417/99‑1 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                      >  2:3032383/1‑99 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                     <  1:2034654/99‑1 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                   >  1:2865956/1‑100 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                   >  1:1125556/1‑100 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTTTTCGTTGTAAAAAAAAAAAAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                  <  2:1628/100‑1 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                  >  1:236124/1‑100 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATACTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                  <  2:1454643/100‑1 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAcc                                                                                >  1:484591/1‑100 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAcc                                                                                >  2:2002038/1‑100 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAcc                                                                                <  1:511287/100‑1 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAcc                                                                                <  2:1629704/100‑1 (MQ=255)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAGATACATGCTTGATATTAATCATGTATTTTAcc                                                                                <  1:3458145/99‑1 (MQ=255)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttaca                                                                                <  1:1137808/100‑5 (MQ=38)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                              <  2:1438549/100‑1 (MQ=255)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                              <  1:1855105/100‑1 (MQ=255)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacag                                                                               <  1:1370906/100‑6 (MQ=37)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAAGTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATCAATCATGGAtttttaca                                                                                <  2:1972507/99‑5 (MQ=255)
                        gTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAAGACATGATTAATATCAATCATGTAtttttacagc                                                                              <  2:1554393/100‑7 (MQ=37)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAg                                                                            <  1:435486/100‑1 (MQ=255)
                           aTGTTGGGGAGTCGGTAAAACATTTTAGTCCTGAATTTGTTAGTGGTTGTCTTGCTGTAACAATAAATGATTGATATTATTCATGTGTTTTAACTaa                                                                             >  2:1195559/1‑97 (MQ=255)
                           aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                          <  2:953793/100‑1 (MQ=255)
                            tGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCCGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacg                                                                         <  1:2101687/100‑1 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacga                                                                        >  2:1988165/1‑100 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacga                                                                        >  2:316503/1‑100 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTCTAATCCTGAATTTGTTTGTGGGTGCCTTGATGTAAAAACACACGACTGATAATCATCAGGTATTTTGCCTAAgacga                                                                        >  2:1918129/1‑100 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATTTATTTTACCTAAGACGACa                                                                      >  2:697112/1‑100 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACa                                                                      >  1:1978538/1‑100 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                     <  2:2431511/100‑1 (MQ=255)
                                       cGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtt                                                              <  1:1011010/100‑1 (MQ=255)
                                        ggcaaaaCATTTTAATCCTGACTTTGTTCGTCGTTGTCTTGCTGTAAAAATACATGATTGATTTTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                >  2:823123/4‑97 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                            <  1:160949/100‑1 (MQ=255)
                                          tctaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTa                                                           <  2:1245369/97‑1 (MQ=255)
                                            aaaCATTTAAATCCTGAATTTTTTCGTGTTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                         <  2:3269830/100‑1 (MQ=255)
                                             aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                        >  1:3232480/1‑100 (MQ=255)
                                              aCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAg                                                       <  2:2075055/100‑1 (MQ=255)
                                                aTTTTAAACCTAAAATTGTTCGTGGTTGTCTTGCTATAAAAAAAAATGATTGATATTAATCATGTATTTTAAAAAAAAAAAAAAAAACCTTTTAACAggg                                                     <  2:900279/100‑1 (MQ=255)
                                                 ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAg                                                       >  2:1513848/1‑97 (MQ=255)
                                                 ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGa                                                    <  1:3615160/100‑1 (MQ=255)
                                                  tttAATCCTGAATTTGTTCGTGGTTGGCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGGCGACAAAAACATTTTAACAGGGAg                                                   <  2:284183/100‑1 (MQ=255)
                                                     aaTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGa                                                <  2:697203/100‑1 (MQ=255)
                                                       tCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAATTACATGATTGGTATTAATCATGTATTTTACCTAAGACGACAAAGACCTTTTAACCGGGCGTgc                                                >  2:1718029/1‑97 (MQ=255)
                                                        ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                             >  2:2137996/1‑100 (MQ=255)
                                                        ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                             >  2:359710/1‑100 (MQ=255)
                                                         cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGatat                                             >  2:3307828/1‑99 (MQ=255)
                                                         cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                            <  1:2878034/100‑1 (MQ=255)
                                                                tGTTCGTGGTTGTCTTGCTGTAAAAATACATGTTTGATATTACTCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTc                                     <  1:1279143/100‑1 (MQ=255)
                                                                 gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCa                                    >  2:2315582/1‑100 (MQ=255)
                                                                   tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                  >  1:1395869/1‑100 (MQ=255)
                                                                       ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                              <  2:1734560/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  1:1275583/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  1:1008134/100‑1 (MQ=255)
                                                                         cttccgatcTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  2:436575/92‑1 (MQ=255)
                                                                          tGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACATAAAAAAAAAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                           <  2:1308657/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          <  1:748500/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          >  2:920263/1‑100 (MQ=255)
                                                                               tGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtctc                      >  2:2038423/1‑100 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                     <  1:2517272/100‑1 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                     >  2:1826724/1‑100 (MQ=255)
                                                                                 taggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                    >  1:1398563/6‑100 (MQ=37)
                                                                                 taggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                    >  1:2306129/6‑100 (MQ=37)
                                                                                 cTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCg                    <  1:328842/100‑1 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   <  2:3438409/100‑1 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   <  2:2745689/100‑1 (MQ=255)
                                                                                   ggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  >  1:411598/4‑100 (MQ=38)
                                                                                   ggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  >  1:2758502/4‑100 (MQ=38)
                                                                                   ggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  >  1:603339/4‑100 (MQ=38)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCaagatcggaag                  >  1:436577/1‑90 (MQ=255)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  <  1:3386872/100‑1 (MQ=255)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAAGAAGTGATATTGATCTTCACTCTGTCATATAACCGGt                  >  2:3476835/1‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:704964/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:1356980/3‑100 (MQ=255)
                                                                                     aaaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTaa                >  1:2764106/1‑100 (MQ=255)
                                                                                       aaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              <  1:1609834/100‑1 (MQ=255)
                                                                                       aaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              <  1:2429826/100‑1 (MQ=255)
                                                                                          tACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgn           <  1:1239494/100‑2 (MQ=255)
                                                                                           aCATGATTGATACTAATCATGTATTTTACCTAAGACTACAAAACCCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGc          <  1:2970844/100‑1 (MQ=255)
                                                                                            cATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg         >  1:3564135/1‑100 (MQ=255)
                                                                                            cATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCACCAGTaccagcggg         >  2:803532/1‑92 (MQ=255)
                                                                                            cATGATTAATATCAATCATGTATTTTTACAGCAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATgg           >  1:1870696/1‑100 (MQ=21)
                                                                                            cATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg         >  1:3203307/1‑100 (MQ=255)
                                                                                              tgattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc       <  1:1808904/100‑1 (MQ=255)
                                                                                              tgattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc       <  1:1816846/100‑1 (MQ=255)
                                                                                              tgattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc       <  1:2270790/100‑1 (MQ=255)
                                                                                              tgattaatATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg         >  2:498955/1‑98 (MQ=255)
                                                                                                  tgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGc   >  1:3506227/1‑100 (MQ=255)
                                                                                                   gatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCt  >  2:2088937/1‑100 (MQ=255)
                                                                                                   gatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCt  <  1:2924208/100‑1 (MQ=255)
                                                                                                   |                                                                                                     
CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  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: