Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,596,781 A→T 7.6% K264N (AAA→AAT lgoD → L‑galactonate oxidoreductase; L‑gulonate oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,596,7810AT 7.6% 2327.9 / 26.0 763K264N (AAA→AATlgoDL‑galactonate oxidoreductase; L‑gulonate oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (428/222);  new base T (4/65);  total (432/331)
Fisher's exact test for biased strand distribution p-value = 1.53e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 789 total aligned reads displayed.

GCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTA  >  NC_000913/4596632‑4596930
                                                                                                                                                     |                                                                                                                                                     
gCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTaaa                                                                                                                                                       <  1:118882/150‑1 (MQ=255)
gCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTaaa                                                                                                                                                       <  1:12726209/150‑1 (MQ=255)
 ccGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAg                                                                                                                                                      <  2:10785206/150‑1 (MQ=255)
 ccGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAg                                                                                                                                                      >  1:13433737/1‑150 (MQ=255)
 ccGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAg                                                                                                                                                      >  1:9556137/1‑150 (MQ=255)
   gAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGt                                                                                                                                                    >  2:13902815/1‑150 (MQ=255)
   gAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGt                                                                                                                                                    >  1:14858413/1‑150 (MQ=255)
    aaGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGTTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTg                                                                                                                                                   >  1:2585698/1‑150 (MQ=255)
    aaGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTg                                                                                                                                                   >  2:12923837/1‑150 (MQ=255)
     aGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGa                                                                                                                                                  <  2:2261249/150‑1 (MQ=255)
     aGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGa                                                                                                                                                  <  1:14175616/150‑1 (MQ=255)
      gATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAg                                                                                                                                                 >  2:779457/1‑150 (MQ=255)
      gATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAg                                                                                                                                                 >  1:3418649/1‑150 (MQ=255)
       aTTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGt                                                                                                                                                >  2:3128872/1‑150 (MQ=255)
          ttGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGc                                                                                                                                             >  1:3024479/1‑150 (MQ=255)
           tGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCa                                                                                                                                            <  2:131387/150‑1 (MQ=255)
              cgcgCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGtt                                                                                                                                         <  1:5786564/150‑1 (MQ=255)
               gcgcAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTtc                                                                                                                                        >  1:12378013/1‑150 (MQ=255)
               gcgcAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTtc                                                                                                                                        >  2:6200305/1‑150 (MQ=255)
               gcgcAGCTACGGGCGCAGTTTGGTGGATAGCTGGCGAAGAAAGTGATCGACGCGACATGTAATCAACAATCGATGAAAAACACCGTGAATTTGATTAGACACGGCGGCACGGTTGTATTTGTCTGCATGTTTAAAGGTGAGTTTCAGtta                                                                                                                                        >  1:6916253/1‑149 (MQ=255)
                cgcAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTtct                                                                                                                                       >  1:9817221/1‑150 (MQ=255)
                 gcAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTtctc                                                                                                                                      >  1:14625513/1‑150 (MQ=255)
                 gcAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTtctc                                                                                                                                      <  2:13449131/150‑1 (MQ=255)
                 gcAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTtctc                                                                                                                                      <  2:5927204/150‑1 (MQ=255)
                 gcAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTtctc                                                                                                                                      <  1:3859969/150‑1 (MQ=255)
                  cAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTcc                                                                                                                                     >  2:7416889/1‑150 (MQ=255)
                  cAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTcc                                                                                                                                     <  2:951931/150‑1 (MQ=255)
                      tACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATc                                                                                                                                 <  2:11855056/150‑1 (MQ=255)
                         gggCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACGTGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCgg                                                                                                                              <  2:11045313/150‑1 (MQ=255)
                           gcgcAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAg                                                                                                                                                 <  1:15298990/129‑1 (MQ=255)
                           gcgcAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGaa                                                                                                                            <  1:10023310/150‑1 (MQ=255)
                           gcgcAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGaa                                                                                                                            >  1:10364248/1‑150 (MQ=255)
                            cgcAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAAt                                                                                                                           >  2:9051154/1‑150 (MQ=255)
                                gTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCa                                                                                                                       <  1:753180/150‑1 (MQ=255)
                                   tggtggTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaa                                                                                                                    <  1:9292213/150‑1 (MQ=255)
                                      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                 <  2:1149544/150‑1 (MQ=255)
                                      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                 >  1:3730843/1‑150 (MQ=255)
                                      tggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaa                                                                                                                 >  1:4974493/1‑150 (MQ=255)
                                       ggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATaaaaaa                                                                                                                <  2:29614/150‑1 (MQ=255)
                                       ggtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGAGTCGTCACGGCGGCACGGTGGTATTTGTCGGCGTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATataaaa                                                                                                                <  2:1850709/150‑1 (MQ=255)
                                        gtTCGCTGGCGCAGAAAGTGATCGACGCGTCAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAg                                                                                                               >  2:4860765/1‑150 (MQ=255)
                                        gtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAg                                                                                                               >  2:13323031/1‑150 (MQ=255)
                                        gtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAg                                                                                                               >  2:427958/1‑150 (MQ=255)
                                        gtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAg                                                                                                               >  1:1387107/1‑150 (MQ=255)
                                        gtTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTAGTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAg                                                                                                               >  2:5623303/1‑150 (MQ=255)
                                           cGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGaaa                                                                                                            <  2:8577708/150‑1 (MQ=255)
                                                    aGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGACTTCCATAAAAAAGAAACGACgagga                                                                                                   >  2:10068559/1‑150 (MQ=255)
                                                      aaaGTTATCGACGCGACTGGTAATCAACTTGCGATTAATAACACGGTGATTTTGATTCGTCACGGCGGCAGGGTGGTTTTTGTGGGCCTTTTTAAAGGTTATTTGAAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGatgatg                                                                                                 <  2:8314029/150‑1 (MQ=255)
                                                        aGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATggg                                                                                               >  1:5090870/1‑150 (MQ=255)
                                                           gATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAg                                                                                            >  2:1264204/1‑150 (MQ=255)
                                                             tCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGcc                                                                                          >  1:7388423/1‑150 (MQ=255)
                                                             tCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGcc                                                                                          >  1:8180189/1‑150 (MQ=255)
                                                              cGACGCGACTGTTAATCAACATGCGATTAATAACACCGTTAATTTGATTCGTCACGGCGGCACGGTGTTATTTTTCGGCTTTTTTAAAGGTGTGTTGCAGTTCTCCGATCCGTAATTCAATAAAAAAGAAACGACGATGATGGGCAGCCg                                                                                         <  2:15494276/150‑1 (MQ=255)
                                                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                        >  1:4700032/1‑150 (MQ=255)
                                                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                        >  1:6339722/1‑150 (MQ=255)
                                                               gACGCGACAGGTAATCAACATGCGATGAATAACACCGAGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTAGCAGTTCTCCGATCCGGACTTCCATAAAAAAGAAACGACGATGATGGGCAGCCGc                                                                                        >  1:113858/1‑150 (MQ=255)
                                                                  gcgACTGTTAATCAACTTGCTTTTAATAACACCTTTATTTTTTTTCGTCTGGGGGGCAGGGTGTTTTTTGTGGGCTTTTTTAATGGTGAGTTGCTGTTCTCCGATCGGTATTTCAATAAAAAATAAACGACGATGATGGGCAGCCGCAac                                                                                     <  2:9248064/150‑1 (MQ=255)
                                                                  gcgACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAac                                                                                     >  2:7583492/1‑150 (MQ=255)
                                                                    gACTGGTAATCAACAGGCGATTATTAACACCGTGATTTTTATTCGTCACGGCGGCAGGGTGGTTTTTGTCGGCCTGTTTAAAGGTGAGTTGCTGTTCTCCGATCCGGATTTCCAAAAAAAAGAAACGACGATGATGGGCAGCCGCAacgc                                                                                   <  2:3658278/150‑1 (MQ=255)
                                                                       aGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccac                                                                                >  1:4224308/1‑150 (MQ=255)
                                                                         gTAATCACCCTGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCAGGGTGGTATTTGTGGGCCTGTTTAAAGGTGTGTTGCAGTTCTCGGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                              <  2:8803301/150‑1 (MQ=255)
                                                                         gTAATCAACTTGCTATGAATAACCCCGTGAATTTTATTGGTCACGGGGGCAGGGTGGTTTTTGGCGGCCTGTTTAATGGTGTTTTGCTGTTCTCCGATCCGGAATACCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                              <  2:336098/150‑1 (MQ=255)
                                                                         gTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAacgccacgc                                                                              >  2:14677023/1‑150 (MQ=255)
                                                                              cAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGaa                                                                         >  2:5245967/1‑150 (MQ=255)
                                                                                  aTGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGAtt                                                                     >  2:14955735/1‑150 (MQ=255)
                                                                                     cGATGAATAACACCGTGATTTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTTTTTAAAGGTGAGTTGCTGTTCTCCGATCCGTAATTCAATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTg                                                                  <  2:5156347/150‑1 (MQ=255)
                                                                                         taaaaaCACCGTGATTTTGATTCGTCACGGCGGCACGGTGTTATTTGTCGGCCTGTTTAATGGTGTGTTGCAGTTCGCCGATCCGGATTTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTaa                                                              <  2:11237436/146‑1 (MQ=255)
                                                                                            tAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGaa                                                                         <  2:12843747/136‑1 (MQ=255)
                                                                                            tAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGAAGCCGCAACGCCACGCCGGAAGAGTTTGCTAAAgt                                                           >  2:7253740/1‑150 (MQ=255)
                                                                                             aaCACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGCGAGTTGCAGTTATCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCTCAACGCCACGCCGGAAGATTTTGCTAAAttc                                                          >  2:2688900/1‑150 (MQ=255)
                                                                                                   gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                    >  1:12981496/1‑150 (MQ=255)
                                                                                                   gTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGa                                                    >  1:1874788/1‑150 (MQ=255)
                                                                                                    tGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGAc                                                   >  2:12211908/1‑150 (MQ=255)
                                                                                                    tGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGAc                                                   >  2:15222581/1‑150 (MQ=255)
                                                                                                        tttGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGAt                                               >  1:11217145/1‑150 (MQ=255)
                                                                                                        tttGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGAt                                               >  1:10590762/1‑150 (MQ=255)
                                                                                                               ggTCCCGGCGGCCGGGTGGTTTTTGTCGGCTTGTTTACTGGGGTGTTGCTGTTCGCCGACCCGGAATCCCAAAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaa                                        <  1:8412513/149‑1 (MQ=255)
                                                                                                                gTCACGGCGGCAGGGTGTTTTTTGTCGGCCTGTTTAATGGTGAGTTGCAGTTCTCCTATCCGTAATTCCATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaag                                       <  1:8759431/150‑1 (MQ=255)
                                                                                                                gTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaag                                       >  2:2584626/1‑150 (MQ=255)
                                                                                                                gTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaag                                       >  1:2764216/1‑150 (MQ=255)
                                                                                                                 tCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaagg                                      >  2:10601355/1‑150 (MQ=255)
                                                                                                                   aCGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCggaaggaa                                    >  1:2424552/1‑150 (MQ=255)
                                                                                                                   aCGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGAAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTTATGGCggaaggaa                                    >  2:6523365/1‑150 (MQ=255)
                                                                                                                      gcggcAGGGTGTTATTTGTCGGCCTTTTTAATGGTGAGTTGCAGTTCTCCGCCCCTTAATACAATAAAAAACACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGaaaaa                                 <  1:918856/150‑1 (MQ=255)
                                                                                                                         gccggttggtTTTTGTCGGCATGTTTAATGGTGAGTTGCAGTTCTCCGATCGGTATTTCCATAAAAAAGACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                              <  1:6587820/144‑1 (MQ=255)
                                                                                                                         gcAGGGTGGTTTTTGGGGGCCTTTTTAATGTTGTTTTTCATTTCACCTACCCGTAATACAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                              <  2:14432786/150‑1 (MQ=255)
                                                                                                                         gcACGGTGGTATTTGTCGGCCTGTTTAAAGGTGTGTTGCTGTCCTCCGAACCGGAATACCAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                              <  2:390974/150‑1 (MQ=255)
                                                                                                                         gcACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCa                              <  2:10364248/150‑1 (MQ=255)
                                                                                                                              gtggtTTTTGTGGGCCTTTTTAAGGGTGTTTTGCTGTTCTCCTATCCGTAATTCAAAAAAAAACAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActgct                         <  2:6895792/150‑1 (MQ=255)
                                                                                                                              gtggtATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCActgct                         >  2:8151637/1‑150 (MQ=255)
                                                                                                                                    tttGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatg                   >  1:6025812/1‑150 (MQ=255)
                                                                                                                                      tGTCGGCCTTTTTAATGGTGTGTTGCTGTTCTCCGTACCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACatgat                 <  1:10214128/150‑1 (MQ=255)
                                                                                                                                         cGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGtt              >  1:530366/1‑150 (MQ=255)
                                                                                                                                            cttttttAATGTTGTGTTGCTGTTCCCCGATGCGTATTCCAATAAAAAATACACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAc           <  1:7249296/146‑1 (MQ=255)
                                                                                                                                            ccTGTTTAATGTTGATTTCCTGTTCTCCTATCCGTAATTACATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAc           <  2:9660417/150‑1 (MQ=255)
                                                                                                                                            ccTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATCCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAc           <  2:1045897/150‑1 (MQ=255)
                                                                                                                                             cTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAAcc          >  2:9932589/1‑150 (MQ=255)
                                                                                                                                                 ttAACGGTTTGTTGCTGTTATCCTATCCTGAATTACATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATc      <  1:10981489/150‑1 (MQ=255)
                                                                                                                                                 ttAAAGGTTAGTTGCTGTTCCCCGATCCGTAATTCAATAAAAAATAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATc      <  2:13661261/150‑1 (MQ=255)
                                                                                                                                                  tAAAGGTGTGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCg     <  1:8929740/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCGAAGATTTTGACGCGCAGCTACGGGCGCAGTTTGGTGGTTCGCTGGCGCAGAAAGTGATCGACGCGACAGGTAATCAACATGCGATGAATAACACCGTGAATTTGATTCGTCACGGCGGCACGGTGGTATTTGTCGGCCTGTTTAAAGGTGAGTTGCAGTTCTCCGATCCGGAATTCCATAAAAAAGAAACGACGATGATGGGCAGCCGCAACGCCACGCCGGAAGATTTTGCTAAAGTCGGTCGACTGATGGCGGAAGGAAAAATCACTGCTGACATGATGTTAACCCATCGCTA  >  NC_000913/4596632‑4596930

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