Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 317,566 G→A 18.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,5660GA18.0% 246.7 / 42.3 111pseudogene (128/129 nt)ykgQpseudogene, putative dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (39/51);  new base A (14/6);  total (54/57)
Fisher's exact test for biased strand distribution p-value = 4.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

Only 100 of 114 total aligned reads displayed.

CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  NC_000913/317467‑317665
                                                                                                   |                                                                                                    
cacaAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACAtgattg                                                                                                      >  1:1481463/1‑100 (MQ=255)
       tAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                               <  2:1083220/100‑1 (MQ=255)
       tAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTGCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTaa                                                                                               >  1:962223/1‑100 (MQ=255)
        aGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCTTGAATTTGTTCGTGGTTGTGTTGCTGTAAAAATACATGATTGATATTAAt                                                                                              >  1:799775/1‑100 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                             >  2:346971/1‑100 (MQ=255)
         gTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATc                                                                                             <  1:1901363/100‑1 (MQ=255)
          ttCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCa                                                                                            <  2:741591/100‑1 (MQ=255)
          ttCGATTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTATAAAAAAAAAAAAATACATGATTGATATTAATc                                                                                             <  1:2035717/99‑1 (MQ=255)
           tCGATTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCAt                                                                                           >  1:1402068/1‑100 (MQ=255)
            cGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTTGTGGTTGTCTTGCTGTAACAATACATGATTGATATTAATCATg                                                                                          >  1:865155/1‑100 (MQ=255)
             gCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGt                                                                                         >  1:1930554/1‑100 (MQ=255)
              cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                        >  2:1860980/1‑100 (MQ=255)
              cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                        <  2:2078588/100‑1 (MQ=255)
              cTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                        <  2:847314/100‑1 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAt                                                                                      <  2:2248043/100‑1 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTCCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                        >  1:2422074/1‑99 (MQ=255)
                tGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtt                                                                                     <  2:124587/100‑1 (MQ=255)
                tGATTTGGATTATGTTGGGCAGTCGGTAAAACATTTTAATCCTGAATATGTTCGTGGTTATCTTGCTGTGGAAATACATAAGTAATAGCAATCATGTa                                                                                        <  1:2425220/98‑1 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                    <  2:1481472/100‑1 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTAttt                                                                                    >  2:56214/1‑100 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAttttt                                                                                  <  2:700835/100‑2 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAcc                                                                                <  2:761801/100‑1 (MQ=255)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttaca                                                                               <  2:1536490/100‑5 (MQ=38)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttaca                                                                               <  2:1885464/100‑5 (MQ=38)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTa                                                                              <  2:1733747/100‑1 (MQ=255)
                        gTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTaa                                                                             >  2:1257592/1‑100 (MQ=255)
                        gTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTAtttttacagc                                                                             <  1:1081871/100‑7 (MQ=37)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAg                                                                            <  2:1633467/100‑1 (MQ=255)
                            tGTTGGGGGGTCGGTAAAACATTTTAATCCTGAATTTGTAAGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTGATCATGTATTTTACCTAAg                                                                            <  1:780019/97‑1 (MQ=255)
                            tGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacg                                                                         >  2:1924808/1‑100 (MQ=255)
                            tGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacg                                                                         <  1:943267/100‑1 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAgacga                                                                        >  1:2201066/1‑100 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACa                                                                      <  2:1085894/100‑1 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                     >  1:2283468/1‑100 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                     >  1:1784866/1‑100 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACaa                                                                     <  2:2026323/100‑1 (MQ=255)
                                  ggatacGGTAAAACATTTTAATCCTGACTTTTTTCGTGGTTGTCTTGCTGTTAGTATTCCTGATTTATATTTATCACGTATTTTGCTTAAGACGACaata                                                                   >  1:1975049/6‑98 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                 >  1:288848/1‑100 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                 <  2:2110318/100‑1 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAc                                                                 <  2:386734/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATTCTGAATTTGTTTGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                <  1:1007148/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                >  2:1098377/1‑100 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                <  2:233541/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                <  2:1185792/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                >  1:1309606/1‑100 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                >  2:1463317/1‑100 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAAcc                                                                <  1:1292454/100‑1 (MQ=255)
                                       cGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtt                                                              >  1:885906/1‑100 (MQ=255)
                                        ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                             <  1:2011361/100‑1 (MQ=255)
                                        ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCttt                                                             >  1:1852569/1‑100 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                            <  2:54261/100‑1 (MQ=255)
                                          tAAAACATTTTAATCCTGAATTTGTGCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCtttt                                                            <  1:1836704/99‑1 (MQ=255)
                                             aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAAc                                                         <  1:1616362/99‑1 (MQ=255)
                                             aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACa                                                        <  1:2396586/100‑1 (MQ=255)
                                               cATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTAAAAAAAAAAAAAAAAACCTTTTAACAgg                                                      <  1:404842/100‑1 (MQ=255)
                                                 ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGa                                                    >  1:1742851/1‑100 (MQ=255)
                                                   ttAATCCTCAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                  <  1:1258533/100‑1 (MQ=255)
                                                     aaTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGa                                                <  2:1964511/100‑1 (MQ=255)
                                                      aTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGat                                               <  1:1259201/100‑1 (MQ=255)
                                                          tGAATTTGTTCGTTGTTGTCTTGATGTAAAAATACATGATTGATATTAATCATGTATGTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                            <  1:2097074/99‑1 (MQ=255)
                                                          tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTg                                           <  2:1159581/100‑1 (MQ=255)
                                                           gAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGa                                          >  2:2309809/1‑100 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                        <  2:1295462/100‑1 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                        <  2:358259/100‑1 (MQ=255)
                                                               ttGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCtt                                      <  2:1480839/100‑1 (MQ=255)
                                                                 gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCa                                    >  1:36500/1‑100 (MQ=255)
                                                                    cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                 >  2:1701740/1‑100 (MQ=255)
                                                                    cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                 >  1:782566/1‑100 (MQ=255)
                                                                     gTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                >  1:14890/1‑100 (MQ=255)
                                                                       ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                              <  2:219636/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  1:124656/100‑1 (MQ=255)
                                                                          tGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                           <  2:444869/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          <  1:965763/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          <  2:1379335/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          >  1:1485711/1‑100 (MQ=255)
                                                                             gtcttaggtaaaaTACATGATTAATATCAATCATGTATTTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGAGATATTGATCTTCAATCTGTCata                          >  1:47979/10‑99 (MQ=18)
                                                                              ttGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtct                       >  1:2103696/1‑100 (MQ=255)
                                                                              ttGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtct                       >  1:987164/1‑100 (MQ=255)
                                                                               tGCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtct                       >  1:2321372/1‑99 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                     >  1:2175967/1‑100 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   <  2:885909/100‑1 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   >  2:1051729/1‑100 (MQ=255)
                                                                                    tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 <  1:779920/100‑1 (MQ=255)
                                                                                    tAAAAATACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 <  2:2379186/100‑1 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTGCgg                   >  1:946274/3‑98 (MQ=38)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:668519/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:613542/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  2:40306/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  2:654378/3‑100 (MQ=255)
                                                                                     aaaaaTACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTaa                >  1:820656/1‑100 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  1:1644183/1‑100 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  2:1077973/1‑100 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  2:1190266/1‑100 (MQ=255)
                                                                                        aaTACATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat             >  1:734123/1‑100 (MQ=255)
                                                                                          tACATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGG           >  1:863098‑M1/1‑100 (MQ=255)
                                                                                           aCATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGATTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGC          >  1:2296257‑M1/1‑100 (MQ=255)
                                                                                           aCATGATTAATATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGC          >  2:1902061‑M1/1‑100 (MQ=255)
                                                                                            cATGATTGATATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCG         >  2:265710‑M1/1‑100 (MQ=255)
                                                                                              tgattaatATCAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTC       >  2:1707775‑M1/1‑100 (MQ=255)
                                                                                                   gatATTAATCATGTATTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCt  >  2:2174726/1‑100 (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 < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: