Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,371,064 A→C 6.9% intergenic (‑28/+20) nanE ← / ← nanT putative N‑acetylmannosamine‑6‑P epimerase/sialic acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,371,0640AC 6.9% 478.9 / 18.8 154intergenic (‑28/+20)nanE/nanTputative N‑acetylmannosamine‑6‑P epimerase/sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base A (62/80);  new base C (4/7);  total (67/87)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 155 total aligned reads displayed.

GCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGG  >  NC_000913/3370915‑3371212
                                                                                                                                                     |                                                                                                                                                    
    aTGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATa                                                                                                                                                  <  2:583/150‑1 (MQ=255)
         aTGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCg                                                                                                                                             <  1:1486841/150‑1 (MQ=255)
          tGGCGGCGACGATTTCGGGTTTCTCTAGCGGGCTGTCCGGAACCGGCGGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCgg                                                                                                                                            <  2:163749/150‑1 (MQ=255)
           ggcggcGTCGATTTCGGGTTTATCGTGCGGGCTGTCCGGAACCGGCTGGGAGGATACAATCAGGCCACCGTTAGCATCGATTTTTTTATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAACGGCCCGGTACATAGACCggg                                                                                                                                           <  2:1059238/150‑1 (MQ=255)
            gcggcgACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGc                                                                                                                                          <  1:1012327/150‑1 (MQ=255)
            gcggcgACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGc                                                                                                                                          <  1:628600/150‑1 (MQ=255)
              ggcgACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGCCCGGGCaa                                                                                                                                        >  2:1286999/1‑150 (MQ=255)
              ggcgACGATTTCGGGTTTATCGAGCGGGCGGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCTTTAGCAGCGATGTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGCTCTATGTACCGGGcct                                                                                                                                        <  1:134633/150‑3 (MQ=255)
                cgacgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACa                                                                                                                                      >  2:868986/1‑150 (MQ=255)
                cgacgaTTTCGGGGTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACACTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGCCCGGGCAACa                                                                                                                                      >  2:844131/1‑150 (MQ=255)
                 gacgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCCGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAg                                                                                                                                     >  1:711812/1‑150 (MQ=255)
                   cgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcctttcgc                                                                                                                                   <  1:853009/150‑8 (MQ=255)
                   cgaTTTCGGGTTTATCGAGCGGGCTATCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTATCGAAAGGCCCGGTCTATGTACCGGGcctttcgc                                                                                                                                   <  1:555101/150‑8 (MQ=255)
                    gaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGAt                                                                                                                                  <  2:1555571/150‑1 (MQ=255)
                      tttCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTa                                                                                                                                <  1:1702177/150‑1 (MQ=255)
                         cGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACt                                                                                                                             >  1:574280/1‑150 (MQ=255)
                           ggTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACttt                                                                                                                           >  1:1600827/1‑150 (MQ=255)
                           ggTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACttt                                                                                                                           >  2:456917/1‑150 (MQ=255)
                           ggTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACttt                                                                                                                           >  2:1019554/1‑150 (MQ=255)
                            gTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACtttt                                                                                                                          <  2:1595105/150‑1 (MQ=255)
                             tttATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTGATTCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTgg                                                                                                                        <  1:805495/150‑1 (MQ=255)
                               tATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGt                                                                                                                       >  2:1149688/1‑150 (MQ=255)
                                     gcgggcggTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGAc                                                                                                                 <  2:238628/150‑1 (MQ=255)
                                     gCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGAc                                                                                                                 <  2:175381/150‑1 (MQ=255)
                                      cGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACt                                                                                                                <  2:1600827/150‑1 (MQ=255)
                                      cGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACt                                                                                                                <  1:98744/150‑1 (MQ=255)
                                           tGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATc                                                                                                           <  2:92959/150‑1 (MQ=255)
                                           tGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATc                                                                                                           >  1:17114/1‑150 (MQ=255)
                                                 gAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGttttt                                                                                                     <  1:1376503/150‑1 (MQ=255)
                                                 gAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGttttt                                                                                                     <  2:1370268/150‑1 (MQ=255)
                                                  aaCCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTg                                                                                                    <  2:193064/150‑1 (MQ=255)
                                                   aCCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTgg                                                                                                   >  2:133553/1‑150 (MQ=255)
                                                   aCCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTgg                                                                                                   <  2:200248/150‑1 (MQ=255)
                                                    ccGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTCACGACATACATCTTCCCATAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGgc                                                                                                  <  2:1089817/150‑1 (MQ=255)
                                                      ggCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTTACGACATACATCTTCCCTTAGCGTAAGGCCCGGTACATTGCGCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGgcgc                                                                                                >  1:1135873/1‑150 (MQ=255)
                                                       gCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGTTCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCt                                                                                               >  2:302639/1‑150 (MQ=255)
                                                          ggCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGcc                                                                                            <  2:955509/150‑1 (MQ=255)
                                                          ggCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGcc                                                                                            >  1:1442924/1‑150 (MQ=255)
                                                                agaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacgg                                                                                      >  1:1702381/1‑150 (MQ=255)
                                                                agaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCAGCGCCaaaaac                                                                                      <  1:415486/150‑4 (MQ=255)
                                                                 gaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggc                                                                                     >  2:759571/1‑150 (MQ=255)
                                                                 gaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCAGCGCCaaaaacg                                                                                     <  1:1019554/150‑5 (MQ=255)
                                                                    aaTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacg                                                                                  <  1:464347/150‑1 (MQ=255)
                                                                     aTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGc                                                                                                                                          >  2:1379431/1‑93 (MQ=255)
                                                                     aTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGc                                                                                                                                          <  1:1379431/93‑1 (MQ=255)
                                                                     aTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacgg                                                                                 >  2:1462389/1‑150 (MQ=255)
                                                                     aTCAGGCCACCGTAAGCAGCGATTATTTGACCCAGTTCGCCAAGTCCCGACATACATCATCCCTTAGCGAAAGGCCCGGCACATAGACCGGGCAACAGGATTAACTTTTGGTATTGACTAAATCGTTTTTGGCGCTGCCAAacggcacgg                                                                                 <  1:326013/150‑1 (MQ=255)
                                                                         ggCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCAcc                                                                             <  2:1526457/150‑1 (MQ=255)
                                                                         ggCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCAcc                                                                             >  2:1286980/1‑150 (MQ=255)
                                                                             aCCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTg                                                                         >  1:1401123/1‑150 (MQ=255)
                                                                                gTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAAt                                                                      >  2:799299/1‑150 (MQ=255)
                                                                                    gcagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGACAAACGGCACGACACCGCTGAATGGtt                                                                  >  1:1290700/1‑150 (MQ=255)
                                                                                     cagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTTCATAGAACGGGCAACAGGATTAACTTTTGGTTTTGACTAAATAGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGttt                                                                 >  1:248015/1‑150 (MQ=255)
                                                                                     cagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGttt                                                                 >  1:1554013/1‑150 (MQ=255)
                                                                                     cagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGttt                                                                 <  1:133553/150‑1 (MQ=255)
                                                                                        cGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTAcc                                                              >  1:944820/1‑150 (MQ=255)
                                                                                         gATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCg                                                             >  1:1654188/1‑150 (MQ=255)
                                                                                            tttttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCg                                                          <  1:204078/150‑1 (MQ=255)
                                                                                             ttttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGa                                                         >  2:757197/1‑150 (MQ=255)
                                                                                               ttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATa                                                       <  1:1198668/150‑1 (MQ=255)
                                                                                               ttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATa                                                       <  2:1482940/150‑1 (MQ=255)
                                                                                                tGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAg                                                      >  1:770597/1‑150 (MQ=255)
                                                                                                  aTCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCg                                                    <  1:458385/150‑1 (MQ=255)
                                                                                                    ccAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTc                                                  <  2:442652/150‑1 (MQ=255)
                                                                                                    ccAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTc                                                  <  1:844131/150‑1 (MQ=255)
                                                                                                     cAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCa                                                 >  2:809579/1‑150 (MQ=255)
                                                                                                     cAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCa                                                 <  2:503864/150‑1 (MQ=255)
                                                                                                     cAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCa                                                 <  2:1567899/150‑1 (MQ=255)
                                                                                                        ttGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGa                                              <  2:1348945/150‑1 (MQ=255)
                                                                                                         tgtgCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAg                                             <  1:1486534/150‑1 (MQ=255)
                                                                                                           tgCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTa                                           >  1:259342/1‑150 (MQ=255)
                                                                                                            gCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTAc                                          <  1:493717/150‑1 (MQ=255)
                                                                                                                 tAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATTGCGTCATGAGTACGCAAc                                     <  2:1635113/150‑1 (MQ=255)
                                                                                                                    cGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCt                                  <  1:309751/150‑1 (MQ=255)
                                                                                                                     gACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATg                                               <  1:513720/136‑1 (MQ=255)
                                                                                                                     gACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATg                                               >  2:513720/1‑136 (MQ=255)
                                                                                                                      acatacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTc                                <  1:365251/150‑1 (MQ=255)
                                                                                                                      acatacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTcc                               >  2:31071/1‑150 (MQ=255)
                                                                                                                       catacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTcc                               >  2:1358130/1‑150 (MQ=255)
                                                                                                                        atacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCg                              <  2:574280/150‑1 (MQ=255)
                                                                                                                            atCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTAAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGgcg                          >  2:1527893/1‑150 (MQ=255)
                                                                                                                             tCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGgcgc                         <  1:1286999/150‑1 (MQ=255)
                                                                                                                                 cccTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaacc                     >  2:411821/1‑150 (MQ=255)
                                                                                                                                    ttAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaaccaac                  <  1:1462389/150‑1 (MQ=255)
                                                                                                                                      aGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGc                <  1:1583058/150‑1 (MQ=255)
                                                                                                                                      aGCGAAAGGCCCGGTACATAGACAGGGCAACAGGATTAACTGTTGGTTTTGACTAAATCGTTTTTGGCGCTGCAAAACGGCACGGCACCGCTGAATGGTTTACCGACGATAGCGTCAAGAGTACGCAACGCTTCCGGGCGCAACCAACGc                >  1:68276/1‑150 (MQ=255)
                                                                                                                                       gCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCt               <  1:759571/150‑1 (MQ=255)
                                                                                                                                       gCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCt               <  2:944820/150‑1 (MQ=255)
                                                                                                                                        cGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCCCCGCTGAATGGTTTACCGTCGCTAGCGTCATGAGTACGCAACGCGTCCGGGCGCAACCAACGCTg              >  2:1006305/1‑150 (MQ=255)
                                                                                                                                        cGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTg              >  2:1155780/1‑150 (MQ=255)
                                                                                                                                          aaaGGCCCGGTACATAGACCGTGCAACAGGATTAACTTTTGGTTGTGACTAAATCGTTTTTGGCGATGACAAACGGCACGGCTCCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGaa            >  2:1630359/1‑150 (MQ=255)
                                                                                                                                           gttgCCCGGTCTATGTACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAc           >  1:1642718/4‑150 (MQ=255)
                                                                                                                                           gttgCCCGGTCTATGTACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAc           >  1:827563/4‑150 (MQ=255)
                                                                                                                                            aGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcg          >  2:1662594/1‑150 (MQ=255)
                                                                                                                                            aGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGATTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcg          <  2:700674/150‑1 (MQ=255)
                                                                                                                                              gCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgcg        <  1:1104246/150‑1 (MQ=255)
                                                                                                                                               cccGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAAGCGTTTTTGGCGCTGCCAAACGGCACGGCCCCGCTGAATGGTTTACCGTCGATAGCGTCAAGAGTACGCAACGCTTCCGGTCGCAACCAAACCTGAACGCga       >  1:88531/1‑150 (MQ=255)
                                                                                                                                                ccGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCgag      >  2:1661227/1‑150 (MQ=255)
                                                                                                                                                   gTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAg   >  1:1245558/1‑150 (MQ=255)
                                                                                                                                                   gTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAgg  <  1:31071/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGG  >  NC_000913/3370915‑3371212

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