Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,666,470 C→A 8.3% L386I (CTA→ATA)  clcB → H(+)/Cl(‑) exchange transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,666,4700CA 8.3% 628.2 / 30.1 175L386I (CTA→ATA) clcBH(+)/Cl(‑) exchange transporter
Reads supporting (aligned to +/- strand):  ref base C (98/62);  new base A (11/4);  total (109/66)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 176 total aligned reads displayed.

GTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGG  >  NC_000913/1666324‑1666619
                                                                                                                                                  |                                                                                                                                                      
gTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTAc                                                                                                                                                     >  2:2966557/1‑150 (MQ=255)
 tCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACt                                                                                                                                                    >  2:427073/1‑150 (MQ=255)
    tAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTccc                                                                                                                                                 >  2:521793/1‑150 (MQ=255)
      gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCg                                                                                                                                               <  2:1743381/150‑1 (MQ=255)
      gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCg                                                                                                                                               <  2:2448028/150‑1 (MQ=255)
      gCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCg                                                                                                                                               <  2:1863023/150‑1 (MQ=255)
        cTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑t                                                                                                                                            >  1:1489106/1‑150 (MQ=255)
          gggATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑ttt                                                                                                                                          >  2:1995991/1‑150 (MQ=255)
            gATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTat                                                                                                                                        >  1:951373/1‑150 (MQ=255)
                 tGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTattgatt                                                                                                                                   >  1:1893745/1‑150 (MQ=255)
                  ggTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTattgattg                                                                                                                                  <  1:331426/150‑1 (MQ=255)
                   gTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGc                                                                                                                                 >  1:2133679/1‑150 (MQ=255)
                    ttCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTCTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGcc                                                                                                                                >  2:470932/1‑150 (MQ=255)
                    ttCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGcc                                                                                                                                >  1:1140579/1‑150 (MQ=255)
                    ttCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTAGTGATTGcc                                                                                                                                >  2:2475853/1‑150 (MQ=255)
                       ccTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGc                                                                                                                             <  2:951373/150‑1 (MQ=255)
                         tGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGt                                                                                                                           >  1:2339564/1‑150 (MQ=255)
                             ggCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAAtt                                                                                                                       <  1:313738/150‑1 (MQ=255)
                               cGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGc                                                                                                                     >  2:1306380/1‑150 (MQ=255)
                               cGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGc                                                                                                                     <  2:652460/150‑1 (MQ=255)
                                gaagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCg                                                                                                                    >  1:2062899/1‑150 (MQ=255)
                                gaagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCg                                                                                                                    >  1:3218671/1‑150 (MQ=255)
                                gaagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCg                                                                                                                    >  1:3202503/1‑150 (MQ=255)
                                 aagaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGt                                                                                                                   >  2:1212374/1‑150 (MQ=255)
                                  agaaATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTc                                                                                                                  <  1:1377167/150‑1 (MQ=255)
                                                   tCGGATTGACCGGGATGGCGACACTGTAGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCACCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTAGCGGACGt                                                                                                 >  1:2214877/1‑150 (MQ=255)
                                                           aCCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGt                                                                                         >  2:1569065/1‑150 (MQ=255)
                                                            ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                                                                        >  1:1033813/1‑150 (MQ=255)
                                                            ccGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTg                                                                                        >  1:1033767/1‑150 (MQ=255)
                                                             cGGGATGGCGACCCTCTTGGCGGCGACCACGCCCGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACACCCCGG‑TTTCTTGAGTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                                                                       <  2:2952848/150‑1 (MQ=255)
                                                             cGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTCGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGa                                                                                       >  1:1358670/1‑150 (MQ=255)
                                                                gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActc                                                                                    <  1:1599438/150‑1 (MQ=255)
                                                                gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActc                                                                                    <  2:2732806/150‑1 (MQ=255)
                                                                gATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGCAATTTCGCGGACGTTACACCGTGActc                                                                                    <  1:3203647/150‑1 (MQ=255)
                                                                 aTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActct                                                                                   <  2:570003/150‑1 (MQ=255)
                                                                 aTGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGActct                                                                                   <  2:2622765/150‑1 (MQ=255)
                                                                   ggCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctat                                                                                 <  2:664030/150‑1 (MQ=255)
                                                                     cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                                                                               >  1:2514100/1‑150 (MQ=255)
                                                                     cGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatct                                                                               >  1:862202/1‑150 (MQ=255)
                                                                      gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatcta                                                                              >  1:1158091/1‑150 (MQ=255)
                                                                      gACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACtctatcta                                                                              >  1:1343383/1‑150 (MQ=255)
                                                                        cacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTAcc                                                                            <  1:682615/150‑1 (MQ=255)
                                                                        cacTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTAGCTGATACTCCCCg                                                                                                                                               <  2:1078410/84‑1 (MQ=39)
                                                                         acTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCg                                                                           <  2:3071455/150‑1 (MQ=255)
                                                                           tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                                                                         <  1:1497652/150‑1 (MQ=255)
                                                                           tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                                                                         <  2:1476604/150‑1 (MQ=255)
                                                                           tGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                                                                         <  1:2776315/150‑1 (MQ=255)
                                                                           tGTTGGCGGCGACCACGCACGCGACGATTATGTCGACGTTGATGATATGTGAATTGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGACTGCGTATTTGGGTCGGTTATTTCGCGGACGTTACACCGTGACTCTATCTACCGcc                                                                         <  1:3232334/150‑1 (MQ=255)
                                                                                  ggcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCacac                                                                  >  1:2558212/1‑150 (MQ=255)
                                                                                   gcgACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACt                                                                 >  1:1927457/1‑150 (MQ=255)
                                                                                     gACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgc                                                               >  2:1295883/1‑150 (MQ=255)
                                                                                      aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg                                                              <  2:2352577/150‑1 (MQ=255)
                                                                                      aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg                                                              >  2:1563596/1‑150 (MQ=255)
                                                                                      aCCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcg                                                              >  2:1141248/1‑150 (MQ=255)
                                                                                       ccACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTgcgc                                                             >  1:986069/1‑150 (MQ=255)
                                                                                         acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag                                                           <  2:372679/150‑1 (MQ=255)
                                                                                         acgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcag                                                           >  2:280719/1‑150 (MQ=255)
                                                                                          cgcacgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCgcagc                                                          <  2:2811331/150‑1 (MQ=255)
                                                                                             acgcGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATa                                                       >  1:1735885/1‑150 (MQ=255)
                                                                                                 gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTa                                                   >  2:905309/1‑150 (MQ=255)
                                                                                                 gcCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTa                                                   >  1:2957650/1‑150 (MQ=255)
                                                                                                 gcCGATTATGTCGACGTTGATGATAAGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑GTTATTGATTGCCAGCGTAATTCCGTCGGTAATTTCGCGGACGTTACCCCGTGACTCAATCGACCGCCAGCACACTGCGCAGCATAGCTa                                                   >  1:2341214/1‑150 (MQ=255)
                                                                                                  ccGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTaa                                                  <  1:757100/150‑1 (MQ=255)
                                                                                                    gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTCATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCCCACTGCGCAGCATAGCTAAGc                                                >  2:2877861/1‑150 (MQ=255)
                                                                                                    gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGc                                                >  1:91071/1‑150 (MQ=255)
                                                                                                    gATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGc                                                >  1:2283493/1‑150 (MQ=255)
                                                                                                      ttATGTCGACGTTGATGTTATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGACTGCGTAATTTCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATATACCGCCAGCACACTGCGCAGCATAGCTAAGCgg                                              >  1:294496/1‑150 (MQ=255)
                                                                                                      ttATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCgg                                              >  2:1607898/1‑150 (MQ=255)
                                                                                                        aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                            <  2:108352/150‑1 (MQ=255)
                                                                                                        aTGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGAt                                            <  2:310794/150‑1 (MQ=255)
                                                                                                          gTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTTCTCCCCGGTTTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATg                                           <  2:2849401/150‑1 (MQ=255)
                                                                                                           tCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTa                                         <  1:2633456/150‑1 (MQ=255)
                                                                                                             gACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTAtt                                       >  1:134727/1‑150 (MQ=255)
                                                                                                             gACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTAtt                                       <  1:2048560/150‑1 (MQ=255)
                                                                                                                  tgatgaTATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGcccc                                  >  1:1617117/1‑150 (MQ=255)
                                                                                                                           gtgAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGc                         <  2:1358670/150‑1 (MQ=255)
                                                                                                                            tgAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCt                        <  2:2213690/150‑1 (MQ=255)
                                                                                                                            tgAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCt                        <  2:3024131/150‑1 (MQ=255)
                                                                                                                               aaTGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGtt                     <  1:846806/150‑1 (MQ=255)
                                                                                                                                taaaCCGGGGAGTAGCTGATACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTc                    >  1:2987025/4‑150 (MQ=255)
                                                                                                                                taaaCCGGGGAGTAGCTGATACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTc                    >  1:52054/4‑150 (MQ=255)
                                                                                                                                aTGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTc                    <  1:673926/150‑1 (MQ=255)
                                                                                                                                 tGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCg                   >  1:1512126/1‑150 (MQ=255)
                                                                                                                                  aaCCGGGGAGTAGCTGATACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCgg                  >  1:2250691/2‑150 (MQ=255)
                                                                                                                                   aCCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCTGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgc                 <  2:31574/150‑1 (MQ=255)
                                                                                                                                   aCCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgc                 >  1:2406705/1‑150 (MQ=255)
                                                                                                                                    ccGGGGAGTAGCTGATACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcg                >  1:2987769/1‑150 (MQ=255)
                                                                                                                                    ccGGGGAGTAGCTGATACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcg                >  2:1799158/1‑150 (MQ=255)
                                                                                                                                     cGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCTGTTCACGCTGTTCGgcgc               <  1:2164107/150‑1 (MQ=255)
                                                                                                                                     cGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgc               >  1:330841/1‑150 (MQ=255)
                                                                                                                                     cGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgc               <  1:2641831/150‑1 (MQ=255)
                                                                                                                                     cGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgc               >  1:2165082/1‑150 (MQ=255)
                                                                                                                                     cGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgc               >  1:2969142/1‑150 (MQ=255)
                                                                                                                                      ggggAGTAGCTGATACTCCCCGG‑TTTATTGATTTCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgcg              >  1:571986/1‑150 (MQ=255)
                                                                                                                                       gggAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGgcgcgc             >  2:941887/1‑150 (MQ=255)
                                                                                                                                              tCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGAt      >  1:2924242/1‑150 (MQ=255)
                                                                                                                                              tCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGAt      >  1:2924241/1‑150 (MQ=255)
                                                                                                                                                aGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATaa    >  1:1533457/1‑150 (MQ=255)
                                                                                                                                                 gCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAg   >  1:647142/1‑150 (MQ=255)
                                                                                                                                                  cTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAgg  >  2:3063379/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
GTCGTAGCCTGGGATTATGGTTCCCTGATGGCGAAGAAATTACACTTTTACTCGGATTGACCGGGATGGCGACACTGTTGGCGGCGACCACGCACGCGCCGATTATGTCGACGTTGATGATATGTGAAATGACCGGGGAGTATCAGCTACTCCCCGG‑TTTATTGATTGCCTGCGTAATTGCGTCGGTAATTTCGCGGACGTTACACCGTGACTCTATCTACCGCCAGCACACTGCGCAGCATAGCTAAGCGGATGTATTGCCCCAGTTCACGCTGTTCGGCGCGCGGCAGATAAGG  >  NC_000913/1666324‑1666619

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