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.

ATGCCGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAA  >  NC_000913/361174‑361472
                                                                                                                                                      |                                                                                                                                                     
   ccGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACa                                                                                                                                                    >  2:6372770/1‑150 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACaa                                                                                                                                                   <  2:4396888/150‑1 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACaa                                                                                                                                                   <  1:5204373/150‑1 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACaa                                                                                                                                                   <  1:5561124/150‑1 (MQ=255)
    cGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAAGAACCCGACaa                                                                                                                                                   <  2:361156/150‑1 (MQ=255)
     ggATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAg                                                                                                                                                  <  1:3341092/150‑1 (MQ=255)
        tGCGGCTAATGTAGATCGCTGATCTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGaa                                                                                                                                               <  1:6602591/150‑1 (MQ=255)
        tGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGaa                                                                                                                                               <  2:5711042/150‑1 (MQ=255)
         gCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAAc                                                                                                                                              <  1:7350805/150‑1 (MQ=255)
          cGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACg                                                                                                                                             <  1:4050486/150‑1 (MQ=255)
           ggCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGc                                                                                                                                            <  1:6236381/150‑1 (MQ=255)
           ggCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGc                                                                                                                                            <  1:10909550/150‑1 (MQ=255)
            gCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTTTTTGGCGTAGACACCGACAATGAACGcc                                                                                                                                           >  2:5923088/1‑150 (MQ=38)
             cTAATGTAGATCGCTGAACTTGTAGGCCTTATAAGCGCAGCGTATCATGAAATTTTTATAATTTAAACTGACGATTAAACTTTATAATCTTTTAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCa                                                                                                                                          <  1:9544353/150‑1 (MQ=255)
                aTGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAgcc                                                                                                                                       <  2:5305872/150‑1 (MQ=255)
                aTGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAgcc                                                                                                                                       <  2:5126273/150‑1 (MQ=255)
                aTGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGAATAACCCGACAAGGAACGCCAgcc                                                                                                                                       >  1:10078164/1‑150 (MQ=255)
                   tAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGAATAACCCCACAAGGAACGCCAgccgcc                                                                                                                                    >  1:10434915/1‑150 (MQ=255)
                         gCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTAGCGAAAAAACCCACCAGGAACGCCAcccccccccacc                                                                                                                              >  1:8466629/1‑144 (MQ=38)
                               cTTGTTGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAAATTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTtggt                                                                                                                        <  1:8203785/150‑1 (MQ=255)
                               cTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTTAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTtggt                                                                                                                        <  1:8877629/150‑1 (MQ=255)
                                 tGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTtggtgg                                                                                                                      <  2:9343043/150‑1 (MQ=255)
                                  gTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGa                                                                                                                     <  1:9954267/150‑1 (MQ=255)
                                  gTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGa                                                                                                                     <  1:8772748/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATACCCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  1:4587804/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  2:3986136/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  2:9032543/150‑1 (MQ=255)
                                   tAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGaa                                                                                                                    <  2:1802705/150‑1 (MQ=255)
                                      gCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTAGCGAATAAAACGACAATTAAAGCCACCCCCCACGACGTTTTTTGGAAttt                                                                                                                 >  1:2462783/1‑148 (MQ=255)
                                        cTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCt                                                                                                               <  1:10977178/150‑1 (MQ=255)
                                         ttaTAAGCGCAGCGTATAAGTAAATTTTTATAATTTAAACTGACTATTCAACTTTATAATCTTTTAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCCACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCtt                                                                                                              <  1:2981865/148‑1 (MQ=255)
                                         tGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTTTTTGGGGGATAACCCGACAATGGACGCCAGCCGCCACGACGGTTTTTTGAATGTGtt                                                                                                              >  1:10225181/1‑150 (MQ=255)
                                         tGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAAAGCCAGCCGCCACGCCGTTTGGTGGAATGTCtt                                                                                                              >  2:666223/1‑150 (MQ=255)
                                         tGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTCGCGAATAACCCGACAAGGAACGCCAGCCCCCACGACGGGTGGTGGGATGTCtt                                                                                                              >  2:539323/1‑150 (MQ=255)
                                           aTAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAAAAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCtttt                                                                                                            >  2:1131467/1‑150 (MQ=255)
                                            tAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGAATAAACCCACAAGGAAACCCAGCCGCCAACACCTTTTGGGGAAAGTCttttt                                                                                                           >  2:1320246/1‑149 (MQ=255)
                                              aGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTtgtg                                                                                                         <  1:9696367/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGTCAATTTTAATAATTTAAACTTACCATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGAAAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  1:2790101/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTTTTTGGGGTATAAACCGACAAGGTAAGCCAGCCGCCAAGACGTTTGGTGGAATTTGtttttttt                                                                                                        >  1:1704215/1‑146 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATACCCCGCCAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  2:3084571/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  2:8787475/150‑1 (MQ=255)
                                               gcgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGa                                                                                                        <  1:7610823/150‑1 (MQ=255)
                                                cgcAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGAc                                                                                                       <  2:6870159/150‑1 (MQ=255)
                                                  cAGCTTATCAGGCAATTTTTAAAAATTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCTAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGa                                                                                                     <  1:5405557/150‑1 (MQ=255)
                                                     cGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGACATAACCCACCAAGAAAGCCACCCCCCCCCACCATTGGTGGGATTTATTTTGTGTCGAtac                                                                                                  >  1:1891057/1‑150 (MQ=255)
                                                       tATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtacta                                                                                                <  2:948482/150‑1 (MQ=255)
                                                        aTCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtactac                                                                                               <  1:1275243/150‑1 (MQ=255)
                                                        aTCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtactac                                                                                               <  2:6149596/150‑1 (MQ=255)
                                                        aTCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGAtactac                                                                                               <  2:10816459/150‑1 (MQ=255)
                                                           aGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAAAAAACCCACAACGCAACCCAGCCCCCACCACCTTTGGTGGAATGTCATTTTTGAAGAtacaaacac                                                                                            >  1:1184625/1‑144 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTGGGGACTAAACCCCCAACGAAAGCCAGCCCCCCCCACCTTTGGTGGGATGTAGTTTGTTACGATAATACCcgcc                                                                                          >  2:2886990/1‑149 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGAATAACCCCACAAGGAACACCAGCCCCCACCACGGTTTGTGGAATGTCTTTTGTGACGATAAAAACcgcg                                                                                          >  1:5844541/1‑150 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGAATAAACCGACCAGGAACGCCAGCCGCCACCACGGTTTGTGGAATGTCTTTTTTGACGATACTACCcccg                                                                                          >  1:6937853/1‑150 (MQ=255)
                                                             gCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTGTATTTGGCGAATAACCCGACAACGAACGCCAGCCGCCCCGCCCTTTGGGGGCATGTCTTTTGTGAAGATAAAAcccccc                                                                                          >  2:4918901/1‑147 (MQ=255)
                                                                 tttttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTATCGAATAAACCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCGTTTGTGACGATACTACCCGCGCCaa                                                                                      >  2:8095976/1‑150 (MQ=255)
                                                                   tttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATa                                                                                    <  1:4957669/150‑1 (MQ=255)
                                                                    ttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATaa                                                                                   <  1:4402856/150‑1 (MQ=255)
                                                                    ttATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGGCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATaa                                                                                   <  1:379251/150‑1 (MQ=255)
                                                                      ataATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACa                                                                                 <  1:5702362/150‑1 (MQ=255)
                                                                      ataATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACa                                                                                 <  2:4514744/150‑1 (MQ=255)
                                                                         aTTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGaa                                                                              <  1:6823779/150‑1 (MQ=255)
                                                                         aTTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGaa                                                                              <  2:7625528/150‑1 (MQ=255)
                                                                          tttAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAAACCCACAAGGAACGCCAGCCGCCACCACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAAt                                                                             >  2:648702/1‑150 (MQ=255)
                                                                            tAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTa                                                                           <  2:7118040/150‑1 (MQ=255)
                                                                             aaaCTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTAt                                                                          <  2:2559988/150‑1 (MQ=255)
                                                                              aaCTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATc                                                                         <  1:1471162/150‑1 (MQ=255)
                                                                                 tGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAAACCGACAAGGAACCCCAGCCCCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACAACCCGCGCCAATAACAGAAATATaccc                                                                      >  2:7893694/1‑150 (MQ=255)
                                                                                       aacAACTTTATAATCTATGAAATAATAGTGCTTATCCCGGTCGATTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGATTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGt                                                                <  2:1379902/148‑1 (MQ=255)
                                                                                        tCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTg                                                               >  2:5492847/1‑150 (MQ=255)
                                                                                             tttATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGAAAAAACCCACAAGGAACGCCAGCCCCCACGACGTTTGGTGGAAAGTATTTTGTGACGATACAAACCGCGCCAAAAACAGAAATATACCCCATTGTGACGcc                                                          >  2:2339589/1‑150 (MQ=255)
                                                                                               tataATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGCCTAACCCGACAAGGAACGCCCGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATACCAGAATTATCCCCGATGGTGACGCCTg                                                        <  2:10129036/150‑1 (MQ=255)
                                                                                                 taATCTTTGAAATAATAGTGCTGATCCCGGTTGTTTATTTCGCGCATCACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACGACCCGCGCCAATTACAGAATTATCCCCGATGGTGACGCCTGGa                                                      <  1:10837873/150‑1 (MQ=255)
                                                                                                      tttGTAATAATAGTGCTTATCCCGGTCGTTTATTTCGAGAATAACCAGACTAGGAACGCCAGCCGCCACTTCGTTTGGTGGAATGTCTTTTGTGACGATACTACACGCGCCAATATAAGAATTATCCCCGATGGTGACGCCTGGATtaat                                                 <  2:733850/150‑1 (MQ=255)
                                                                                                      tttGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAAACCCACAACGAACGCCAGCCGCCACCACGTTTGGTGGAATGTCTTTTGTGACGATACAACCCGCGCCAATAAAACAAATATACCCCATGGTGACGCCGGGATtaaa                                                 >  2:3795266/1‑149 (MQ=255)
                                                                                                      tttGAAATAATAGTGCTTATCCCGGTCGTTTATTTAGCGAATAAACCGACAAAGAAACCCAGCCCCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTAACCCCGCCAAAAACAGAATTATACCCGATGGTGAAGCCTGGATtaat                                                 >  2:1095407/1‑150 (MQ=255)
                                                                                                        tGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATtaataa                                               <  1:2256844/150‑1 (MQ=255)
                                                                                                         gAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAAc                                              <  2:10740597/150‑1 (MQ=255)
                                                                                                              aataGTGCTTATCCCGGTCGTTTATTTAGCGAAAAAACCCACAAGGAACGCCAGCCGCCACCCCGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAAATATCCCCCATGGTGACGCCTGGATTAATAACcacaa                                         >  2:3354108/1‑149 (MQ=255)
                                                                                                                  gTGCTTATCCCGGTCGTTTATTTCGCGAATAAACCGACAAGGAAACCCAACCCCCACGACGTTTGGTGGGATGTCTTTTTTTACGATACAAACCCCGCCAAAAACAGAATTATACCCGATGGTGACGCCGGGGTTAATAACcacaagtct                                     >  2:1964321/1‑145 (MQ=255)
                                                                                                                   tGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACtt                                    <  1:3423981/150‑1 (MQ=255)
                                                                                                                     cTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTcc                                  <  1:309299/150‑1 (MQ=255)
                                                                                                                      ttATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACCCCAGCCGCCACGCCGGTTGGTGGAAAGTCTTTTGTGACGATACAACCCCCGCCAAAAACACAAATATACCCCATGGTGACGCCAGGATTAATAAACACATGACATCCg                                 >  2:1504909/1‑150 (MQ=255)
                                                                                                                        aTCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGTATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTTATAACCACATGACTTCCGAt                               <  2:9742135/150‑1 (MQ=255)
                                                                                                                           ccGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCa                            <  1:9186124/150‑1 (MQ=255)
                                                                                                                           ccGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCa                            <  2:11095194/150‑1 (MQ=255)
                                                                                                                                    tATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTAtt                   <  2:680744/150‑1 (MQ=255)
                                                                                                                                     aTTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTg                  <  2:10420756/150‑1 (MQ=255)
                                                                                                                                        tCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCa               >  1:7203540/1‑150 (MQ=255)
                                                                                                                                        tCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCa               >  2:10328836/1‑150 (MQ=255)
                                                                                                                                         cgcgAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCaa              <  2:6954226/150‑1 (MQ=255)
                                                                                                                                          gcgAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAAt             >  1:11030405/1‑150 (MQ=255)
                                                                                                                                           cgAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATc            >  2:2120635/1‑150 (MQ=255)
                                                                                                                                           cgAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATc            <  2:6312384/150‑1 (MQ=255)
                                                                                                                                            gAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCg           <  1:8579985/150‑1 (MQ=255)
                                                                                                                                            gAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCg           <  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)
                                                                                                                                                      |                                                                                                                                                     
ATGCCGGATGCGGCTAATGTAGATCGCTGAACTTGTAGGCCTGATAAGCGCAGCGTATCAGGCAATTTTTATAATTTAAACTGACGATTCAACTTTATAATCTTTGAAATAATAGTGCTTATCCCGGTCGTTTATTTCGCGAATAACCCGACAAGGAACGCCAGCCGCCACGACGTTTGGTGGAATGTCTTTTGTGACGATACTACCCGCGCCAATAACAGAATTATCCCCGATGGTGACGCCTGGATTAATAACCACATGACTTCCGATCCAGACGTTATTGCCAATCGTTATCGGAA  >  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: