Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,323 G→C 10.6% R180G (CGG→GGG)  lacA ← thiogalactoside acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,3230GC10.6% 1608.2 / 31.0 513R180G (CGG→GGG) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (135/313);  new base C (59/4);  total (196/317)
Fisher's exact test for biased strand distribution p-value = 3.79e-23
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 520 total aligned reads displayed.

ATGCCGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAA  >  NC_000913/361174‑361472
                                                                                                                                                      |                                                                                                                                                     
   ccGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACa                                                                                                                                                    >  2:6372770/1‑150 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACaa                                                                                                                                                   <  2:4396888/150‑1 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACaa                                                                                                                                                   <  1:5204373/150‑1 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACaa                                                                                                                                                   <  1:5561124/150‑1 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑AGAACCCGACaa                                                                                                                                                   <  2:361156/150‑1 (MQ=255)
     ggATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAg                                                                                                                                                  <  1:3341092/150‑1 (MQ=255)
        tGCGGCTAATGTAGATCGCTGATCTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGaa                                                                                                                                               <  1:6602591/150‑1 (MQ=255)
        tGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGaa                                                                                                                                               <  2:5711042/150‑1 (MQ=255)
         gCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAAc                                                                                                                                              <  1:7350805/150‑1 (MQ=255)
          cGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACg                                                                                                                                             <  1:4050486/150‑1 (MQ=255)
           ggCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGc                                                                                                                                            <  1:6236381/150‑1 (MQ=255)
           ggCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGc                                                                                                                                            <  1:10909550/150‑1 (MQ=255)
            gCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTTTTTGGCGT‑AGACACCGACAATGAACGcc                                                                                                                                           >  2:5923088/1‑150 (MQ=38)
             cTAATGTAGATCGCTGAACTTGTAGGCCTTATAAGCGCAGCGTATCATGAAATTTTTATAATTTAAACTGACGATTAAACTTTATAATCTTTTAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCa                                                                                                                                          <  1:9544353/150‑1 (MQ=255)
                aTGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAgcc                                                                                                                                       <  2:5305872/150‑1 (MQ=255)
                aTGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAgcc                                                                                                                                       <  2:5126273/150‑1 (MQ=255)
                aTGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGA‑ATAACCCGACAAGGAACGCCAgcc                                                                                                                                       >  1:10078164/1‑150 (MQ=255)
                   tAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGA‑ATAACCCCACAAGGAACGCCAgccgcc                                                                                                                                    >  1:10434915/1‑150 (MQ=255)
                         gCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTAGCGA‑AAAAACCCACCAGGAACGCCAcccccccccacc                                                                                                                              >  1:8466629/1‑144 (MQ=38)
                               cTTGTTGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAAATTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTtggt                                                                                                                        <  1:8203785/150‑1 (MQ=255)
                               cTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTTAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTtggt                                                                                                                        <  1:8877629/150‑1 (MQ=255)
                                 tGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTtggtgg                                                                                                                      <  2:9343043/150‑1 (MQ=255)
                                  gTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGa                                                                                                                     <  1:9954267/150‑1 (MQ=255)
                                  gTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGa                                                                                                                     <  1:8772748/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATACCCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  1:4587804/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  2:3986136/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  2:9032543/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  2:1802705/150‑1 (MQ=255)
                                      gCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTAGCGA‑ATAAAACGACAATTAAAGCCACCCCCCACGACGTTTTTTGGAAttt                                                                                                                 >  1:2462783/1‑148 (MQ=255)
                                        cTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCt                                                                                                               <  1:10977178/150‑1 (MQ=255)
                                         ttaTAAGCGCAGCGTATAAGTAAATTTTTATAATTTAAACTGACTATTCAACTTTATAATCTTTTAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCCACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCtt                                                                                                              <  1:2981865/148‑1 (MQ=255)
                                         tGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTTTTTGGGGG‑ATAACCCGACAATGGACGCCAGCCGCCACGACGGTTTTTTGAATGTGtt                                                                                                              >  1:10225181/1‑150 (MQ=255)
                                         tGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAAAGCCAGCCGCCACGCCGTTTGGTGGAATGTCtt                                                                                                              >  2:666223/1‑150 (MQ=255)
                                         tGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCCCCACGACGGGTGGTGGGATGTCtt                                                                                                              >  2:539323/1‑150 (MQ=255)
                                           aTAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑AAAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCtttt                                                                                                            >  2:1131467/1‑150 (MQ=255)
                                            tAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGA‑ATAAACCCACAAGGAAACCCAGCCGCCAACACCTTTTGGGGAAAGTCttttt                                                                                                           >  2:1320246/1‑149 (MQ=255)
                                              aGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTtgtg                                                                                                         <  1:9696367/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGTCAATTTTAATAATTTAAACTTACCATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGAAAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  1:2790101/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTTTTTGGGGT‑ATAAACCGACAAGGTAAGCCAGCCGCCAAGACGTTTGGTGGAATTTGtttttttt                                                                                                        >  1:1704215/1‑146 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATACCCCGCCAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  2:3084571/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  2:8787475/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  1:7610823/150‑1 (MQ=255)
                                                cgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGAc                                                                                                       <  2:6870159/150‑1 (MQ=255)
                                                  cAGCTTATCAGGCAATTTTTAAAAATTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCTA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGa                                                                                                     <  1:5405557/150‑1 (MQ=255)
                                                     cGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGACATAACCC‑ACCAAGAAAGCCACCCCCCCCCACCATTGGTGGGATTTATTTTGTGTCGAtac                                                                                                  >  1:1891057/1‑150 (MQ=255)
                                                       tATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtacta                                                                                                <  2:948482/150‑1 (MQ=255)
                                                        aTCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtactac                                                                                               <  1:1275243/150‑1 (MQ=255)
                                                        aTCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtactac                                                                                               <  2:6149596/150‑1 (MQ=255)
                                                        aTCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtactac                                                                                               <  2:10816459/150‑1 (MQ=255)
                                                           aGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑AAAAACCCACAACGCAACCCAGCCCCCACCACCTTTGGTGGAATGTCATTTTTGAAGAtacaaacac                                                                                            >  1:1184625/1‑144 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTGGGGA‑CTAAACCCCCAACGAAAGCCAGCCCCCCCCACCTTTGGTGGGATGTAGTTTGTTACGATAATACCcgcc                                                                                          >  2:2886990/1‑149 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGA‑ATAACCCCACAAGGAACACCAGCCCCCACCACGGTTTGTGGAATGTCTTTTGTGACGATAAAAACcgcg                                                                                          >  1:5844541/1‑150 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGA‑ATAAACCGACCAGGAACGCCAGCCGCCACCACGGTTTGTGGAATGTCTTTTTTGACGATACTACCcccg                                                                                          >  1:6937853/1‑150 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTGGCGA‑ATAACCCGACAACGAACGCCAGCCGCCCCGCCCTTTGGGGGCATGTCTTTTGTGAAGATAAAAcccccc                                                                                          >  2:4918901/1‑147 (MQ=255)
                                                                 tttttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTATCGA‑ATAAACCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCGTTTGTGACGATACTACCCGCGCCaa                                                                                      >  2:8095976/1‑150 (MQ=255)
                                                                   tttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATa                                                                                    <  1:4957669/150‑1 (MQ=255)
                                                                    ttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATaa                                                                                   <  1:4402856/150‑1 (MQ=255)
                                                                    ttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGGCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATaa                                                                                   <  1:379251/150‑1 (MQ=255)
                                                                      ataATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACa                                                                                 <  1:5702362/150‑1 (MQ=255)
                                                                      ataATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACa                                                                                 <  2:4514744/150‑1 (MQ=255)
                                                                         aTTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGaa                                                                              <  1:6823779/150‑1 (MQ=255)
                                                                         aTTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGaa                                                                              <  2:7625528/150‑1 (MQ=255)
                                                                          tttAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAAACCCACAAGGAACGCCAGCCGCCACCACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAAt                                                                             >  2:648702/1‑150 (MQ=255)
                                                                            tAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTa                                                                           <  2:7118040/150‑1 (MQ=255)
                                                                             aaaCTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTAt                                                                          <  2:2559988/150‑1 (MQ=255)
                                                                              aaCTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATc                                                                         <  1:1471162/150‑1 (MQ=255)
                                                                                 tGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAAACCGACAAGGAACCCCAGCCCCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACAACCCGCGCCAATAACAGAAATATaccc                                                                      >  2:7893694/1‑150 (MQ=255)
                                                                                       aacAACTTTATAATCTATGAAATAATAGTGCTTATCCCGGTCGATTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGATTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGt                                                                <  2:1379902/148‑1 (MQ=255)
                                                                                        tCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTg                                                               >  2:5492847/1‑150 (MQ=255)
                                                                                             tttATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGA‑AAAAACCCACAAGGAACGCCAGCCCCCACGACGTTTGGTGGAAAGTATTTTGTGACGATACAAACCGCGCCAAAAACAGAAATATACCCCATTGTGACGcc                                                          >  2:2339589/1‑150 (MQ=255)
                                                                                               tataATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGC‑CTAACCCGACAAGGAACGCCCGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATACCAGAATTATCCCCGATGGTGACGCCTg                                                        <  2:10129036/150‑1 (MQ=255)
                                                                                                 taATCTTTGAAATAATAGTGCTGATCCCGGTTGTTTATTTCGCGC‑ATCACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACGACCCGCGCCAATTACAGAATTATCCCCGATGGTGACGCCTGGa                                                      <  1:10837873/150‑1 (MQ=255)
                                                                                                      tttGTAATAATAGTGCTTATCCCGGTCGTTTATTTCGAGA‑ATAACCAGACTAGGAACGCCAGCCGCCACTTCGTTTGGTGGAATGTCTTTTGTGACGATACTACACGCGCCAATATAAGAATTATCCCCGATGGTGACGCCTGGATtaat                                                 <  2:733850/150‑1 (MQ=255)
                                                                                                      tttGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAAACCCACAACGAACGCCAGCCGCCACCACGTTTGGTGGAATGTCTTTTGTGACGATACAACCCGCGCCAATAAAACAAATATACCCCATGGTGACGCCGGGATtaaa                                                 >  2:3795266/1‑149 (MQ=255)
                                                                                                      tttGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGA‑ATAAACCGACAAAGAAACCCAGCCCCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTAACCCCGCCAAAAACAGAATTATACCCGATGGTGAAGCCTGGATtaat                                                 >  2:1095407/1‑150 (MQ=255)
                                                                                                        tGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATtaataa                                               <  1:2256844/150‑1 (MQ=255)
                                                                                                         gAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAAc                                              <  2:10740597/150‑1 (MQ=255)
                                                                                                              aataGTGCTTATCCCGGTCGTTTATTTAGCGA‑AAAAACCCACAAGGAACGCCAGCCGCCACCCCGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAAATATCCCCCATGGTGACGCCTGGATTAATAACcacaa                                         >  2:3354108/1‑149 (MQ=255)
                                                                                                                  gTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAAACCGACAAGGAAACCCAACCCCCACGACGTTTGGTGGGATGTCTTTTTTTACGATACAAACCCCGCCAAAAACAGAATTATACCCGATGGTGACGCCGGGGTTAATAACcacaagtct                                     >  2:1964321/1‑145 (MQ=255)
                                                                                                                   tGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACtt                                    <  1:3423981/150‑1 (MQ=255)
                                                                                                                     cTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTcc                                  <  1:309299/150‑1 (MQ=255)
                                                                                                                      ttATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACCCCAGCCGCCACGCCGGTTGGTGGAAAGTCTTTTGTGACGATACAACCCCCGCCAAAAACACAAATATACCCCATGGTGACGCCAGGATTAATAAACACATGACATCCg                                 >  2:1504909/1‑150 (MQ=255)
                                                                                                                        aTCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGTATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTTATAACCACATGACTTCCGAt                               <  2:9742135/150‑1 (MQ=255)
                                                                                                                           ccGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCa                            <  1:9186124/150‑1 (MQ=255)
                                                                                                                           ccGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCa                            <  2:11095194/150‑1 (MQ=255)
                                                                                                                                    tATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTAtt                   <  2:680744/150‑1 (MQ=255)
                                                                                                                                     aTTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTg                  <  2:10420756/150‑1 (MQ=255)
                                                                                                                                        tCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCa               >  1:7203540/1‑150 (MQ=255)
                                                                                                                                        tCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCa               >  2:10328836/1‑150 (MQ=255)
                                                                                                                                         cgcgA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCaa              <  2:6954226/150‑1 (MQ=255)
                                                                                                                                          gcgA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAAt             >  1:11030405/1‑150 (MQ=255)
                                                                                                                                           cgA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATc            >  2:2120635/1‑150 (MQ=255)
                                                                                                                                           cgA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATc            <  2:6312384/150‑1 (MQ=255)
                                                                                                                                            gA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCg           <  1:8579985/150‑1 (MQ=255)
                                                                                                                                            gA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCg           <  1:6470124/150‑1 (MQ=255)
                                                                                                                                                tAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCCTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTa        <  1:4221237/150‑1 (MQ=255)
                                                                                                                                                  aCCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATc      <  1:2099836/150‑1 (MQ=255)
                                                                                                                                                  aCCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATc      >  1:7074714/1‑150 (MQ=255)
                                                                                                                                                   cccGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCg     <  1:9438366/150‑1 (MQ=255)
                                                                                                                                                    ccGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCgg    <  2:8597240/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATGCCGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGA‑ATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAA  >  NC_000913/361174‑361472

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