Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,573,528 A→G 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.2% 450.2 / ‑5.4 125pseudogene (1115/1503 nt)yjiTpseudogene
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (51/73);  total (51/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 126 total aligned reads displayed.

AGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCA  >  NC_000913/4573434‑4573623
                                                                                              |                                                                                               
aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTg                                                                                            >  2:1294352/1‑100 (MQ=255)
 gCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGa                                                                                           <  2:135928/100‑1 (MQ=255)
 gCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGa                                                                                           >  1:409844/1‑100 (MQ=255)
  cGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaa                                                                                          <  1:2130641/100‑1 (MQ=255)
   gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                         <  1:2618344/100‑1 (MQ=255)
   gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                         <  2:2842270/100‑1 (MQ=255)
   gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGaaa                                                                                         <  1:2931311/100‑1 (MQ=255)
    aaGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAAt                                                                                        >  1:1301845/1‑100 (MQ=255)
      gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                      <  1:386786/100‑1 (MQ=255)
      gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                      <  1:192428/100‑1 (MQ=255)
      gTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTg                                                                                      >  1:316087/1‑100 (MQ=255)
        gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAAATTGATGGCAGTGAAATTGcc                                                                                    >  1:475556/1‑100 (MQ=255)
          gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                  >  1:955287/1‑100 (MQ=255)
          gaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAt                                                                                  <  1:131962/100‑1 (MQ=255)
           aTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAtt                                                                                 <  2:2499782/100‑1 (MQ=255)
           aTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCAtt                                                                                 <  2:1572303/100‑1 (MQ=255)
              tGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGt                                                                              <  1:2727712/100‑1 (MQ=255)
               ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                             <  2:1700714/100‑1 (MQ=255)
               ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                             >  1:488676/1‑100 (MQ=255)
               ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTg                                                                             >  2:1199326/1‑100 (MQ=255)
                 caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTggg                                                                           >  1:1770815/1‑98 (MQ=255)
                 caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAg                                                                           <  2:1494667/100‑1 (MQ=255)
                 caggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAg                                                                           >  2:905961/1‑100 (MQ=255)
                  aggcaAAATCCATCGGTAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAgg                                                                          <  2:2816787/100‑1 (MQ=255)
                  aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAgg                                                                          <  1:2879135/100‑1 (MQ=255)
                    gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                        <  2:2123059/100‑1 (MQ=255)
                    gcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTg                                                                        <  1:3032395/100‑1 (MQ=255)
                     caAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGttc                                                                       >  1:289045/1‑98 (MQ=255)
                      aaaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGcc                                                                      <  1:1727703/100‑1 (MQ=255)
                        aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGGGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGt                                                                    >  1:1018269/1‑100 (MQ=255)
                        aaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGt                                                                    >  1:516293/1‑100 (MQ=255)
                          tCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTa                                                                  <  1:1931370/100‑1 (MQ=255)
                            cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAc                                                                <  2:855070/100‑1 (MQ=255)
                            cATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTAGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAc                                                                >  1:2896177/1‑100 (MQ=255)
                             aTCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                               <  1:2496913/100‑1 (MQ=255)
                             aTCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                               <  2:1182500/100‑1 (MQ=255)
                             aTCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAAcc                                                               <  2:2868484/100‑1 (MQ=255)
                              tCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCt                                                              <  2:1305466/100‑1 (MQ=255)
                               cGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCtt                                                             <  1:6754/100‑1 (MQ=255)
                                ggCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGGGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTc                                                            >  2:2827392/1‑100 (MQ=255)
                                 gCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTc                                                            >  1:3050988/1‑99 (MQ=255)
                                   aGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGGGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATc                                                         <  2:1787380/100‑1 (MQ=255)
                                     gTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGt                                                       >  1:2023352/1‑100 (MQ=255)
                                      tttGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGtt                                                      <  1:1106619/100‑1 (MQ=255)
                                       ttGTCGTTAGTGGCTCGTGCTGGAGGGATTATTGTTGGGAGCATTAAACTTGATGGCATTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTg                                                     <  2:289038/100‑1 (MQ=255)
                                         gTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGAt                                                   >  2:836454/1‑100 (MQ=255)
                                         gTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGAt                                                   <  1:1427340/100‑1 (MQ=255)
                                          tCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATg                                                  <  2:406474/100‑1 (MQ=255)
                                            gTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCa                                                >  2:509562/1‑100 (MQ=255)
                                             gatctGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAg                                               <  1:113076/96‑1 (MQ=255)
                                              tAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTTAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGa                                              <  1:388270/100‑1 (MQ=255)
                                              tAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGa                                              <  1:1631782/100‑1 (MQ=255)
                                                gTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGGGGTGCCGTTAACCTTCATCGTTGATGCAGATc                                            <  2:1770795/100‑1 (MQ=255)
                                                gTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCAGTTAACCTTCATCGTTGATGCAGATc                                            >  2:1242757/1‑100 (MQ=255)
                                                   gCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAAt                                         <  1:661015/100‑1 (MQ=255)
                                                    cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                        <  1:166177/100‑1 (MQ=255)
                                                    cTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                        <  2:179169/100‑1 (MQ=255)
                                                     tCGTGCTGGAGGGATGATTGTTGGGAGCATTAGACTTGACGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATgg                                       <  2:2147180/100‑1 (MQ=255)
                                                      cGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                      <  1:528392/100‑1 (MQ=255)
                                                        tGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTg                                    >  2:119568/1‑100 (MQ=255)
                                                        gatctgAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTg                                    <  2:2305360/95‑1 (MQ=255)
                                                          cTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                  <  2:1428597/100‑1 (MQ=255)
                                                            ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                                >  1:2772384/1‑100 (MQ=255)
                                                            ggAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTAc                                >  1:1554432/1‑100 (MQ=255)
                                                             gAGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTgagatc                               >  1:2305390/1‑95 (MQ=255)
                                                              aGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAg                              <  2:2022152/100‑1 (MQ=255)
                                                              aGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAg                              >  1:588875/1‑100 (MQ=255)
                                                              aGGGATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAg                              >  2:1597446/1‑100 (MQ=255)
                                                               gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                             >  1:2400633/1‑100 (MQ=255)
                                                               gggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAgg                             >  1:1442777/1‑100 (MQ=255)
                                                                ggATGATTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAggg                            <  1:2495988/100‑1 (MQ=255)
                                                                 gatgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGa                           >  1:2340571/1‑100 (MQ=255)
                                                                   tgatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACa                         >  1:2633177/1‑100 (MQ=255)
                                                                   tgatTGTTGGGAGCATTAAACTTGATGGCAGGGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACa                         <  1:1873007/100‑1 (MQ=255)
                                                                    gatTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAg                        >  1:2758299/1‑100 (MQ=255)
                                                                     atTGTTGGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAgg                       <  2:270779/100‑1 (MQ=255)
                                                                      ttgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGc                      >  1:1838073/1‑100 (MQ=255)
                                                                       tgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTTCAGGGACAGGct                     <  1:2006192/100‑1 (MQ=255)
                                                                       tgttgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGct                     <  1:2979522/100‑1 (MQ=255)
                                                                          tgGGAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttct                  >  2:2762442/1‑100 (MQ=255)
                                                                           gggAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGcttctt                 >  2:3038691/1‑100 (MQ=255)
                                                                             gAGCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc               <  2:1401108/100‑1 (MQ=255)
                                                                               gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg             <  2:1008029/100‑1 (MQ=255)
                                                                               gCATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAg             >  2:2385308/1‑100 (MQ=255)
                                                                                cATTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgt            >  1:507867/1‑100 (MQ=255)
                                                                                 aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg           <  1:2015/100‑1 (MQ=255)
                                                                                 aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg           >  2:1874261/1‑100 (MQ=255)
                                                                                 aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg           <  2:1411664/100‑1 (MQ=255)
                                                                                 aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg           <  1:1270678/100‑1 (MQ=255)
                                                                                 aTTAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtg           <  2:1018258/100‑1 (MQ=255)
                                                                                  ttAAACTTGATGGCAGTGGAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgt            <  2:1715869/98‑1 (MQ=255)
                                                                                  ttAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAgtgt          >  1:2736047/1‑100 (MQ=255)
                                                                                   tAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGtgtg         <  1:871162/100‑1 (MQ=255)
                                                                                   tAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGtgtg         <  2:1197385/100‑1 (MQ=255)
                                                                                   tAAACTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGtgtg         <  1:1461317/100‑1 (MQ=255)
                                                                                    aaaCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGc        <  2:2896137/100‑1 (MQ=255)
                                                                                    aaaCTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGc        >  1:1325678/1‑100 (MQ=255)
                                                                                      aCTTGATTACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCGTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGa      <  2:287729/100‑1 (MQ=255)
                                                                                       cTTGATGGCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGaa     >  1:1347397/1‑100 (MQ=255)
                                                                                          gatggCAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagca  >  2:506537/1‑100 (MQ=255)
                                                                                              |                                                                                               
AGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGCTGGAGGGATGATTGTTGGGAGCATTAAACTTGATGACAGTGAAATTGCCATTGGTGAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCA  >  NC_000913/4573434‑4573623

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