Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,573,528 A→G 100% pseudogene (1115/1503 nt) yjiT → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,573,5280AG99.8% 572.4 / ‑6.7 171pseudogene (1115/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (75/94);  total (75/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

AAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGC  >  NC_000913/4573431‑4573622
                                                                                                 |                                                                                              
   aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTg                                                                                           <  1:3548607/100‑1 (MQ=255)
     cGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaa                                                                                         >  1:1324592/1‑100 (MQ=255)
      gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                        >  1:3247319/1‑100 (MQ=255)
        aGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGAATTAAACTTGATGGCAGTGAAAtt                                                                                      >  1:951106/1‑100 (MQ=255)
         gTGAGATTTGTCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                     <  1:3918032/100‑1 (MQ=255)
         gTGAGATTTGGCAGGCAAAATTCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTAGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                     <  1:3647752/100‑1 (MQ=255)
         gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                     >  1:688224/1‑100 (MQ=255)
           gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                   <  1:1653386/100‑1 (MQ=255)
           gagaTTTGGCAGGCAAAATCCATCGGCAGGTTGGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGcc                                                                                   >  2:3275682/1‑100 (MQ=255)
             gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                 >  2:187692/1‑100 (MQ=255)
             gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                 >  1:2180444/1‑100 (MQ=255)
             gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                 >  2:1628381/1‑100 (MQ=255)
              aTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAtt                                                                                <  1:3381115/100‑1 (MQ=255)
                ttGGCATGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGTCTGTGAAATTGCCATTgg                                                                              <  2:285581/100‑1 (MQ=255)
                ttGGCAGGCAAAATCCATCGGCAGGTTTGTTGTTAGTGGTTCGTTCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTgg                                                                              <  1:3283908/100‑1 (MQ=255)
                ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTgg                                                                              <  2:125669/100‑1 (MQ=255)
                ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTgg                                                                              <  1:3681578/100‑1 (MQ=255)
                ttGGAAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTAGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTgg                                                                              <  1:2128561/100‑1 (MQ=255)
                  ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                            <  2:2858231/100‑1 (MQ=255)
                   gcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGa                                                                           <  1:1484780/100‑1 (MQ=255)
                   gcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGa                                                                           <  1:2604585/100‑1 (MQ=255)
                    caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTTCTGGGGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                            >  1:2607368/1‑98 (MQ=255)
                    caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAg                                                                          <  2:1397367/100‑1 (MQ=255)
                    caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAg                                                                          <  1:2112092/100‑1 (MQ=255)
                      ggcaAAATCCATCGGCAGGTTTGTCGTTAGTGTCTCGTGCTGGAGAGATCATTATTGGGAGCATTAATCTTGATGGCTGTGAAATTGCCAGTGGTGAgg                                                                         <  2:147006/99‑1 (MQ=255)
                       gcaAACTCCATCGGCAGGTTTGTCGTTAGCGGCTCGTGCTGGCGGGATGATTGTTCGCAGCACTAAACTCGATGGCCGTGAAATTGCCATTGCTGacg                                                                         >  2:684434/1‑96 (MQ=255)
                       gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                       >  1:1175099/1‑100 (MQ=255)
                       gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                       <  1:3147008/100‑1 (MQ=255)
                        ccaactCCATCGGCAGGTATGTCGTTAGTGGCTCGTGCTGGAGAGATAATTTTTGGGAGTATTAAACTTGATTGCAGTGAAATTGCCATTGTTGAgg                                                                         <  2:3692233/92‑1 (MQ=255)
                           aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGt                                                                   <  1:3527685/100‑1 (MQ=255)
                           aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTACCATTGGTGAGGTGCCGt                                                                   <  1:1863031/100‑1 (MQ=255)
                            aTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGtt                                                                  >  2:2232651/1‑100 (MQ=255)
                            aTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGtt                                                                  >  1:1972667/1‑100 (MQ=255)
                              ccATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTaa                                                                <  2:2195103/100‑1 (MQ=255)
                               cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAc                                                               <  1:2835187/100‑1 (MQ=255)
                                  cGGCAGTTTTGTTGTTAGTTGTTCGTTTTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCtt                                                            <  1:220618/100‑1 (MQ=255)
                                       ggTCTGTCGTTAGTGGCTCGTGCTGGAGCGATGCTTGTTGGGAGCATTAAACTTGATGGCCGTGACATTGCCCTTGGCGCGGTGCCGTTACCCTTCATCg                                                       >  2:1970232/1‑100 (MQ=255)
                                            gTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTTTTTGTAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGAt                                                  <  1:8518/100‑1 (MQ=255)
                                              cGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                <  2:3607151/100‑1 (MQ=255)
                                              cGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                >  2:3871555/1‑100 (MQ=255)
                                              cGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                >  2:3314591/1‑100 (MQ=255)
                                               gTTAGTGGCTCGTTCTGGAGGTATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                <  1:1717358/99‑1 (MQ=255)
                                               gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGAAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGc                                                <  1:490814/99‑1 (MQ=255)
                                                ttAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAg                                              <  2:137232/100‑1 (MQ=255)
                                                  aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                            >  1:3143635/1‑100 (MQ=255)
                                                  aGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGAt                                            <  1:1270504/100‑1 (MQ=255)
                                                   gTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATc                                           <  1:579855/100‑1 (MQ=255)
                                                   gTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATc                                           <  1:2243757/100‑1 (MQ=255)
                                                     ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCGATTGGTGAGGTGGCGTTAACCTTCATCGTTGATGCAGATCa                                          >  2:2470400/1‑99 (MQ=255)
                                                     ggCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCaa                                         >  2:2521577/1‑100 (MQ=255)
                                                      gCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAAt                                        <  2:384041/100‑1 (MQ=255)
                                                      gCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAAt                                        <  1:2878927/100‑1 (MQ=255)
                                                       cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       <  2:1223461/100‑1 (MQ=255)
                                                       cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       <  1:2984886/100‑1 (MQ=255)
                                                       cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       <  1:2053354/100‑1 (MQ=255)
                                                       cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       <  2:496923/100‑1 (MQ=255)
                                                       cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       <  2:3905844/100‑1 (MQ=255)
                                                       cTCGTGCTGGAGGGATGATTATTGGGAGCATTAAACTTGATGGCAGTGAAATTGCTATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       >  2:3506405/1‑100 (MQ=255)
                                                        tCGTTCTGGATGGATGATTGTTGGGAGCATTAAACTTTATGGCAGTTAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                       <  1:1299541/99‑1 (MQ=255)
                                                         cGTGCTGGAGGGATGATTGTTGGGATTATTAAACTTTATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                     <  1:108329/100‑1 (MQ=255)
                                                         cGTGCTGGAGGGATGATTGTTGGGAGCTTTAAACTTGTTTGCAGTGTATTTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                     <  1:372893/100‑1 (MQ=255)
                                                         cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGATGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                     <  2:1521667/100‑1 (MQ=255)
                                                         cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                     <  2:763159/100‑1 (MQ=255)
                                                         cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAAATGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                     <  1:3778602/100‑1 (MQ=255)
                                                          gTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCt                                    >  2:1991544/1‑100 (MQ=255)
                                                            gCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGATGt                                  >  1:1689078/1‑100 (MQ=255)
                                                              tGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGttn                                <  2:3006713/100‑2 (MQ=255)
                                                              tGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTa                                <  2:208155/100‑1 (MQ=255)
                                                               ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                               >  1:3119002/1‑100 (MQ=255)
                                                               ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                               <  1:1618467/100‑1 (MQ=255)
                                                               ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGACGTTAACATTAATCGTTGATGCAGATCAATGACTGTTAc                               >  1:903305/1‑100 (MQ=255)
                                                                  gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                            >  2:572503/1‑100 (MQ=255)
                                                                  gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                            >  1:3598838/1‑100 (MQ=255)
                                                                      tgatTGTTGGGAGCATTAAACTTGATGGCTGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCAGCGTTGATGCAGATCAATGGCTGTTACAGGGACa                        <  2:2347932/100‑1 (MQ=255)
                                                                        atTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAg                       >  1:1106041/1‑99 (MQ=255)
                                                                        atTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAAAGGGACAg                       >  2:899037/1‑99 (MQ=255)
                                                                         ttgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGc                     >  1:2387187/1‑100 (MQ=255)
                                                                         ttgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGc                     <  1:2838208/100‑1 (MQ=255)
                                                                             tgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttct                 >  1:1631264/1‑100 (MQ=255)
                                                                              gggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt                >  2:215507/1‑100 (MQ=255)
                                                                              gggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt                >  2:147993/1‑100 (MQ=255)
                                                                              gggAGCATTAAACTTGATGGAAGTGAAATTTCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATAAATGGCTGTTACAGGGACAGGcttctt                <  2:1647164/100‑1 (MQ=255)
                                                                              gggAGCATCAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAg                       >  2:258089/1‑93 (MQ=255)
                                                                                gAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc              <  2:444611/100‑1 (MQ=255)
                                                                                  gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg            <  2:2722587/100‑1 (MQ=255)
                                                                                  gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg            >  1:2812032/1‑100 (MQ=255)
                                                                                   cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgt           >  1:407568/1‑100 (MQ=255)
                                                                                   cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgt           >  2:3309005/1‑100 (MQ=255)
                                                                                   cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgt           >  2:2496270/1‑100 (MQ=255)
                                                                                   cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg            <  1:544455/99‑1 (MQ=255)
                                                                                    aTTAAACTTGATGGCAGTGAAATTGCGATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgt           >  1:2386834/1‑99 (MQ=255)
                                                                                      tAAACTTGATGGCAGTGAAATTGCTATTGGTGAGTTGCTGTTAATCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGtgtg        <  1:13946/100‑1 (MQ=255)
                                                                                      tAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGtgtg        <  1:2564641/100‑1 (MQ=255)
                                                                                        aaCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCg      >  1:1547511/1‑100 (MQ=255)
                                                                                         aCTTGATTGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGa     >  2:1620987/1‑100 (MQ=255)
                                                                                          cTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGaa    <  1:533194/100‑1 (MQ=255)
                                                                                          cTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGaa    >  1:3297740/1‑100 (MQ=255)
                                                                                           ttGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAag   >  2:3388766/1‑100 (MQ=255)
                                                                                            tgatggCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagc  >  1:2746070/1‑100 (MQ=255)
                                                                                            tgatggCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagc  >  1:397219/1‑100 (MQ=255)
                                                                                                 |                                                                                              
AAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGC  >  NC_000913/4573431‑4573622

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: