Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,371,068 A→G 7.1% intergenic (‑32/+16) 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,0680AG 7.1% 452.2 / 18.8 157intergenic (‑32/+16)nanE/nanTputative N‑acetylmannosamine‑6‑P epimerase/sialic acid transporter
Reads supporting (aligned to +/- strand):  ref base A (66/79);  new base G (4/7);  total (71/86)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 159 total aligned reads displayed.

ATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATAT  >  NC_000913/3370919‑3371217
                                                                                                                                                     |                                                                                                                                                     
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)
         cggcgACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcc                                                                                                                                              <  2:1055815/150‑2 (MQ=255)
          ggcgACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcct                                                                                                                                             <  2:1642101/150‑3 (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)
             gacgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGAGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcctttc                                                                                                                                          <  2:70141/150‑6 (MQ=255)
               cgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcctttcgc                                                                                                                                        <  1:853009/150‑8 (MQ=255)
                  tttCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTa                                                                                                                                     <  1:1702177/150‑1 (MQ=255)
                     cGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACt                                                                                                                                  >  1:574280/1‑150 (MQ=255)
                       ggTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACttt                                                                                                                                >  1:1600827/1‑150 (MQ=255)
                        gTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACtttt                                                                                                                               <  2:1308215/150‑1 (MQ=255)
                        gTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACtttt                                                                                                                               <  2:1633470/150‑1 (MQ=255)
                           tATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGt                                                                                                                            >  2:1149688/1‑150 (MQ=255)
                               gAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTg                                                                                                                        <  2:1240714/150‑1 (MQ=255)
                                 gcgggcggTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGAc                                                                                                                      <  2:238628/150‑1 (MQ=255)
                                 gCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGAc                                                                                                                      <  2:175381/150‑1 (MQ=255)
                                  cGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACt                                                                                                                     <  1:98744/150‑1 (MQ=255)
                                  cGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACt                                                                                                                     <  2:1600827/150‑1 (MQ=255)
                                      cTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAAt                                                                                                                 <  1:548278/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)
                                               aCCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTgg                                                                                                        >  2:133553/1‑150 (MQ=255)
                                               aCCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTgg                                                                                                        >  1:1583416/1‑150 (MQ=255)
                                               aCCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTgg                                                                                                        <  2:200248/150‑1 (MQ=255)
                                                  ggCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTTACGACATACATCTTCCCTTAGCGTAAGGCCCGGTACATTGCGCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGgcgc                                                                                                     >  1:1135873/1‑150 (MQ=255)
                                                   gCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGTTCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCt                                                                                                    >  2:302639/1‑150 (MQ=255)
                                                     tGGCAGGAGACAATCCGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTATATCGTTTTTGGCGCTGc                                                                                                  <  2:196411/150‑1 (MQ=255)
                                                     tGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGc                                                                                                  >  2:1583058/1‑150 (MQ=255)
                                                     tGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGc                                                                                                  >  2:1486534/1‑150 (MQ=255)
                                                      ggCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGcc                                                                                                 <  1:625131/150‑1 (MQ=255)
                                                      ggCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGcc                                                                                                 <  2:955509/150‑1 (MQ=255)
                                                      ggCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGcc                                                                                                 >  1:1442924/1‑150 (MQ=255)
                                                            agaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacgg                                                                                           >  1:1702381/1‑150 (MQ=255)
                                                             gaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggc                                                                                          >  2:759571/1‑150 (MQ=255)
                                                             gaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCAGCGCCaaaaacg                                                                                          <  1:1019554/150‑5 (MQ=255)
                                                             gaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCAGCGCCaaaaacg                                                                                          <  1:456917/150‑5 (MQ=255)
                                                                 aTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacgg                                                                                      >  2:1462389/1‑150 (MQ=255)
                                                                 aTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacgg                                                                                      <  1:485191/150‑1 (MQ=255)
                                                                     ggCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCAcc                                                                                  <  2:1526457/150‑1 (MQ=255)
                                                                     ggCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCAcc                                                                                  >  2:792315/1‑150 (MQ=255)
                                                                     ggCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCAcc                                                                                  >  2:1286980/1‑150 (MQ=255)
                                                                         aCCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTg                                                                              >  1:1401123/1‑150 (MQ=255)
                                                                            gTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAAt                                                                           <  1:986807/150‑1 (MQ=255)
                                                                                 cagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTTCATAGAACGGGCAACAGGATTAACTTTTGGTTTTGACTAAATAGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGttt                                                                      >  1:248015/1‑150 (MQ=255)
                                                                                 cagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGttt                                                                      <  1:133553/150‑1 (MQ=255)
                                                                                 cagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGttt                                                                      >  1:1554013/1‑150 (MQ=255)
                                                                                     gATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCg                                                                  >  1:1654188/1‑150 (MQ=255)
                                                                                           ttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATa                                                            <  1:1198668/150‑1 (MQ=255)
                                                                                            tGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAg                                                           >  1:770597/1‑150 (MQ=255)
                                                                                              aTCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCg                                                         <  1:458385/150‑1 (MQ=255)
                                                                                               tCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCgg                                                        >  2:537210/1‑149 (MQ=255)
                                                                                                ccAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTc                                                       <  2:442652/150‑1 (MQ=255)
                                                                                                 cAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCa                                                      <  2:503864/150‑1 (MQ=255)
                                                                                                 cAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCa                                                      <  2:1567899/150‑1 (MQ=255)
                                                                                                 cAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCa                                                      >  2:809579/1‑150 (MQ=255)
                                                                                                  aGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCTGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCAt                                                     >  1:314781/1‑150 (MQ=255)
                                                                                                   gTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATg                                                    <  2:1554576/150‑1 (MQ=255)
                                                                                                    ttGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGa                                                   <  2:1348945/150‑1 (MQ=255)
                                                                                                    ttGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGa                                                   <  2:1598942/150‑1 (MQ=255)
                                                                                                     tgtgCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAg                                                  <  1:1486534/150‑1 (MQ=255)
                                                                                                           aGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCa                                            >  1:246578/1‑150 (MQ=255)
                                                                                                               aCGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGc                                        <  2:1401123/150‑1 (MQ=255)
                                                                                                               aCGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGc                                        <  1:829644/150‑1 (MQ=255)
                                                                                                                cGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCt                                       <  1:309751/150‑1 (MQ=255)
                                                                                                                 gACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATg                                                    <  1:513720/136‑1 (MQ=255)
                                                                                                                  acatacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTc                                     <  1:365251/150‑1 (MQ=255)
                                                                                                                  acatacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTcc                                    >  2:31071/1‑150 (MQ=255)
                                                                                                                   catacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTcc                                    >  2:1358130/1‑150 (MQ=255)
                                                                                                                      acatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCggg                                 <  2:1702381/150‑1 (MQ=255)
                                                                                                                       catCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTACGGgc                                >  1:267023/1‑150 (MQ=255)
                                                                                                                         tCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGgcgc                              <  1:1286999/150‑1 (MQ=255)
                                                                                                                             cccTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGAGTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaacc                          >  1:734069/1‑150 (MQ=255)
                                                                                                                             cccTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaacc                          >  2:411821/1‑150 (MQ=255)
                                                                                                                                ttAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaaccaac                       <  1:1462389/150‑1 (MQ=255)
                                                                                                                                   gCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCt                    <  2:944820/150‑1 (MQ=255)
                                                                                                                                   gCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCt                    <  1:799299/150‑1 (MQ=255)
                                                                                                                                    cGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCCCCGCTGAATGGTTTACCGTCGCTAGCGTCATGAGTACGCAACGCGTCCGGGCGCAACCAACGCTg                   >  2:1006305/1‑150 (MQ=255)
                                                                                                                                    cGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTg                   >  2:1155780/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)
                                                                                                                                        aGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcg               <  1:1377814/150‑1 (MQ=255)
                                                                                                                                        aGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGATTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcg               <  2:700674/150‑1 (MQ=255)
                                                                                                                                         tgCCCGGTCTATGTACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCATCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgc              >  1:671319/2‑150 (MQ=255)
                                                                                                                                         ggCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgc              >  2:1104246/1‑150 (MQ=255)
                                                                                                                                          gCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgcg             <  1:1104246/150‑1 (MQ=255)
                                                                                                                                           cccGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAAGCGTTTTTGGCGCTGCCAAACGGCACGGCCCCGCTGAATGGTTTACCGTCGATAGCGTCAAGAGTACGCAACGCTTCCGGTCGCAACCAAACCTGAACGCga            >  1:88531/1‑150 (MQ=255)
                                                                                                                                            ccGGTACATAGACCGGGCAACAGTATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCATCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCgag           <  2:471464/150‑1 (MQ=255)
                                                                                                                                              ggTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCAGGGCGCAACCAACGCTGAACGCGAGaa         >  1:346041/1‑150 (MQ=255)
                                                                                                                                               gTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAg        >  1:1245558/1‑150 (MQ=255)
                                                                                                                                                    tAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCata   >  1:903603/1‑150 (MQ=255)
                                                                                                                                                     aGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCatat  <  2:555309/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGGCATAT  >  NC_000913/3370919‑3371217

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