Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,371,064 A→C 9.3% intergenic (‑28/+20) nanE ← / ← nanT putative N‑acetylmannosamine‑6‑phosphate epimerase/N‑acetylneuraminate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,371,0640AC 9.3% 243.9 / 14.7 85intergenic (‑28/+20)nanE/nanTputative N‑acetylmannosamine‑6‑phosphate epimerase/N‑acetylneuraminate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base A (41/36);  new base C (4/4);  total (45/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAGG  >  NC_000913/3370915‑3371212
                                                                                                                                                     |                                                                                                                                                    
gCTAATGCCATGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTa                                                                                                                                                      >  1:1424757/1‑150 (MQ=255)
      gCCATGGTGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGTCCCGGTACATAGa                                                                                                                                                >  1:1586409/1‑150 (MQ=255)
         aTGGCGGCGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCg                                                                                                                                             <  1:984333/150‑1 (MQ=255)
           ggcggcGACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCggg                                                                                                                                           <  2:464754/150‑1 (MQ=255)
             cggcgACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCa                                                                                                                                         <  1:1337178/150‑1 (MQ=255)
              ggcgACGATTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcct                                                                                                                                        <  2:781351/150‑3 (MQ=255)
                cgacgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACa                                                                                                                                      <  2:1445019/150‑1 (MQ=255)
                cgacgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACa                                                                                                                                      >  1:1350642/1‑150 (MQ=255)
                  acgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcctttcg                                                                                                                                    <  2:50515/150‑7 (MQ=255)
                  acgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAgg                                                                                                                                    <  2:1621298/150‑1 (MQ=255)
                   cgaTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTCTATGTACCGGGcctttcgc                                                                                                                                   <  2:846701/150‑8 (MQ=255)
                     ttttCTGTTTTATCGAGCGGGCTGTCCGGAACCGTCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTTATCCATTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGACCGGTACATAGACCGGGCAACAGGAtt                                                                                                                                 <  2:1586409/149‑1 (MQ=255)
                     aTTTCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGAtt                                                                                                                                 <  2:505428/150‑1 (MQ=255)
                      tttCGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAATTAACGACATACATCTTCCCTTAGCGAAATGCCCTTTACATAGATCGGGCAACAGTATTa                                                                                                                                <  2:1523733/150‑1 (MQ=255)
                        tCGGGTTCATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAAc                                                                                                                              <  1:1522438/150‑1 (MQ=255)
                         cGGGTTTATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACt                                                                                                                             <  2:494823/150‑1 (MQ=255)
                               tATCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGt                                                                                                                       >  2:894085/1‑150 (MQ=255)
                                 tCGAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGttt                                                                                                                     <  1:1447584/150‑1 (MQ=255)
                                   gAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTg                                                                                                                   <  2:959356/150‑1 (MQ=255)
                                   gAGCGGGCTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTg                                                                                                                   <  2:1330825/150‑1 (MQ=255)
                                          cTGTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAAt                                                                                                            >  2:171466/1‑150 (MQ=255)
                                            gTCCGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCg                                                                                                          >  1:743397/1‑150 (MQ=255)
                                              ccGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTg                                                                                                                         >  1:926413/1‑133 (MQ=255)
                                              ccGGAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTg                                                                                                                         <  2:926413/133‑1 (MQ=255)
                                                ggAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGtttt                                                                                                      <  1:407550/150‑1 (MQ=255)
                                                 gAACCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGttttt                                                                                                     <  1:607660/150‑1 (MQ=255)
                                                   aCCGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTgg                                                                                                   <  2:1006471/150‑1 (MQ=255)
                                                    ccGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGgc                                                                                                  <  2:539469/150‑1 (MQ=255)
                                                     cGGCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGgcg                                                                                                 >  2:1612930/1‑150 (MQ=255)
                                                      ggCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGgcgc                                                                                                <  2:736139/150‑1 (MQ=255)
                                                       gCTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCt                                                                                               >  1:1323842/1‑150 (MQ=255)
                                                        cTGGCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTg                                                                                              <  1:502688/150‑1 (MQ=255)
                                                          ggCAGGAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGcc                                                                                            <  1:382217/150‑1 (MQ=255)
                                                             aGGACACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACAACATACATCTTACCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCaaa                                                                                         >  2:705189/1‑150 (MQ=255)
                                                              ggAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAac                                                                                        >  1:1569818/1‑150 (MQ=255)
                                                              ggAGACAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTACCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAac                                                                                        >  2:1336211/1‑150 (MQ=255)
                                                                agaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacgg                                                                                      <  1:120866/150‑1 (MQ=255)
                                                                 gaCAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggc                                                                                     >  1:252372/1‑150 (MQ=255)
                                                                   cAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcac                                                                                   >  2:669907/1‑150 (MQ=255)
                                                                   cAATCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcac                                                                                   >  1:1119143/1‑150 (MQ=255)
                                                                    aaTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacg                                                                                  <  1:1612930/150‑1 (MQ=255)
                                                                    aaTCAGGCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAacggcacg                                                                                  >  1:1146737/1‑150 (MQ=255)
                                                                         ggCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCAcc                                                                             >  1:206566/1‑150 (MQ=255)
                                                                         ggCCACCGTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCAcc                                                                             >  2:213980/1‑150 (MQ=255)
                                                                            cACCGTTAGCAGCGATTTTTTGCTCCAGTGGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCt                                                                          <  2:1578211/150‑1 (MQ=255)
                                                                                gTTAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAAt                                                                      <  2:130632/150‑1 (MQ=255)
                                                                                  tAGCAGCGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATgg                                                                    <  1:639774/150‑1 (MQ=255)
                                                                                     cagcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGTCGCTGCCAAACGGCACGGCACCGCTTATTGGttt                                                                 <  2:743397/150‑1 (MQ=255)
                                                                                       gcGATTTTTTGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTAc                                                               >  1:1517510/1‑150 (MQ=255)
                                                                                           ttttttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTc                                                           >  1:474929/1‑150 (MQ=255)
                                                                                             ttttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGa                                                         >  2:871722/1‑150 (MQ=255)
                                                                                              tttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGAt                                                        <  2:1350642/150‑1 (MQ=255)
                                                                                               ttGATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATa                                                       <  2:1424757/150‑1 (MQ=255)
                                                                                                 gATCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACAGTCGATAGc                                                     >  2:1014606/1‑150 (MQ=255)
                                                                                                  aTCCAGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCg                                                    >  1:475460/1‑150 (MQ=255)
                                                                                                   tCCAGTTGTGCAAGTAACGACATACATCCTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGt                                                   <  1:39737/150‑1 (MQ=255)
                                                                                                      aGTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCAt                                                >  2:1471139/1‑150 (MQ=255)
                                                                                                       gTTGTGCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATg                                               <  1:600172/150‑1 (MQ=255)
                                                                                                          gtgCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGt                                            >  2:997570/1‑150 (MQ=255)
                                                                                                           tgCAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTa                                           >  2:549381/1‑150 (MQ=255)
                                                                                                             cAAGTAACGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACg                                         >  1:30310/1‑150 (MQ=255)
                                                                                                                    cGACATACATCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGATTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCt                                  <  1:661409/150‑1 (MQ=255)
                                                                                                                       catacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTcc                               >  2:732133/1‑150 (MQ=255)
                                                                                                                         tacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCCCCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCgg                             >  2:662825/1‑150 (MQ=255)
                                                                                                                         tacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCgg                             >  1:1031499/1‑150 (MQ=255)
                                                                                                                         tacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCgg                             >  1:593049/1‑150 (MQ=255)
                                                                                                                         tacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCgg                             <  2:1323842/150‑1 (MQ=255)
                                                                                                                         tacatCTTCCCTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCgg                             >  2:601728/1‑150 (MQ=255)
                                                                                                                                  ccTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaacca                    >  1:1612008/1‑150 (MQ=255)
                                                                                                                                   cTTAGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGcaaccaa                   >  1:250366/1‑150 (MQ=255)
                                                                                                                                      aGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGc                >  1:689389/1‑150 (MQ=255)
                                                                                                                                      aGCGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGAAACCCTTCCGCGCGCAACCAACGc                >  2:1322095/1‑150 (MQ=255)
                                                                                                                                       tcctgttgCCCGGTCTATGTACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCt               >  2:80895/8‑150 (MQ=255)
                                                                                                                                        cctgttgCCCGGTCTATGTACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTg              >  2:1503377/7‑150 (MQ=255)
                                                                                                                                        cGAAAGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTg              >  2:306776/1‑150 (MQ=255)
                                                                                                                                         ctgttgCCCGGTCTATGTACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCCACCAACGCTGa             >  1:131048/6‑150 (MQ=255)
                                                                                                                                          tgttgCCCGGTCTATGTACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGATTCCGGGCGCAACCAACGCTGaa            >  1:1463157/5‑150 (MQ=255)
                                                                                                                                           aaGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAc           <  1:306271/150‑1 (MQ=255)
                                                                                                                                           aaGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAc           <  1:853916/150‑1 (MQ=255)
                                                                                                                                            aGGCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTATTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTAACCCTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGATGAAcc          >  2:1294189/1‑149 (MQ=255)
                                                                                                                                              gCCCGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAAcgcg        >  2:1534997/1‑150 (MQ=255)
                                                                                                                                               cccGGTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCga       >  1:315627/1‑150 (MQ=255)
                                                                                                                                               cccGGTACATAGACCGGGCAACACGATTAACTTTTGGTATTGAATAAATCGTTTTTGGCGCTGCCAAACCGCACGGCACCGCTGAATGGATTACCGTCGATAGAGTCATCACTAAGCAACGCATCCAGGCGCAACCAACGATCAACGCga       >  1:731669/1‑150 (MQ=255)
                                                                                                                                                  ggTACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACCGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGaa    >  2:864804/1‑150 (MQ=255)
                                                                                                                                                    tACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAgg  <  1:1561074/150‑1 (MQ=255)
                                                                                                                                                    tACATAGACCGGGCAACAGGATTAACTTTTGGTTTTGACTAAATCGTTTTTGGCGCTGCCAAACGGCACGGCACCGCTGAATGGTTTACCGTCGATAGCGTCATGAGTACGCAACGCTTCCGGGCGCAACCAACGCTGAACGCGAGAAgg  <  1:1555448/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 < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: