Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,566 G→A 16.0% 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,5660GA16.0% 345.8 / 45.8 149pseudogene (128/129 nt)ykgQpseudogene, putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (57/65);  new base A (12/12);  total (70/79)
Fisher's exact test for biased strand distribution p-value = 8.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

Only 100 of 163 total aligned reads displayed.

CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  NC_000913/317467‑317665
                                                                                                   |                                                                                                   
cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattg                                                                                                     <  1:340433/100‑1 (MQ=255)
cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattg                                                                                                     >  1:1455159/1‑100 (MQ=255)
cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattg                                                                                                     <  1:553149/100‑1 (MQ=255)
 acaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattga                                                                                                    >  1:111417/1‑100 (MQ=255)
   aaaTTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGata                                                                                                  <  2:1122226/100‑1 (MQ=255)
   aaaTTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGata                                                                                                  <  1:1618641/100‑1 (MQ=255)
     aTTAGTTCGCTTGATTTTTGTTATGTTGGGTAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATAtt                                                                                                <  1:2438077/100‑1 (MQ=255)
       tAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATACTGAATTTGTTCGTGGTTGGCTTGCTGTAAAAATACATGATTAatatc                                                                                                <  1:3048414/98‑2 (MQ=255)
        aGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAAt                                                                                             <  1:2821620/100‑1 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                            >  1:2342291/1‑100 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                            <  1:1592085/100‑1 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                            <  1:1991378/100‑1 (MQ=255)
          ttCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                           <  1:2954979/100‑1 (MQ=255)
           tCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATGAATCa                                                                                           <  2:2137698/99‑1 (MQ=255)
            cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATg                                                                                         <  1:644650/100‑1 (MQ=255)
             gtttGATTTGGGTTATGTTGGTGAGTCGGTAAAACATTTTAATCCTGAATGTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGt                                                                                        <  1:3396026/98‑1 (MQ=255)
             gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTACTCATg                                                                                         >  2:2379617/1‑99 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAt                                                                                      >  1:1363319/1‑100 (MQ=255)
                tGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                     <  2:3843326/100‑1 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                     >  1:2331541/1‑99 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtt                                                                                     >  2:3882163/1‑99 (MQ=255)
                  aTTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAtttt                                                                                   <  2:3755727/100‑1 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAATATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTa                                                                                  <  1:1666013/100‑1 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTTGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAttttt                                                                                  <  1:1024244/100‑2 (MQ=38)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAttttt                                                                                  <  2:344764/100‑2 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATTCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTTATCATGTATTTTAcc                                                                                <  2:3776854/100‑1 (MQ=255)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttaca                                                                               <  1:2156155/100‑5 (MQ=38)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTTTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacag                                                                              <  1:113710/100‑6 (MQ=37)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                              <  2:3275914/100‑1 (MQ=255)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATAAATGATTAATATTAATCATGTATTTTACCTa                                                                              >  2:370435/1‑100 (MQ=255)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAATACATGATTAATATCAATCATGTAtttttacagc                                                                             <  1:1278788/100‑7 (MQ=25)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTCAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                              <  1:2748008/100‑1 (MQ=255)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAg                                                                            <  1:2879104/100‑1 (MQ=255)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAg                                                                            <  2:3662592/100‑1 (MQ=255)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacagc                                                                             <  2:1886689/99‑7 (MQ=37)
                           aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                          <  1:2605847/100‑1 (MQ=255)
                           aTGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgac                                                                          <  2:689694/100‑1 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACa                                                                      >  1:3740523/1‑100 (MQ=255)
                                  ggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaaaa                                                                   >  1:2447385/1‑100 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaaaaa                                                                  >  2:272784/1‑99 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                 >  2:558569/1‑100 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                <  2:3494285/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                >  2:759955/1‑100 (MQ=255)
                                       cGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTATCTAAGACGACAAAAACCtt                                                              >  2:1792333/1‑100 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                            <  1:3516688/100‑1 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                            <  1:300082/100‑1 (MQ=255)
                                          tAAAACCTTTTAATCCTGAATTGGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACATAAAACAACAAAAACCttt                                                             <  2:1085074/98‑1 (MQ=255)
                                           aaaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTaa                                                          <  2:1566616/100‑1 (MQ=255)
                                            aaaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                         <  2:352071/100‑1 (MQ=255)
                                                 ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGa                                                    >  1:1761900/1‑100 (MQ=255)
                                                    tAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTg                                                 <  2:81786/100‑1 (MQ=255)
                                                    tAATCCTGAATCTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTg                                                 >  1:2674858/1‑100 (MQ=255)
                                                         cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                            >  1:3742789/1‑100 (MQ=255)
                                                         cTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                            >  1:3564306/1‑100 (MQ=255)
                                                          tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAAACGACAAAAACATTTTAACAGGGAATGATAtt                                            >  1:3462137/1‑99 (MQ=255)
                                                           gAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGa                                          >  2:1354441/1‑100 (MQ=255)
                                                           gAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGa                                          <  1:2497688/100‑1 (MQ=255)
                                                           gAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGa                                          >  1:917274/1‑100 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTTATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                        <  2:2784112/100‑1 (MQ=255)
                                                             aTTTGTTAGTGGTTGTCTTGATGTAAAAATATATGATTGATATTAAGCATGTATTTTAGATAAGACGAAAAAAACCTTTTAACAGGGAGTGATATTGATc                                        <  2:202023/100‑1 (MQ=255)
                                                              tttGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCt                                       >  2:1258566/1‑100 (MQ=255)
                                                                tGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTc                                     <  1:1496154/100‑1 (MQ=255)
                                                                   tCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCAct                                  <  2:2849858/100‑1 (MQ=255)
                                                                    cGTTGTTGTCTTGCTGTAAAAATACTTGATTTATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                 <  1:3117538/100‑1 (MQ=255)
                                                                    cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                 <  1:1089322/100‑1 (MQ=255)
                                                                      tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTg                               <  2:111413/100‑1 (MQ=255)
                                                                      tGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTg                               >  2:3121305/1‑100 (MQ=255)
                                                                       ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                              >  2:1914978/1‑100 (MQ=255)
                                                                        gTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTc                             <  2:2439202/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  2:3183411/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  2:654318/100‑1 (MQ=255)
                                                                          tGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                           >  2:1737623/1‑100 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          <  1:3651327/100‑1 (MQ=255)
                                                                            tCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCatat                         <  2:2674834/100‑1 (MQ=255)
                                                                            tCTTGCTGTAAAAATACATGATTGATATTAATAATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGAAGTGATATTGATCTTCACTCTGTCatat                         >  1:3725612/1‑100 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                     <  1:3155169/100‑1 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   <  1:558576/100‑1 (MQ=255)
                                                                                  aggtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   >  2:134436/5‑100 (MQ=38)
                                                                                   gTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAACTTTTAACAGGAAGTGATATTGATCTTCACTCTGTCATATCTCCGGt                  >  1:465665/1‑100 (MQ=255)
                                                                                    tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 <  1:2264433/100‑1 (MQ=255)
                                                                                    tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:358054/1‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  2:3442439/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  2:3198768/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:3930265/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCAGGTa                 >  1:399049/3‑100 (MQ=38)
                                                                                      aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat               >  1:1079213/1‑100 (MQ=255)
                                                                                      aaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAat               <  1:2191131/100‑1 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  2:766833/1‑100 (MQ=255)
                                                                                        aaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat             >  2:733583/1‑100 (MQ=255)
                                                                                        aaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  1:3784725/1‑99 (MQ=255)
                                                                                          tACATGTTTGTTTTTAATCATGTATTTTATTTATGATTAATATTACTTTTTAAAAGGGAGTGGTATGGATCTTCACTCTGTCATATCTCCGGTAATATgg           <  1:1231052/100‑1 (MQ=255)
                                                                                            cATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCg         <  2:1664042/100‑1 (MQ=255)
                                                                                             aTGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGt        <  2:3482957/100‑1 (MQ=255)
                                                                                             aTGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGt        <  2:2017810/100‑1 (MQ=255)
                                                                                              tgattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTc       <  1:1835615/100‑1 (MQ=255)
                                                                                               gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa      >  1:3598182/1‑100 (MQ=255)
                                                                                               gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa      <  1:2972001/100‑1 (MQ=255)
                                                                                               gattgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCa      <  1:355465/100‑1 (MQ=255)
                                                                                                  tgatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGc   >  2:2571529/1‑100 (MQ=255)
                                                                                                   gatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGCAGTGATATTGATCTTCACTCTGTAATATCTTCGGTAATATGGCGTCAGACt  <  2:3626746/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 < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: