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.

CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  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)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑t                                                                                      <  2:2248043/100‑1 (MQ=255)
               ttGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTCCTGTAAAAATACATGATTGATATTAATCATGTa                                                                                        >  1:2422074/1‑99 (MQ=255)
                tGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTA‑tt                                                                                     <  2:124587/100‑1 (MQ=255)
                tGATTTGGATTATGTTGGGCAGTCGGTAAAACATTTTAATCCTGAATATGTTCGTGGTTATCTTGCTGTGGAAATACATAAGTAATAGCAATCATGTa                                                                                        <  1:2425220/98‑1 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑ttt                                                                                    <  2:1481472/100‑1 (MQ=255)
                 gATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑ttt                                                                                    >  2:56214/1‑100 (MQ=255)
                   tttGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTA‑ttttt                                                                                  <  2:700835/100‑2 (MQ=255)
                     tGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAcc                                                                                <  2:761801/100‑1 (MQ=255)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTA‑tttttaca                                                                               <  2:1536490/100‑5 (MQ=38)
                      gggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTA‑tttttaca                                                                               <  2:1885464/100‑5 (MQ=38)
                       ggTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTa                                                                              <  2:1733747/100‑1 (MQ=255)
                        gTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTaa                                                                             >  2:1257592/1‑100 (MQ=255)
                        gTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTAATATCAATCATGTA‑tttttacagc                                                                             <  1:1081871/100‑7 (MQ=37)
                         ttATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAg                                                                            <  2:1633467/100‑1 (MQ=255)
                            tGTTGGGGGGTCGGTAAAACATTTTAATCCTGAATTTGTAAGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTGATCATGTA‑TTTTACCTAAg                                                                            <  1:780019/97‑1 (MQ=255)
                            tGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAgacg                                                                         >  2:1924808/1‑100 (MQ=255)
                            tGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAgacg                                                                         <  1:943267/100‑1 (MQ=255)
                             gTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAgacga                                                                        >  1:2201066/1‑100 (MQ=255)
                               tGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACa                                                                      <  2:1085894/100‑1 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACaa                                                                     >  1:2283468/1‑100 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACaa                                                                     >  1:1784866/1‑100 (MQ=255)
                                ggggAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACaa                                                                     <  2:2026323/100‑1 (MQ=255)
                                  ggatacGGTAAAACATTTTAATCCTGACTTTTTTCGTGGTTGTCTTGCTGTTAGTATTCCTGATTTATATTTATCACGTA‑TTTTGCTTAAGACGACaata                                                                   >  1:1975049/6‑98 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAc                                                                 >  1:288848/1‑100 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAc                                                                 <  2:2110318/100‑1 (MQ=255)
                                    aGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAc                                                                 <  2:386734/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATTCTGAATTTGTTTGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAcc                                                                <  1:1007148/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAcc                                                                >  2:1098377/1‑100 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAcc                                                                <  2:233541/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAcc                                                                <  2:1185792/100‑1 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAcc                                                                >  1:1309606/1‑100 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAcc                                                                >  2:1463317/1‑100 (MQ=255)
                                     gTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAAcc                                                                <  1:1292454/100‑1 (MQ=255)
                                       cGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCtt                                                              >  1:885906/1‑100 (MQ=255)
                                        ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCttt                                                             <  1:2011361/100‑1 (MQ=255)
                                        ggTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCttt                                                             >  1:1852569/1‑100 (MQ=255)
                                         gTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCtttt                                                            <  2:54261/100‑1 (MQ=255)
                                          tAAAACATTTTAATCCTGAATTTGTGCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCtttt                                                            <  1:1836704/99‑1 (MQ=255)
                                             aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAAc                                                         <  1:1616362/99‑1 (MQ=255)
                                             aaCATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACa                                                        <  1:2396586/100‑1 (MQ=255)
                                               cATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTAAAAAAAAAAAAAAAAACCTTTTAACAgg                                                      <  1:404842/100‑1 (MQ=255)
                                                 ttttAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGa                                                    >  1:1742851/1‑100 (MQ=255)
                                                   ttAATCCTCAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGt                                                  <  1:1258533/100‑1 (MQ=255)
                                                     aaTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGa                                                <  2:1964511/100‑1 (MQ=255)
                                                      aTCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGat                                               <  1:1259201/100‑1 (MQ=255)
                                                          tGAATTTGTTCGTTGTTGTCTTGATGTAAAAATACATGATTGATATTAATCATGTA‑TGTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATAtt                                            <  1:2097074/99‑1 (MQ=255)
                                                          tGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTg                                           <  2:1159581/100‑1 (MQ=255)
                                                           gAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGa                                          >  2:2309809/1‑100 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                        <  2:1295462/100‑1 (MQ=255)
                                                             aTTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATc                                        <  2:358259/100‑1 (MQ=255)
                                                               ttGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCtt                                      <  2:1480839/100‑1 (MQ=255)
                                                                 gTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCa                                    >  1:36500/1‑100 (MQ=255)
                                                                    cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                 >  2:1701740/1‑100 (MQ=255)
                                                                    cGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActc                                 >  1:782566/1‑100 (MQ=255)
                                                                     gTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCActct                                >  1:14890/1‑100 (MQ=255)
                                                                       ggTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGt                              <  2:219636/100‑1 (MQ=255)
                                                                         ttGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCa                            <  1:124656/100‑1 (MQ=255)
                                                                          tGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCat                           <  2:444869/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          <  1:965763/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          <  2:1379335/100‑1 (MQ=255)
                                                                           gTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCata                          >  1:1485711/1‑100 (MQ=255)
                                                                             gtcttaggtaaaaTACATGATTAATATCAATCATGTATTTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGAGATATTGATCTTCAATCTGTCata                          >  1:47979/10‑99 (MQ=18)
                                                                              ttGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtct                       >  1:2103696/1‑100 (MQ=255)
                                                                              ttGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtct                       >  1:987164/1‑100 (MQ=255)
                                                                               tGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATAtct                       >  1:2321372/1‑99 (MQ=255)
                                                                                gCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTcc                     >  1:2175967/1‑100 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   <  2:885909/100‑1 (MQ=255)
                                                                                  tGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCgg                   >  2:1051729/1‑100 (MQ=255)
                                                                                    tAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 <  1:779920/100‑1 (MQ=255)
                                                                                    tAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 <  2:2379186/100‑1 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTGCgg                   >  1:946274/3‑98 (MQ=38)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:668519/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  1:613542/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  2:40306/3‑100 (MQ=255)
                                                                                    gtaaaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTa                 >  2:654378/3‑100 (MQ=255)
                                                                                     aaaaaTACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTaa                >  1:820656/1‑100 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  1:1644183/1‑100 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  2:1077973/1‑100 (MQ=255)
                                                                                       aaaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAata              >  2:1190266/1‑100 (MQ=255)
                                                                                        aaTACATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAatat             >  1:734123/1‑100 (MQ=255)
                                                                                          tACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGG           >  1:863098‑M1/1‑100 (MQ=255)
                                                                                           aCATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGATTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGC          >  1:2296257‑M1/1‑100 (MQ=255)
                                                                                           aCATGATTAATATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGC          >  2:1902061‑M1/1‑100 (MQ=255)
                                                                                            cATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCG         >  2:265710‑M1/1‑100 (MQ=255)
                                                                                              tgattaatATCAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTC       >  2:1707775‑M1/1‑100 (MQ=255)
                                                                                                   gatATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCt  >  2:2174726/1‑100 (MQ=255)
                                                                                                   |                                                                                                    
CACAAATTAGTTCGCTTGATTTGGGTTATGTTGGGGAGTCGGTAAAACATTTTAATCCTGAATTTGTTCGTGGTTGTCTTGCTGTAAAAATACATGATTGATATTAATCATGTA‑TTTTACCTAAGACGACAAAAACCTTTTAACAGGGAGTGATATTGATCTTCACTCTGTCATATCTCCGGTAATATGGCGTCAGGCT  >  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: ‑