Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,257,774 G→A 22.5% C198C (TGC→TGT psuG ← pseudouridine 5'‑phosphate glycosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,257,7740GA22.5% 238.0 / 65.9 120C198C (TGC→TGTpsuGpseudouridine 5'‑phosphate glycosidase
Reads supporting (aligned to +/- strand):  ref base G (42/51);  new base A (10/17);  total (52/68)
Fisher's exact test for biased strand distribution p-value = 5.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 125 total aligned reads displayed.

AGGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGATTTC  >  NC_000913/2257676‑2257873
                                                                                                  |                                                                                                   
aGGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGc                                                                                                    <  1:1062464/100‑1 (MQ=255)
aGGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGc                                                                                                    <  1:2200725/100‑1 (MQ=255)
 ggCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCa                                                                                                   <  1:1246735/100‑1 (MQ=255)
  gCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAg                                                                                                  <  1:68423/100‑1 (MQ=255)
   ccACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGa                                                                                                 >  2:899055/1‑100 (MQ=255)
    cACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGa                                                                                                 <  2:2846523/99‑1 (MQ=255)
     aCCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCa                                                                                                   <  2:2604091/96‑1 (MQ=255)
      ccGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGaaaa                                                                                              >  2:2696474/1‑100 (MQ=255)
         atctgACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACg                                                                                           <  2:3301533/96‑1 (MQ=255)
           cAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACgca                                                                                         <  2:889265/100‑1 (MQ=255)
             gACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACgagat                                                                                       >  1:3301509/1‑96 (MQ=255)
              aCCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACgcaggc                                                                                      <  1:434813/100‑1 (MQ=255)
               ccGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACgcaggca                                                                                     >  1:1674714/1‑100 (MQ=255)
               ccGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACgcaggca                                                                                     <  2:2443864/100‑1 (MQ=255)
                 gCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAgc                                                                                   >  1:1635455/1‑100 (MQ=255)
                     tGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCtt                                                                               >  2:277799/1‑100 (MQ=255)
                      gCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCttt                                                                              <  2:1499042/100‑1 (MQ=255)
                      gCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCttt                                                                              <  1:1155008/100‑1 (MQ=255)
                      gCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGGCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCttt                                                                              <  1:781097/100‑1 (MQ=255)
                       ccATTTCACCACCATTGCGCGGGCAATTTCGCTGGCGCTGTCTAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTa                                                                             >  2:1428503/1‑100 (MQ=255)
                       ccATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTa                                                                             <  1:598096/100‑1 (MQ=255)
                        cATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAg                                                                            <  2:620384/100‑1 (MQ=255)
                        cATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAg                                                                            >  2:387059/1‑100 (MQ=255)
                         aTTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGt                                                                           >  2:1342656/1‑100 (MQ=255)
                          tttCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTc                                                                          <  2:710313/100‑1 (MQ=255)
                           ttCACCACCATTGTACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCt                                                                         <  2:3412849/100‑1 (MQ=255)
                           ttCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCt                                                                         >  2:1609137/1‑100 (MQ=255)
                           ttCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCt                                                                         <  1:1548071/100‑1 (MQ=255)
                           ttCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCt                                                                         <  1:275506/100‑1 (MQ=255)
                            tCACCACCATTGCACGGGCAATTTCGATGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTg                                                                        >  1:1006881/1‑100 (MQ=255)
                             caccacCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGa                                                                       >  1:2536441/1‑100 (MQ=255)
                              accaccATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGAt                                                                      <  1:681959/100‑1 (MQ=255)
                               ccaccaTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATa                                                                     >  2:1423233/1‑100 (MQ=255)
                               ccaccaTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATa                                                                     >  1:1719370/1‑100 (MQ=255)
                                  ccaTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGcc                                                                  <  1:1054696/100‑1 (MQ=255)
                                   caTTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCa                                                                 <  2:3302807/100‑1 (MQ=255)
                                     ttGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAAt                                                               >  2:3260022/1‑100 (MQ=255)
                                     ttGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAAt                                                               <  2:3359528/100‑1 (MQ=255)
                                      tGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAAtt                                                              >  2:106794/1‑100 (MQ=255)
                                       gCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTa                                                             >  2:53002/1‑100 (MQ=255)
                                        cACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTaa                                                            <  2:2637124/100‑1 (MQ=255)
                                          cGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACg                                                          >  2:2560904/1‑100 (MQ=255)
                                          cGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACg                                                          <  2:1901293/100‑1 (MQ=255)
                                            ggCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGc                                                        <  2:1635471/100‑1 (MQ=255)
                                            ggCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGc                                                        >  1:2082608/1‑100 (MQ=255)
                                              cAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGcac                                                      <  1:2696453/100‑1 (MQ=255)
                                                aTTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCacac                                                    >  1:2615310/1‑100 (MQ=255)
                                                 tttCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACAcc                                                   >  2:3095803/1‑100 (MQ=255)
                                                    cGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGaa                                                >  1:2474320/1‑100 (MQ=255)
                                                    cGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGaa                                                >  2:1859223/1‑100 (MQ=255)
                                                     gCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAg                                               >  2:17655/1‑100 (MQ=255)
                                                      cTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAgg                                              >  2:568608/1‑100 (MQ=255)
                                                        ggCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGtt                                            >  2:317576/1‑100 (MQ=255)
                                                          cgcTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTc                                          >  2:1100495/1‑100 (MQ=255)
                                                            cTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTa                                        <  1:3385337/100‑1 (MQ=255)
                                                            cTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTa                                        <  2:1540960/100‑1 (MQ=255)
                                                            cTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTa                                        >  2:859155/1‑100 (MQ=255)
                                                               tCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAAt                                     <  1:2157832/100‑1 (MQ=255)
                                                               tCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAAt                                     <  1:1710456/100‑1 (MQ=255)
                                                                cGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGAAGGCAGGGATTTAGTCTGATAGCCAATTAACGGAACACCGAAGGTTTCTAAAt                                     >  1:686380/1‑99 (MQ=255)
                                                                 gagaCGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATAc                                   <  2:178013/100‑1 (MQ=255)
                                                                 gagaCGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATAc                                   <  1:2719263/100‑1 (MQ=255)
                                                                  agaCGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACt                                  <  1:2480642/100‑1 (MQ=255)
                                                                  agaCGAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACt                                  <  1:1140932/100‑1 (MQ=255)
                                                                   gaCGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTc                                 <  2:2355178/100‑1 (MQ=255)
                                                                    aagAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCa                                <  1:1609121/98‑1 (MQ=255)
                                                                     cGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAg                               >  2:2669140/1‑100 (MQ=255)
                                                                      gAATGCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGt                              >  1:868136/1‑100 (MQ=255)
                                                                        aTGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTgg                            <  2:862765/100‑1 (MQ=255)
                                                                          gCTGACGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGtt                          <  1:3277236/100‑1 (MQ=255)
                                                                              aCGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATc                      <  2:3390672/100‑1 (MQ=255)
                                                                              aCGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATc                      <  2:708362/100‑1 (MQ=255)
                                                                               cGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATcc                     <  2:2260906/100‑1 (MQ=255)
                                                                               cGTCAAACGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATcc                     <  2:3249432/100‑1 (MQ=255)
                                                                                 tCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTa                   >  1:2285218/1‑100 (MQ=255)
                                                                                  ctaaCGGGCTGGTACGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTaa                  >  2:2344668/3‑100 (MQ=255)
                                                                                  cAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTaa                  <  2:1082035/100‑1 (MQ=255)
                                                                                  cAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAATCTGATAGCCAATTAACGGGACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTaa                  >  1:1893929/1‑100 (MQ=255)
                                                                                   aaaCGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTaaa                 <  1:1859207/100‑1 (MQ=255)
                                                                                    aaCGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAAt                <  2:1326491/100‑1 (MQ=255)
                                                                                      cGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCg              <  1:1406169/100‑1 (MQ=255)
                                                                                      cGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCg              >  1:428331/1‑100 (MQ=255)
                                                                                      cGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCg              <  1:2956118/100‑1 (MQ=255)
                                                                                      cGGGCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCg              <  1:2037873/100‑1 (MQ=255)
                                                                                        ggCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgag            <  2:2615329/100‑1 (MQ=255)
                                                                                        ggCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgag            <  1:1428491/100‑1 (MQ=255)
                                                                                        ggCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgag            <  1:2777291/100‑1 (MQ=255)
                                                                                        ggCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGGACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgag            <  1:473765/100‑1 (MQ=255)
                                                                                         gCTGGTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgaga           <  2:548685/100‑1 (MQ=255)
                                                                                         gCTCGTGCAGCAGAAACACGCAGGCGGCGCTTTAGTCTGATAGCCAATTAACGGAACACCGAAGGTGTCTAAACACTCAGTGGTTAATCCTAAAt                <  1:899048/95‑1 (MQ=255)
                                                                                          cTGGTGCGGCAGAAAAAAGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGAACACCGAAGGTTTCTAAATACACAGTGGTTTATCCTAAATCGAGaa          >  1:1101014/1‑100 (MQ=255)
                                                                                           tGGTGCGGGAAAAAAACGCAGGCAGAGCTTTAGTCTGATAGCCAATTAAAGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCgag            <  1:770533/97‑1 (MQ=255)
                                                                                            ggTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATa        >  2:300602/1‑100 (MQ=255)
                                                                                             gTGCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAg       <  2:3122536/100‑1 (MQ=255)
                                                                                              tGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGa      >  1:66344/1‑100 (MQ=255)
                                                                                               gCGACAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGAt     >  2:624677/1‑100 (MQ=255)
                                                                                                 ggCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGAttt   >  1:1146026/1‑100 (MQ=255)
                                                                                                 ggCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGAttt   >  1:1962526/1‑100 (MQ=255)
                                                                                                  gCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGATTTc  <  2:1146037/100‑1 (MQ=255)
                                                                                                  gCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGATTTc  >  2:479379/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
AGGCCACCGTTCAGACCGCTTTGCCATTTCACCACCATTGCACGGGCAATTTCGCTGGCGCTGTCGAGACGAATGCTGACGTCAAACGGGCTGGTGCGGCAGAAAAACGCAGGCAGCGCTTTAGTCTGATAGCCAATTAACGGCACACCGAAGGTTTCTAAATACTCAGTGGTTAATCCTAAATCGAGAATAGATTTC  >  NC_000913/2257676‑2257873

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