Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,743,165 C→G 15.1% V26V (GTG→GTC ribC ← riboflavin synthase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,743,1650CG15.1% 426.7 / 58.1 151V26V (GTG→GTCribCriboflavin synthase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (61/67);  new base G (6/17);  total (67/84)
Fisher's exact test for biased strand distribution p-value = 6.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

Only 100 of 164 total aligned reads displayed.

GACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAG  >  NC_000913/1743066‑1743263
                                                                                                   |                                                                                                   
gACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTcc                                                                                                     >  1:2262147/1‑100 (MQ=255)
    tGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTccacca                                                                                                 >  2:2641149/1‑100 (MQ=255)
      gTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACcaca                                                                                               <  1:961400/100‑1 (MQ=255)
       ttGCCGCTAATTTCCGTCACGGTCAGGCAGCAACCGTTATCCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCcaccac                                                                                                <  1:4325348/98‑1 (MQ=255)
         gCCGTTAATTTCCGTCGCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGc                                                                                            <  1:1416952/100‑1 (MQ=37)
         gCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGc                                                                                            <  2:1716685/100‑1 (MQ=38)
          ccGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGct                                                                                           <  2:2620377/100‑2 (MQ=37)
          ccGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGct                                                                                           <  1:2949298/100‑2 (MQ=37)
          ccGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGct                                                                                           <  1:1377948/100‑2 (MQ=37)
          ccGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCg                                                                                          <  1:3207232/100‑1 (MQ=255)
           cGTTGATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGct                                                                                           <  1:997270/99‑2 (MQ=37)
           cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGctg                                                                                          <  2:2864203/100‑3 (MQ=37)
           cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGctg                                                                                          <  2:909155/100‑3 (MQ=37)
           cGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGt                                                                                         >  1:1536751/1‑100 (MQ=255)
            gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTa                                                                                        >  2:3706930/1‑100 (MQ=255)
            gTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTa                                                                                        >  2:172983/1‑100 (MQ=255)
            gTTAATTTCCGTCACGGCCAGGCTGCACGCGTTATGCGCCACGGACGCACCGGTTTCCAGGCCGTCCAGCATGTGGTGGAGTTACCCGACCACATGctg                                                                                          <  1:3720004/99‑3 (MQ=21)
             tgattttGCGTCACGGGCAGACAGGAACCGTTATCGGCCACGGAAGCACCGGGTTCCAGGCGGTCCAGTATGTGGTGGAGTTACCCGACCACATGctg                                                                                          <  1:2177897/94‑3 (MQ=12)
              tAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACg                                                                                      <  1:3615378/100‑1 (MQ=255)
               aaTTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACg                                                                                      >  1:1195499/1‑99 (MQ=255)
                aTTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGaa                                                                                    >  1:1311751/1‑100 (MQ=255)
                 tttCCGTCACGGTCAGGCAGCAGCCGTGATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGaa                                                                                    <  2:3913079/99‑1 (MQ=255)
                 tttCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGaaa                                                                                   >  2:70920/1‑100 (MQ=255)
                    ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAtt                                                                                <  1:108048/100‑1 (MQ=255)
                    ccGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAtt                                                                                <  1:1859289/100‑1 (MQ=255)
                     cGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAAttt                                                                               >  2:1565966/1‑100 (MQ=255)
                      gTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTg                                                                              <  2:125508/100‑1 (MQ=255)
                       tCACGTTCAGGCAGCAACCGTTATGCGCCACGGAAGCGCCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTgg                                                                             <  1:451939/100‑1 (MQ=255)
                       tCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTgg                                                                             >  2:2233768/1‑100 (MQ=255)
                       tCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTgg                                                                             >  1:4252402/1‑100 (MQ=255)
                       tCACGGTCAAGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTgg                                                                             <  1:4024551/100‑1 (MQ=255)
                         aCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGtt                                                                           <  1:2630971/100‑1 (MQ=255)
                         aCGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGt                                                                            >  1:3986110/1‑99 (MQ=255)
                          cGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGttt                                                                          >  1:2476137/1‑100 (MQ=255)
                          cGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGttt                                                                          >  2:1790882/1‑100 (MQ=255)
                             tCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTtctc                                                                       <  2:1344895/100‑1 (MQ=255)
                              cAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCg                                                                      <  2:2271371/100‑1 (MQ=255)
                              cAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGAAACGCCACCACATGCGTACGAAAATTTGGTTTCTCg                                                                      >  1:2975210/1‑100 (MQ=255)
                               aGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGt                                                                     >  2:2122247/1‑100 (MQ=255)
                                ggCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTc                                                                    <  1:3863457/100‑1 (MQ=255)
                                 gcagcaACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTCCTCGTCa                                                                   >  1:2366584/1‑100 (MQ=255)
                                  ctggaaCCGATATGCGCGGCGGAAGCCCCGGTTCCCAGGCCGTCCAGCCTGTGGGCGGGTGCCTCCACCACATGCGTACGAAAATGTCGtt                                                                           >  1:2495740/5‑91 (MQ=255)
                                    gcaACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATc                                                                <  1:2273658/100‑1 (MQ=255)
                                     caACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCg                                                               >  2:485741/1‑100 (MQ=255)
                                      aaCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGa                                                              <  2:3383942/100‑1 (MQ=255)
                                       aCCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGATCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTGCTCGTCAATCGac                                                             <  1:2185297/100‑1 (MQ=255)
                                          gTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACAcc                                                          <  2:3142642/100‑1 (MQ=255)
                                          gTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACAcc                                                          >  1:2045162/1‑100 (MQ=255)
                                            tATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAg                                                        <  2:108217/100‑1 (MQ=255)
                                                 gcCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTg                                                   <  1:1704872/100‑1 (MQ=255)
                                                  ccACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGc                                                  <  2:2900792/100‑1 (MQ=255)
                                                    aCGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                <  2:3329926/100‑1 (MQ=255)
                                                    aCGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                >  1:88958/1‑100 (MQ=255)
                                                    aCGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCgg                                                <  2:3034196/100‑1 (MQ=255)
                                                     cGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGt                                               <  2:3291073/100‑1 (MQ=255)
                                                          gCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCt                                          <  1:950211/100‑1 (MQ=255)
                                                          gCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCt                                          >  2:3540217/1‑100 (MQ=255)
                                                           cACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTg                                         >  1:3280926/1‑100 (MQ=255)
                                                           cACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTg                                         <  2:1948730/100‑1 (MQ=255)
                                                           cACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTg                                         >  2:1775416/1‑100 (MQ=255)
                                                           cACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTg                                         <  2:3115538/100‑1 (MQ=255)
                                                           cACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTg                                         >  1:3885108/1‑100 (MQ=255)
                                                            aCCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGt                                        <  2:2539138/100‑1 (MQ=255)
                                                             ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGt                                        >  2:270105/1‑99 (MQ=255)
                                                             ccGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTa                                       >  2:797725/1‑100 (MQ=255)
                                                                gTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACaa                                    <  1:70920/100‑1 (MQ=255)
                                                                 tttCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAAt                                   <  2:2905456/100‑1 (MQ=255)
                                                                   tCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcc                                 <  1:1843327/100‑1 (MQ=255)
                                                                   tCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcc                                 <  1:2238987/100‑1 (MQ=255)
                                                                     cAGGCCGTCCCGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGATACCAGTTTTGCGGTGCCCTGTACAATcccc                               <  1:2899750/100‑1 (MQ=255)
                                                                     cAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATcccc                               <  1:3395840/100‑1 (MQ=255)
                                                                      aGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                              <  1:797724/100‑1 (MQ=255)
                                                                      aGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                              >  2:2453216/1‑100 (MQ=255)
                                                                      aGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCg                              <  2:4290696/100‑1 (MQ=255)
                                                                        gCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTa                            >  1:544937/1‑100 (MQ=255)
                                                                         ccGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTaa                           >  1:98216/1‑100 (MQ=255)
                                                                            tCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACa                        >  2:1293283/1‑100 (MQ=255)
                                                                            gtacgCATGTGGTGGAGTTACCCGACCACATGCTGGACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAAc                         >  2:1870862/5‑100 (MQ=18)
                                                                             tacgCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGCCACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACAt                       >  1:1422367/4‑100 (MQ=25)
                                                                              cAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATa                      >  2:3340852/1‑100 (MQ=255)
                                                                              cAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATa                      <  2:1162033/100‑1 (MQ=255)
                                                                               aGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATAc                     >  2:2114392/1‑100 (MQ=255)
                                                                                gCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCACTCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACt                    >  1:578746/1‑100 (MQ=37)
                                                                                gCATGTGGTGGAGTTACCCGACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACt                    >  2:640503/1‑100 (MQ=38)
                                                                                 cATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACtt                   >  1:765001/1‑100 (MQ=255)
                                                                                  aTGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACtt                   <  1:1657503/99‑1 (MQ=255)
                                                                                  aTGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTa                  <  1:4237457/100‑1 (MQ=255)
                                                                                   tgtgGTCGGGTAACTCCACCACATGCGTACTAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTAc                 <  2:2577288/100‑1 (MQ=255)
                                                                                    gtgGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACa                <  2:2716714/100‑1 (MQ=255)
                                                                                      ggTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAAc              >  2:3910093/1‑100 (MQ=255)
                                                                                         cGGGTGACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTcc           >  1:410732/1‑100 (MQ=255)
                                                                                         cGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTcc           >  1:3231611/1‑100 (MQ=255)
                                                                                            gTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGa        <  1:3801035/100‑1 (MQ=255)
                                                                                            gTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTCTTGCGGTGCCCTGTACAATCACCCTAAACATACTTCCACCTCCTGa        >  1:3669913/1‑100 (MQ=255)
                                                                                             tAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGaa       >  1:1009582/1‑100 (MQ=255)
                                                                                              aaCTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGaaa      >  2:75903/1‑100 (MQ=255)
                                                                                               aCTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAAt     >  2:1435703/1‑100 (MQ=255)
                                                                                                cTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAAt     >  1:4318091/1‑99 (MQ=255)
                                                                                                  ccaccaCATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAg  <  2:2228807/100‑1 (MQ=255)
                                                                                                  ccaccaCATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAg  >  1:3126217/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GACATGGTTGCCGTTAATTTCCGTCACGGTCAGGCAGCAACCGTTATGCGCCACGGAAGCACCGGTTTCCAGGCCGTCCAGCATGTGGTCGGGTAACTCCACCACATGCGTACGAAAATTTGGTTTCTCGTCAATCGACACCAGTTTTGCGGTGCCCTGTACAATCCCCGTAAACATACTTACAACTCCTGAAATCAG  >  NC_000913/1743066‑1743263

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