Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,074,695 A→G intergenic (+37/‑84) flu → / → yeeR CP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,074,6950AG99.2% 450.7 / ‑4.4 119intergenic (+37/‑84)flu/yeeRCP4‑44 prophage; antigen 43 (Ag43) phase‑variable biofilm formation autotransporter/CP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (60/58);  total (61/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.93e-01

Only 100 of 119 total aligned reads displayed.

CCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTT  >  NC_000913/2074547‑2074840
                                                                                                                                                    |                                                                                                                                                 
ccGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGcc                                                                                                                                                 >  1:40820/1‑151 (MQ=255)
ccGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGcc                                                                                                                                                 >  1:210915/1‑151 (MQ=255)
ccGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGcc                                                                                                                                                 >  1:1951094/1‑151 (MQ=255)
  gtgtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCgg                                                                                                                                               <  1:344609/151‑1 (MQ=255)
  gtgtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCgg                                                                                                                                               >  1:2056538/1‑151 (MQ=255)
    gtCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCggg                                                                                                                                              >  2:634325/1‑150 (MQ=255)
      ccGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGAc                                                                                                                                            <  2:907718/150‑1 (MQ=255)
        gggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCg                                                                                                                                         >  2:1622266/1‑151 (MQ=255)
        gggAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCg                                                                                                                                         >  1:1322553/1‑151 (MQ=255)
        gggAAAATACCCCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCg                                                                                                                                         <  1:114311/151‑1 (MQ=255)
            aaaTATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGAcc                                                                                                                                     >  1:1541667/1‑151 (MQ=255)
                atCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGgcg                                                                                                                                 <  2:1680011/151‑1 (MQ=255)
                atCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGgc                                                                                                                                  <  1:713955/150‑1 (MQ=255)
                   aCCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCa                                                                                                                               >  1:1385861/1‑150 (MQ=255)
                      cTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACgg                                                                                                                           >  2:2128961/1‑151 (MQ=255)
                       tGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACgg                                                                                                                           >  1:1663011/1‑150 (MQ=255)
                        gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGa                                                                                                                          <  2:673464/150‑1 (MQ=255)
                        gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGAt                                                                                                                         >  2:1966010/1‑151 (MQ=255)
                        gggCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGAt                                                                                                                         <  2:85960/151‑1 (MQ=255)
                          gCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCt                                                                                                                       <  2:418944/151‑1 (MQ=255)
                           cGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCtt                                                                                                                      >  1:2009001/1‑151 (MQ=255)
                            gTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTc                                                                                                                     >  2:879837/1‑151 (MQ=255)
                            gTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTc                                                                                                                     >  1:1663457/1‑151 (MQ=255)
                            gTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTc                                                                                                                     >  1:1458479/1‑151 (MQ=255)
                             ttCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTc                                                                                                                     >  2:584515/1‑150 (MQ=255)
                              tCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCa                                                                                                                    <  1:1875816/150‑1 (MQ=255)
                               cAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCaa                                                                                                                   <  1:278539/150‑1 (MQ=255)
                                 ggCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGc                                                                                                                <  1:404239/151‑1 (MQ=255)
                                    cGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCcac                                                                                                             >  2:1031437/1‑151 (MQ=255)
                                      gTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACAt                                                                                                           >  1:569633/1‑151 (MQ=255)
                                        tATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACAtt                                                                                                          <  2:1281297/150‑1 (MQ=255)
                                         aTGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                        <  2:190870/151‑1 (MQ=255)
                                         aTGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCg                                                                                                        >  1:1159848/1‑151 (MQ=255)
                                           gCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGc                                                                                                       <  1:43237/150‑1 (MQ=255)
                                           gCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCt                                                                                                      <  2:1253429/151‑1 (MQ=255)
                                             ccACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTgg                                                                                                    <  1:409573/151‑1 (MQ=255)
                                              cacaGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGc                                                                                                   <  2:1439535/151‑1 (MQ=255)
                                                caGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAt                                                                                                 <  2:758845/151‑1 (MQ=255)
                                                 aGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCAtt                                                                                                >  2:1358914/1‑151 (MQ=255)
                                                  gCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTa                                                                                               <  2:1849319/151‑1 (MQ=255)
                                                   cGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTaa                                                                                              <  1:1823034/151‑1 (MQ=255)
                                                      cAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACa                                                                                            >  1:944492/1‑150 (MQ=255)
                                                          ggCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAAc                                                                                       >  2:1050705/1‑151 (MQ=255)
                                                           gcagcaGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAAc                                                                                       >  2:1239859/1‑150 (MQ=255)
                                                           gcagcaGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACa                                                                                      <  1:1419300/151‑1 (MQ=255)
                                                            cagcagCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACAt                                                                                     >  2:2155264/1‑151 (MQ=255)
                                                              gcagcGCTGAAGGGTATAATGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGa                                                                                   <  1:973699/151‑1 (MQ=255)
                                                               cagcGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGat                                                                                  <  2:601618/151‑1 (MQ=255)
                                                                  cgcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATAtt                                                                               <  1:1241434/151‑1 (MQ=255)
                                                                   gcTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtc                                                                              <  1:1139225/151‑1 (MQ=255)
                                                                     tGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcat                                                                            >  2:2043042/1‑151 (MQ=255)
                                                                      gAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATtcatc                                                                           <  1:921383/151‑1 (MQ=255)
                                                                          ggTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACgg                                                                       >  2:665306/1‑151 (MQ=255)
                                                                             ataACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGt                                                                    <  2:148636/151‑1 (MQ=255)
                                                                             ataACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGt                                                                    <  2:2124384/151‑1 (MQ=255)
                                                                               aaCGGTCAGGCCAGACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTg                                                                   >  2:168020/1‑150 (MQ=255)
                                                                               aaCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTg                                                                   >  2:176898/1‑150 (MQ=255)
                                                                               aaCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGa                                                                  <  2:1109201/151‑1 (MQ=255)
                                                                                aCGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGa                                                                  <  1:405399/150‑1 (MQ=255)
                                                                                    tCAGGCCACACTGAATGTGACCTTCTGACAGAACAATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTAt                                                              >  1:924851/1‑150 (MQ=255)
                                                                                      aGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGtt                                                           <  2:50453/151‑1 (MQ=255)
                                                                                      aGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGtt                                                           >  2:368416/1‑151 (MQ=255)
                                                                                       ggCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGtt                                                           >  2:1615377/1‑150 (MQ=255)
                                                                                        gCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                         <  2:1172600/151‑1 (MQ=255)
                                                                                        gCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                         <  2:2113977/151‑1 (MQ=255)
                                                                                        gCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                         <  1:128955/151‑1 (MQ=255)
                                                                                         ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACa                                                        <  1:1982839/151‑1 (MQ=255)
                                                                                         ccACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACa                                                        <  1:1426111/151‑1 (MQ=255)
                                                                                          cacaCTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTAc                                                         >  2:256989/1‑149 (MQ=255)
                                                                                            cacTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGAt                                                     <  2:318027/151‑1 (MQ=255)
                                                                                              cTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATa                                                    <  1:1571824/150‑1 (MQ=255)
                                                                                                      gACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcg                                           <  2:210914/151‑1 (MQ=255)
                                                                                                      gACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcg                                           >  2:2169820/1‑151 (MQ=255)
                                                                                                       aCCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGgcgcgc                                          <  2:1015011/151‑1 (MQ=255)
                                                                                                        ccTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCt                                         <  2:1146762/151‑1 (MQ=255)
                                                                                                         cTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTg                                        >  2:198275/1‑151 (MQ=255)
                                                                                                         cTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTg                                        >  1:1720781/1‑151 (MQ=255)
                                                                                                           tCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGAt                                      >  1:1380203/1‑151 (MQ=255)
                                                                                                            cTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATc                                     >  2:2020598/1‑151 (MQ=255)
                                                                                                               aCAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctg                                  >  2:654860/1‑151 (MQ=255)
                                                                                                                   aaCCATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCctgctg                               >  2:741000/1‑150 (MQ=255)
                                                                                                                     ccATCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATc                            <  2:1835427/151‑1 (MQ=255)
                                                                                                                        tCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGc                          >  1:55106/1‑150 (MQ=255)
                                                                                                                        tCGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGc                          >  2:1649573/1‑150 (MQ=255)
                                                                                                                         cGCCTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCa                         <  2:724147/150‑1 (MQ=255)
                                                                                                                           ccTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAgg                       >  1:1187025/1‑150 (MQ=255)
                                                                                                                           ccTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAgg                       >  2:1802274/1‑150 (MQ=255)
                                                                                                                           ccTCTCTGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGa                      <  2:1047643/151‑1 (MQ=255)
                                                                                                                             tctctGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGAtt                    >  2:1781826/1‑151 (MQ=255)
                                                                                                                                ctGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGAttt                   <  1:589740/149‑1 (MQ=255)
                                                                                                                                ctGTGGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGc                 >  1:298162/1‑151 (MQ=255)
                                                                                                                                 tgtgGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGcc                >  2:1503962/1‑151 (MQ=255)
                                                                                                                                 tgtgGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGcc                >  2:84250/1‑151 (MQ=255)
                                                                                                                                  gtgGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTg                  <  1:489366/148‑1 (MQ=255)
                                                                                                                                   tgGTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAt              >  2:2136969/1‑151 (MQ=255)
                                                                                                                                    ggTCCCGGTCATCATGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCAtt             >  2:807612/1‑151 (MQ=255)
                                                                                                                                            tcatcaTGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCt     >  1:486683/1‑151 (MQ=255)
                                                                                                                                             catcatGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCtt    <  1:900587/151‑1 (MQ=255)
                                                                                                                                             catcatGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCt     <  2:1580254/150‑1 (MQ=255)
                                                                                                                                               tcatGGCCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCtttt  >  2:134631/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CCGTGTCCGGGAAAATATCACCCTGGGCGTTCAGGCCGGTTATGCCCACAGCGTCAGCGGCAGCAGCGCTGAAGGGTATAACGGTCAGGCCACACTGAATGTGACCTTCTGACAGAACCATCGCCTCTCTGTGGTCCCGGTCATCATGACCGGGACCCGGACCGGCGCAACGGATCTTCAACGCCACATTCGCTGGCATTAACAATAACATGATATTCATCACGGAGTGACTATGTTACAGATAGTCGGCGCGCTGATCCTGCTGATCGCAGGATTTGCCATTCTTCGCCTTTT  >  NC_000913/2074547‑2074840

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: