Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,210,811 C→T 100% pseudogene (1712/2001 nt) yehQ → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,210,8110CT98.0% 336.5 / ‑2.3 102pseudogene (1712/2001 nt)yehQpseudogene
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (40/59);  total (41/61)
Fisher's exact test for biased strand distribution p-value = 5.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

Only 100 of 102 total aligned reads displayed.

GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCT  >  NC_000913/2210714‑2210906
                                                                                                 |                                                                                               
gTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGc                                                                                               >  1:172463/1‑100 (MQ=255)
  cAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTa                                                                                             <  2:494005/100‑1 (MQ=255)
  cAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTa                                                                                             <  2:2390875/100‑1 (MQ=255)
    gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGCAACTACGCCTGTGGCACTGCTAcc                                                                                           >  1:1516472/1‑100 (MQ=255)
    gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                           <  2:1773586/100‑1 (MQ=255)
    gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                           <  2:410955/100‑1 (MQ=255)
    gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                           <  2:1639597/100‑1 (MQ=255)
    gTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTAcc                                                                                           >  2:262471/1‑100 (MQ=255)
     ttttGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCt                                                                                          <  1:872185/100‑1 (MQ=255)
       ttgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTc                                                                                        <  1:262470/100‑1 (MQ=255)
       ttgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTc                                                                                        <  2:1275820/100‑1 (MQ=255)
       ttgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTc                                                                                        <  2:1393655/100‑1 (MQ=255)
        tgttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCt                                                                                       >  2:2692959/1‑100 (MQ=255)
          ttgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGc                                                                                     >  2:1565338/1‑100 (MQ=255)
           tgAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGcc                                                                                    <  2:1852487/100‑1 (MQ=255)
              aGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGc                                                                                 >  2:248018/1‑100 (MQ=255)
              aGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGc                                                                                 >  1:526450/1‑100 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGGTGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTCCCAg                                                                                  >  1:1595414/1‑98 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                >  2:2233358/1‑100 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                >  2:345350/1‑100 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                >  2:2177698/1‑100 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                >  2:1906591/1‑100 (MQ=255)
               gggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCCGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCg                                                                                >  2:1139292/1‑100 (MQ=255)
                ggCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGt                                                                               <  2:1858803/100‑1 (MQ=255)
                  cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGttt                                                                             >  2:2748666/1‑100 (MQ=255)
                  cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGttt                                                                             <  2:485930/100‑1 (MQ=255)
                  cAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGttt                                                                             <  1:2844494/100‑1 (MQ=255)
                     ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                          <  2:1347327/100‑1 (MQ=255)
                     ttGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGc                                                                          <  1:710650/100‑1 (MQ=255)
                      tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                         <  1:2743040/100‑1 (MQ=255)
                      tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                         <  2:1709595/100‑1 (MQ=255)
                      tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                         <  2:2522631/100‑1 (MQ=255)
                      tGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGcc                                                                         >  1:95656/1‑100 (MQ=255)
                       gACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCt                                                                        <  2:2362793/100‑1 (MQ=255)
                         cacTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgt                                                                      >  1:3323989/1‑100 (MQ=255)
                          acTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                     <  1:2608878/100‑1 (MQ=255)
                          acTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                     <  2:2409991/100‑1 (MQ=255)
                          acTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                     <  2:678411/100‑1 (MQ=255)
                          acTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCtgtg                                                                     <  2:790567/100‑1 (MQ=255)
                          acTCGAACCCCGGGTGATGATGACAAAAACCCGTGGCTGGGCGCTGGACGCAGAAACTACGCCTGTTGCACTGCTAACTTCTGGTAGCGTTTTGTCtgtg                                                                     <  1:1948865/100‑1 (MQ=255)
                            tGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGcc                                                                   <  2:874508/100‑1 (MQ=255)
                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                  >  2:1019269/1‑100 (MQ=255)
                             ggAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCg                                                                  <  2:171861/100‑1 (MQ=255)
                                aCCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCt                                                               >  1:766683/1‑100 (MQ=255)
                                aCCACGGGCGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCt                                                               >  1:746781/1‑100 (MQ=255)
                                  cACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTAc                                                             <  1:2904290/100‑1 (MQ=255)
                                   aCGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACt                                                            <  2:768940/100‑1 (MQ=255)
                                     gggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGc                                                          <  2:2254626/100‑1 (MQ=255)
                                     gggTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGc                                                          <  2:389573/100‑1 (MQ=255)
                                          atgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatc                                                     <  2:3088280/100‑1 (MQ=255)
                                          atgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatc                                                     <  2:3198084/100‑1 (MQ=255)
                                           tgatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCtcatca                                                    >  1:706485/1‑100 (MQ=255)
                                            gatgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAg                                                   >  1:605225/1‑100 (MQ=255)
                                              tgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGtt                                                 >  1:1034222/1‑100 (MQ=255)
                                              tgaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGtt                                                 >  2:1105404/1‑100 (MQ=255)
                                               gaCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTg                                                >  2:1042530/1‑100 (MQ=255)
                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                               <  1:1906578/100‑1 (MQ=255)
                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                               >  1:35434/1‑100 (MQ=255)
                                                aCAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGc                                               <  2:507791/100‑1 (MQ=255)
                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                             <  1:2748643/100‑1 (MQ=255)
                                                  aaaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTg                                             <  1:2149466/100‑1 (MQ=255)
                                                   aaaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGa                                            >  2:2979584/1‑100 (MQ=255)
                                                    aaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGAt                                           >  2:1676227/1‑100 (MQ=255)
                                                    aaaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGAt                                           <  2:3277981/100‑1 (MQ=255)
                                                     aaCCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACAGCTACCTTCTGGCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATa                                          >  1:858956/1‑100 (MQ=255)
                                                       cccGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACg                                        >  1:3196808/1‑100 (MQ=255)
                                                       cccGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACg                                        >  2:1055925/1‑100 (MQ=255)
                                                          gTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTg                                     >  1:1496344/1‑100 (MQ=255)
                                                          gTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTg                                     <  2:301590/100‑1 (MQ=255)
                                                            gCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGcc                                   <  1:1115912/100‑1 (MQ=255)
                                                              cTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAg                                 >  2:553467/1‑100 (MQ=255)
                                                                 ggCGCTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCg                              >  2:519385/1‑100 (MQ=255)
                                                                   cgcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCg                              <  1:2217496/98‑1 (MQ=255)
                                                                   cgcTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGtt                            <  2:2118316/100‑1 (MQ=255)
                                                                     cTGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTAc                          <  1:1792303/100‑1 (MQ=255)
                                                                      tGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACt                         <  1:1105401/100‑1 (MQ=255)
                                                                      tGGACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACt                         <  2:20917/100‑1 (MQ=255)
                                                                       ggACGCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACtt                        >  1:1792090/1‑100 (MQ=255)
                                                                        gACGCAGAAACTACGCCTGTGGCACTGCTTCCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTTATCAGTTGCTGATACGCTGCCAGGCGTTACTTa                       <  1:2033996/100‑1 (MQ=255)
                                                                           gCAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTc                    <  1:1505561/100‑1 (MQ=255)
                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                   <  1:3136313/100‑1 (MQ=255)
                                                                            cAGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                   >  1:3212611/1‑100 (MQ=255)
                                                                             aGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTc                    >  1:2107665/1‑98 (MQ=255)
                                                                             aGAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCa                   <  2:768038/99‑1 (MQ=255)
                                                                              gAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAAc                 <  2:2546715/100‑1 (MQ=255)
                                                                              gAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAAc                 <  1:982311/100‑1 (MQ=255)
                                                                              gAAACTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAAc                 <  1:1410021/100‑1 (MQ=255)
                                                                               aaaCTACGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACt                <  2:2054557/100‑1 (MQ=255)
                                                                                    aCGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTcc           <  1:1144454/100‑1 (MQ=255)
                                                                                     cGCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCa          >  1:195366/1‑100 (MQ=255)
                                                                                      gCCTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCAt         <  1:54703/100‑1 (MQ=255)
                                                                                       ccTGTGGCACTGCTACCTTCTGCCATCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa        <  1:3038481/100‑1 (MQ=255)
                                                                                       ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa        >  1:2511976/1‑100 (MQ=255)
                                                                                       ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa        >  2:1389303/1‑100 (MQ=255)
                                                                                       ccTGTGGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATa        <  1:519383/100‑1 (MQ=255)
                                                                                         tgtgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAAc      <  2:254022/100‑1 (MQ=255)
                                                                                         tgtgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAAc      >  1:895948/1‑100 (MQ=255)
                                                                                          gtgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACg     <  2:2413651/100‑1 (MQ=255)
                                                                                           tgGCACTGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACgg    <  2:1487157/100‑1 (MQ=255)
                                                                                             gCACTGCTACCTTCTGTCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCt  <  2:514770/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GTCAGTTTTGTTGAAGGGCAATTGACACTGGAACCACGGGTGATGATGACAAAAACCCGTGCCTGGGCGCTGGACGCAGAAACTACGCCTGTGGCACCGCTACCTTCTGCCAGCGTTTTGCCTGTGCCGTCTACTGCTCATCAGTTGCTGATACGCTGCCAGGCGTTACTTATTCAACTGCTCCATAACGGCT  >  NC_000913/2210714‑2210906

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: