Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,566 G→A 18.3% 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,5660GA18.3% 408.7 / 71.4 180pseudogene (128/129 nt)ykgQpseudogene, putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (57/89);  new base A (16/17);  total (73/107)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 198 total aligned reads displayed.

CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGC  >  NC_000913/317467‑317664
                                                                                                   |                                                                                                  
cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattg                                                                                                    >  1:2532118/1‑100 (MQ=255)
  caAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattgat                                                                                                  <  1:110511/100‑1 (MQ=255)
   tccgatcttCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGata                                                                                                 <  2:2026738/93‑1 (MQ=255)
   aaaTTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAACTACATGATTGata                                                                                                 <  2:2584008/100‑1 (MQ=255)
   aaaTTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGata                                                                                                 <  2:3283845/100‑1 (MQ=255)
    aaTTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGCATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATCGatat                                                                                                <  2:620060/100‑1 (MQ=255)
      ttAGTTCGCTTGATTTGAGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTCATGTAAAAATAGAtgattgat                                                                                                  >  2:1590333/1‑96 (MQ=255)
         ttccgatctaTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                           <  2:2136979/91‑1 (MQ=255)
           tCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCAt                                                                                         >  1:113065/1‑100 (MQ=255)
           tCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAAATTGTTCGGGGTTGTCTTGCTGTAAAAATACATGATTGataagatcgga                                                                                         >  1:2026741/1‑92 (MQ=255)
            cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                        <  1:408307/100‑1 (MQ=255)
            cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                        >  1:3747954/1‑100 (MQ=255)
             gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGt                                                                                       >  2:2489788/1‑100 (MQ=255)
              cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                      >  2:2423189/1‑100 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCTGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAt                                                                                     <  2:2111/100‑1 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAt                                                                                     >  1:1785761/1‑100 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                   >  1:1726272/1‑100 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                  <  2:4129057/100‑1 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                  >  1:4037421/1‑100 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                  <  2:3091671/100‑1 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                  <  1:2673684/100‑1 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttt                                                                                  <  2:3861919/100‑1 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                 >  1:26242/1‑100 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                 <  1:1964063/100‑1 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                 <  1:1397161/100‑1 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                <  1:2143411/100‑1 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                <  1:2293326/100‑1 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                >  1:2135720/1‑100 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                <  1:296177/100‑1 (MQ=255)
                    ttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAc                                                                                <  1:2803307/100‑1 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttac                                                                               <  2:4057700/100‑4 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATACTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttac                                                                               <  2:402407/100‑4 (MQ=38)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAAGCCGGAATTTGTCCGTGGTTGTCTTGCTCTAAACATACATGCTGAATATCAGTCATGTAtttttaaa                                                                              <  2:3376537/100‑5 (MQ=255)
                        gTTATGTTGGGTAGTCTGTAAAACATTGTAATCCTGAATTTTTTCTTTGTTGTCTTGCTGTAAAAATACATGGTTAATATCAATCATGTAtttttacag                                                                             <  2:1855261/99‑6 (MQ=21)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacagca                                                                           <  2:1677842/100‑8 (MQ=37)
                          tATGTTGGTGAGTCGGTAAAATATTTTAATACTGAATTTATTCGTGGTTGTCTTGCTATAAAAATACATTATTAATATCAATCATGTAtttttacag                                                                             <  2:188320/97‑6 (MQ=21)
                           aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                         >  2:29671/1‑100 (MQ=255)
                              ttGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacgac                                                                      <  2:440147/100‑1 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACa                                                                     <  2:2911786/100‑1 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                <  2:73264/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                               <  1:3388649/100‑1 (MQ=255)
                                        ggTAAAAGATTTTAAGCCTGAATTTGCTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATCAATCATGTATTTTACCTAATACGACAAAAACCttt                                                            <  2:1516118/100‑1 (MQ=38)
                                        ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                            >  1:3722764/1‑100 (MQ=255)
                                        ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                            >  1:746466/1‑100 (MQ=255)
                                        ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                            <  1:2423194/100‑1 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                           <  1:853299/100‑1 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                            <  1:3167248/99‑1 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAAAAAAAAAAAAAAAAACCttt                                                            <  2:116044/99‑1 (MQ=255)
                                             aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                       >  1:1123766/1‑100 (MQ=255)
                                                aTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAggg                                                    >  1:211716/1‑100 (MQ=255)
                                                aTTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAAAAAAAAAAAAAAAAACCTTTTAACAggg                                                    <  2:2840420/100‑1 (MQ=255)
                                                  tttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAg                                                  <  1:3007495/100‑1 (MQ=255)
                                                  tttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAg                                                  <  1:2493307/100‑1 (MQ=255)
                                                   ttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                 <  2:2683360/100‑1 (MQ=255)
                                                   ttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                 <  1:2930403/100‑1 (MQ=255)
                                                    tAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                 >  2:3440221/1‑99 (MQ=255)
                                                        ccTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATTATGTATTTTACCTAAGACGACAAGAACCTTTTAATAGGGTTTGatat                                            >  2:3704546/1‑100 (MQ=255)
                                                         cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                           >  1:276633/1‑100 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                       >  1:2388442/1‑100 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                       <  2:1997976/100‑1 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                       <  2:2216148/100‑1 (MQ=255)
                                                                 gTTCGTGGTTGTCTTGCTGTAAAAATAGATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAAAGGGAGTGATATTGCTATTCa                                   >  2:1251341/1‑100 (MQ=255)
                                                                  ttCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAc                                  <  1:1511752/100‑1 (MQ=255)
                                                                  ttCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAc                                  >  2:3572147/1‑100 (MQ=255)
                                                                   tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                 >  1:3366600/1‑100 (MQ=255)
                                                                   tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                 >  1:2888450/1‑100 (MQ=255)
                                                                    cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                >  2:864069/1‑100 (MQ=255)
                                                                    cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                <  1:2661402/100‑1 (MQ=255)
                                                                      tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTg                              >  1:1021110/1‑100 (MQ=255)
                                                                       ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                             <  1:3259433/100‑1 (MQ=255)
                                                                       ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                             >  1:999596/1‑100 (MQ=255)
                                                                        gTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                            <  2:922859/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                           <  1:3628869/100‑1 (MQ=255)
                                                                          tGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                          <  1:728182/100‑1 (MQ=255)
                                                                          tGTATTGCTTTAAAAATACATGATTGATATTAATTATGTATTTTAAATAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                             <  2:2881703/97‑1 (MQ=255)
                                                                             gtcttaggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtc                       >  1:278049/10‑100 (MQ=25)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                    <  2:4098326/100‑1 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                    <  2:2325397/100‑1 (MQ=255)
                                                                                 taggtaaaaTACATGATTAATAACAATCATGTGTTTTACCTAAGACGACAAAAAACTTTTAACAGGGATTGATATCGATCATCACTCTATCATAtgtc                     >  2:4184897/6‑98 (MQ=21)
                                                                                 cTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTAACTCTGTCATAtctac                    >  2:3915990/1‑97 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                  >  2:983931/1‑100 (MQ=255)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                 >  2:768535/1‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                 >  1:1266988/3‑99 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                >  1:4169300/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCAAAGACGACAAAAACCTTTTAACAGGGTGTGATATTGATCTTCACTCTGTCATATCTCCg                   >  2:2248830/3‑97 (MQ=37)
                                                                                    ggtaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCGTTTAACAGGGAGTGATATTTATCTTCACTCTGTCATATCTCCGGTa                >  2:1336161/4‑100 (MQ=38)
                                                                                      aaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat              >  1:1298576/1‑100 (MQ=255)
                                                                                      aaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat              >  1:1322453/1‑100 (MQ=255)
                                                                                      aaaaTACATGATAGATATTAATCATGAATTTAACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat              <  1:1184031/100‑1 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata             >  2:3602631/1‑100 (MQ=255)
                                                                                        aaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat            >  2:1771267/1‑100 (MQ=255)
                                                                                         aTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATG           >  1:3897535‑M1/1‑100 (MQ=255)
                                                                                           aCATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGC         <  2:2698548‑M1/100‑1 (MQ=255)
                                                                                            cATGATTTATACTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCG        <  2:531451‑M1/100‑1 (MQ=255)
                                                                                              tgattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTC      >  2:282336‑M1/1‑100 (MQ=255)
                                                                                               gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGTGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGT       >  2:2954454‑M1/1‑98 (MQ=255)
                                                                                               gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCA     <  1:3767434‑M1/100‑1 (MQ=255)
                                                                                               gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCA     >  1:3196154‑M1/1‑100 (MQ=255)
                                                                                               gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCA     <  2:1785757‑M1/100‑1 (MQ=255)
                                                                                                attgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAG    >  2:655556‑M1/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGC  >  NC_000913/317467‑317664

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: