Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 529,063 A→G 17.8% pseudogene (424/491 nt) ybbD → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913529,0630AG17.8% 295.5 / 52.7 142pseudogene (424/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base A (48/65);  new base G (12/12);  total (65/77)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 144 total aligned reads displayed.

AAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGA  >  NC_000913/528964‑529161
                                                                                                   |                                                                                                  
   tCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAAt                                                                                                 <  1:1561266/100‑1 (MQ=255)
     gTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTa                                                                                               >  1:368490/1‑100 (MQ=255)
      tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTaa                                                                                              <  1:416408/100‑1 (MQ=255)
       aTCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAAt                                                                                             <  2:789232/100‑1 (MQ=255)
       aTCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAAt                                                                                             <  2:2162000/100‑1 (MQ=255)
        tCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATg                                                                                            <  2:1202882/100‑1 (MQ=255)
        tCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATg                                                                                            >  1:98599/1‑100 (MQ=255)
         cAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGt                                                                                           >  2:2129285/1‑100 (MQ=255)
          aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGtt                                                                                          >  1:166979/1‑100 (MQ=255)
          aaaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGt                                                                                           <  1:1886773/99‑1 (MQ=255)
           aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTa                                                                                         <  1:144900/100‑1 (MQ=255)
           aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                         >  2:1028325/1‑100 (MQ=255)
           aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                         <  1:537231/100‑1 (MQ=255)
            ataaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTAc                                                                                        <  2:2256617/98‑1 (MQ=255)
             gAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACg                                                                                       >  1:2211629/1‑100 (MQ=255)
                 gATCTATATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCa                                                                                     <  2:529322/98‑1 (MQ=255)
                  aTCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTCTAAACCTTTTTTTTGTAAATAAAATTAAAATTACACatta                                                                                  >  1:14902/1‑100 (MQ=255)
                   tCCAGATTCAAAAAAAGTGTGGTATAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACTTTTTTTTTGTTGGTTTAATTAATATTACGCattat                                                                                 >  2:31641/1‑100 (MQ=255)
                     cAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTAttt                                                                               >  1:181181/1‑100 (MQ=255)
                     cAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTAttt                                                                               <  1:298396/100‑1 (MQ=255)
                     cAGATACAAAAAAAGTGTGGCGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTAttt                                                                               >  2:1461559/1‑100 (MQ=255)
                      aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGTGGATTGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCattat                                                                                 <  2:20700/97‑1 (MQ=255)
                      aGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTAtttt                                                                              <  1:281601/100‑1 (MQ=255)
                       gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTta                                                                             <  1:1697882/100‑1 (MQ=255)
                           cAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTTTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                                                         >  2:2388902/1‑100 (MQ=255)
                           cAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                                                         >  2:1503214/1‑100 (MQ=255)
                            aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACTTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                                                         <  2:779301/99‑1 (MQ=255)
                               aaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAat                                                                     <  1:616224/100‑1 (MQ=255)
                               aaaaGTGTAGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAat                                                                     <  1:1520280/100‑1 (MQ=255)
                                    tgtgGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                                                                <  1:1993166/100‑1 (MQ=255)
                                     gtggtgTAAATATGGTAAGATACCAGGGCAAGAGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGa                                                               <  2:176760/100‑1 (MQ=255)
                                      tggtgTAAATATGGTAAGTTATCAGGGAAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAg                                                              <  1:173382/100‑1 (MQ=255)
                                       ggtgTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGc                                                             >  1:1894081/1‑100 (MQ=255)
                                         tgtAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTg                                                           <  2:617810/100‑1 (MQ=255)
                                          gtAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTTTTTTTTTTATTAATGTTACGCATTATTTTATAACAAATATTGTAGCTg                                                           >  2:422431/1‑99 (MQ=255)
                                          gtAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTTTTTGTTTAATTAATGTTACGCATTATTTTATAACAAATATTGTAGCTg                                                           >  2:348806/1‑99 (MQ=255)
                                          gtAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTgg                                                          <  2:1901107/100‑1 (MQ=255)
                                           tAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGa                                                         <  2:1447208/100‑1 (MQ=255)
                                            aaaTATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                        <  1:1086435/100‑1 (MQ=255)
                                            aaaTATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                        >  1:815397/1‑100 (MQ=255)
                                            aaaTATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                        <  1:6835/100‑1 (MQ=255)
                                              atatGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTg                                                      <  1:2217423/100‑1 (MQ=255)
                                                 tGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                                                   <  1:1088836/100‑1 (MQ=255)
                                                  ggTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                  >  1:2071251/1‑100 (MQ=255)
                                                   gTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGa                                                 <  1:2383287/100‑1 (MQ=255)
                                                   gTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGa                                                 >  1:1463276/1‑100 (MQ=255)
                                                   gTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTAAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                                                    <  2:17839/97‑1 (MQ=255)
                                                    tAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGAt                                                <  2:2038376/100‑1 (MQ=255)
                                                    tAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGAt                                                <  2:2037687/100‑1 (MQ=255)
                                                     aaGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATg                                               <  2:1322875/100‑1 (MQ=255)
                                                      aGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGc                                              <  2:1143161/100‑1 (MQ=255)
                                                         tACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGAAGCtt                                            >  1:1195608/1‑99 (MQ=255)
                                                           ccAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTTGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg                                         <  2:1496926/100‑1 (MQ=255)
                                                           ccAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg                                         <  2:2319870/100‑1 (MQ=255)
                                                           ccAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg                                         >  1:352280/1‑100 (MQ=255)
                                                           ccAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg                                         >  1:329093/1‑100 (MQ=255)
                                                           ccAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg                                         >  1:1291640/1‑100 (MQ=255)
                                                           ccAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg                                         >  1:1524070/1‑100 (MQ=255)
                                                            cAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGc                                        <  2:2050077/100‑1 (MQ=255)
                                                             aGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCa                                       >  1:318230/1‑100 (MQ=255)
                                                              gggCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                      >  2:1918750/1‑100 (MQ=255)
                                                              gggCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                      >  1:1671847/1‑100 (MQ=255)
                                                              gggCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                      <  1:1515039/100‑1 (MQ=255)
                                                               ggCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCaga                                     >  1:1063161/1‑100 (MQ=255)
                                                               ggCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                      <  1:677434/99‑1 (MQ=255)
                                                                 cAAGGGGATGGTGTAAACCTTTTTTTTTTTTTTTAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTTGATTGCCTGATGCTTTTGCaga                                     >  2:405597/1‑98 (MQ=255)
                                                                     gggATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAtt                               <  1:1718433/100‑1 (MQ=255)
                                                                     gggATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAtt                               <  2:2062957/100‑1 (MQ=255)
                                                                      ggATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTTGAGCTGGATTGCCTGATGCTTTTGCAGCATAATTg                              >  2:94671/1‑100 (MQ=255)
                                                                       gATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGc                             >  1:902259/1‑100 (MQ=255)
                                                                       gATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGc                             <  2:1984016/100‑1 (MQ=255)
                                                                        aTGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCt                            <  2:269550/100‑1 (MQ=255)
                                                                         tGGTGTAAACCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCtt                           >  2:1631180/1‑100 (MQ=255)
                                                                         tGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCtt                           <  2:1585062/100‑1 (MQ=255)
                                                                          ggTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGGGCTGGATTGCCTGATGCTTGTACAGAGaaa                                 >  1:1253769/1‑93 (MQ=255)
                                                                           gtgtAAACCTTTTTTTTTTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGa                         >  2:1121720/1‑100 (MQ=255)
                                                                            tgtAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGaa                        >  1:2328853/1‑100 (MQ=255)
                                                                            tgtAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGa                         <  2:1110950/99‑1 (MQ=255)
                                                                            tgtAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCAGTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGaa                        <  1:1077641/100‑1 (MQ=255)
                                                                              tAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGaa                        <  1:910043/98‑1 (MQ=255)
                                                                               aaaCCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTa                     <  2:524067/100‑1 (MQ=255)
                                                                               aaaCCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTa                     <  1:1938383/100‑1 (MQ=255)
                                                                                aaCCTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAAAAAATATTGGAGCTGGATTGCTTGATGGTTGTGAAGAGTAATTGCTTGACtt                      >  1:1034296/1‑98 (MQ=255)
                                                                                   cTTTTTTTTGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTAAATTAagag                 >  1:1880590/1‑100 (MQ=255)
                                                                                     tttttttGTTGGTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTc               <  1:58915/100‑1 (MQ=255)
                                                                                       tttttGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTa             <  2:17627/100‑1 (MQ=255)
                                                                                         tttGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATc           <  1:1723950/100‑1 (MQ=255)
                                                                                           tgttgGTGAAATTAATGTTACTCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATc           <  2:2283613/98‑1 (MQ=255)
                                                                                           tgttgGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATccc         >  2:925146/1‑100 (MQ=255)
                                                                                           tgttgGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATcc          >  2:924894/1‑99 (MQ=255)
                                                                                            gttgGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCa        <  1:2071110/100‑1 (MQ=255)
                                                                                              tgGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCata      >  1:1834587/1‑100 (MQ=255)
                                                                                               ggTGGAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCatat     >  2:963796/1‑100 (MQ=255)
                                                                                               ggTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCatat     <  1:560970/100‑1 (MQ=255)
                                                                                               ggTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCatat     >  2:1168729/1‑100 (MQ=255)
                                                                                                gTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATc    >  1:117107/1‑100 (MQ=255)
                                                                                                gTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATc    <  1:1355906/100‑1 (MQ=255)
                                                                                                 tGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCg   >  2:1167904/1‑100 (MQ=255)
                                                                                                 tGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGAAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATc    >  2:1612135/1‑99 (MQ=255)
                                                                                                  gAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGa  >  1:2219208/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
AAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGA  >  NC_000913/528964‑529161

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: ‑