Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 700,980 G→C Y205* (TAC→TAG)  nagC ← N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913700,9800GC98.8% 310.0 / NA 81Y205* (TAC→TAG) nagCN‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (44/36);  total (44/37)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

CACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAAC  >  NC_000913/700832‑701123
                                                                                                                                                    |                                                                                                                                               
cACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAacta                                                                                                                                               <  2:17636/151‑1 (MQ=255)
       cGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTcc                                                                                                                                        >  1:1013250/1‑151 (MQ=255)
        ggTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGTAATCCTGACTTGCACCGAACTAGTGCTCcg                                                                                                                                       >  2:708940/1‑151 (MQ=255)
        ggTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCcg                                                                                                                                       >  2:4134/1‑151 (MQ=255)
        ggTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCcg                                                                                                                                       >  2:862740/1‑151 (MQ=255)
         gTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTcc                                                                                                                                        <  1:1051592/149‑1 (MQ=255)
         gTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCcgc                                                                                                                                      <  1:122839/151‑1 (MQ=255)
               ccTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCcgccagcgc                                                                                                                                >  2:424956/1‑151 (MQ=255)
                      atGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACt                                                                                                                          >  2:779622/1‑150 (MQ=255)
                      aaggCCAATATCACCGACGATGCCGTTGCGGTCAATAAAAATGTGCCCGTTAGAGATAATCTCGGCCCCGGTTGCGCGAGGGACACGCACCCAAATGGCGTCTTCGTCATCCCGACCTGCACCCAACAAGTGCTCCGCCAGCGCCAGACt                                                                                                                          >  1:442808/3‑150 (MQ=255)
                        ggCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTAc                                                                                                                        <  2:936142/150‑1 (MQ=255)
                           cAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGat                                                                                                                    <  2:879646/151‑1 (MQ=255)
                           cAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGat                                                                                                                    <  2:512014/151‑1 (MQ=255)
                               ttCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCg                                                                                                                >  1:223630/1‑151 (MQ=255)
                                   ccGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAc                                                                                                            >  2:146486/1‑151 (MQ=255)
                                    cgacgaTGACGTTGCGGCCAATAAAAATGCGAACGGTAGAGATAATACCGGCCCCGGTTGCGCGATGGACAAGCACGAGAAAGGAGTCTTCGCAATCCTGACTTGCACCGAACAAGTGCTCCGCCAGAGCCAGACAAAGGATATCgaga                                                                                                             >  1:777500/1‑149 (MQ=255)
                                     gACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCa                                                                                                          <  2:1103468/151‑1 (MQ=255)
                                            ccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAAc                                                                                                   >  1:768279/1‑151 (MQ=255)
                                                tgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGt                                                                                               >  2:433798/1‑151 (MQ=255)
                                                tgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGt                                                                                               >  1:39965/1‑151 (MQ=255)
                                                 gcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGt                                                                                               <  1:1224007/150‑1 (MQ=255)
                                                 gcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACAGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTc                                                                                              >  2:530707/1‑151 (MQ=255)
                                                 gcGGCCAATAAAAATGCGCCCGCTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTACGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTc                                                                                              >  1:273266/1‑151 (MQ=255)
                                                  cGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTc                                                                                              >  1:646337/1‑150 (MQ=255)
                                                   ggCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCa                                                                                             <  2:508600/150‑1 (MQ=255)
                                                        aTAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTaa                                                                                       <  1:986619/151‑1 (MQ=255)
                                                           aaaaTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGGCTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAAc                                                                                    <  1:962708/151‑1 (MQ=255)
                                                             aaTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGt                                                                                  >  1:1238240/1‑151 (MQ=255)
                                                             aaTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGt                                                                                  <  2:501138/151‑1 (MQ=255)
                                                                 cgcCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttctt                                                                              >  2:916065/1‑151 (MQ=255)
                                                                  gcCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGTCTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttcttc                                                                             <  1:1186949/151‑1 (MQ=255)
                                                                  gcCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGtt                                                                                 <  2:24365/147‑1 (MQ=255)
                                                                    ccGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTcc                                                                            >  2:793140/1‑150 (MQ=255)
                                                                     cGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCag                                                                          <  2:732914/151‑1 (MQ=255)
                                                                     cGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCag                                                                          <  1:4528/151‑1 (MQ=255)
                                                                         agagATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGc                                                                       <  1:929150/150‑1 (MQ=255)
                                                                         agagATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCt                                                                      <  1:698820/151‑1 (MQ=255)
                                                                          gagaTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCtt                                                                     >  1:995385/1‑151 (MQ=255)
                                                                            gaTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCt                                                                   >  2:133645/1‑151 (MQ=255)
                                                                             aTAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCt                                                                   >  1:96311/1‑150 (MQ=255)
                                                                                 tCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTAcc                                                                >  1:794283/1‑149 (MQ=255)
                                                                                  cccGGCCCCGGTTCCGCGATGGTCACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGc                                                             >  1:450015/1‑151 (MQ=255)
                                                                                   ccGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGcc                                                            <  2:495200/151‑1 (MQ=255)
                                                                                   ccGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGcc                                                            >  2:980691/1‑151 (MQ=255)
                                                                                     ggCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGcccc                                                          >  2:51185/1‑151 (MQ=255)
                                                                                      gCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGccccc                                                         <  1:256468/151‑1 (MQ=255)
                                                                                          cGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGtt                                                     <  1:133316/151‑1 (MQ=255)
                                                                                           ggatCCGCGATGGACACGCACCAGAATGGAGTCCTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACAAACGAAACAGGTCACTTTAAAACGCTCCTCCAGAGCTTCTACCAGCCCCCAGctt                                                    >  1:70641/4‑151 (MQ=255)
                                                                                           ggTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGttt                                                    >  2:765000/1‑151 (MQ=255)
                                                                                           ggTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGttt                                                    <  1:979856/151‑1 (MQ=255)
                                                                                            gTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGtttt                                                   <  2:1141768/151‑1 (MQ=255)
                                                                                              tCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTc                                                  <  2:155791/150‑1 (MQ=255)
                                                                                              tCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCt                                                 >  2:1163421/1‑151 (MQ=255)
                                                                                                cgcgATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTAc                                               >  2:421083/1‑151 (MQ=255)
                                                                                                        acacGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAatat                                       >  2:924933/1‑151 (MQ=255)
                                                                                                           cGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCg                                    <  1:661018/151‑1 (MQ=255)
                                                                                                             cACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCTAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCgg                                   <  1:856211/150‑1 (MQ=255)
                                                                                                              aCCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGc                                  >  2:981131/1‑150 (MQ=255)
                                                                                                              aCCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGc                                  >  1:461359/1‑150 (MQ=255)
                                                                                                               ccAGAAAGAAGTCTGCGCAATCCTGGGTAGCGCCAAAGGAGTGCTCCGCCAGCGCCAGACGGCGGATATCGTGGCCATCGAAGCAGGTCACTTTAAATCGTTCTTCCATTACTTCTACCAGCCCCCAGTTTTCTACCTCAATAtgctgcgt                                <  1:839816/151‑3 (MQ=255)
                                                                                                                cAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATg                               >  1:559027/1‑151 (MQ=255)
                                                                                                                    aTGGAGTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAAt                           <  2:200248/151‑1 (MQ=255)
                                                                                                                         gTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAAtt                      <  2:790810/151‑1 (MQ=255)
                                                                                                                         gTCTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAAt                       <  2:906949/150‑1 (MQ=255)
                                                                                                                           cTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAAtttt                    >  2:49630/1‑151 (MQ=255)
                                                                                                                           cTTCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAAtttt                    >  1:616699/1‑151 (MQ=255)
                                                                                                                             tCGCAATCCTGACTTGCGCCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAAttt                     <  1:593309/148‑1 (MQ=255)
                                                                                                                             tCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAAtttt                    <  1:464325/149‑1 (MQ=255)
                                                                                                                             tCGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAAttt                     <  2:189570/148‑1 (MQ=255)
                                                                                                                              cGCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGcc                 >  1:443115/1‑151 (MQ=255)
                                                                                                                               gCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCg                <  2:425968/151‑1 (MQ=255)
                                                                                                                               gCAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCg                <  2:383423/151‑1 (MQ=255)
                                                                                                                                cAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGc               >  2:734323/1‑151 (MQ=255)
                                                                                                                                cAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGc               >  1:206416/1‑151 (MQ=255)
                                                                                                                                cAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGc               >  1:1129032/1‑151 (MQ=255)
                                                                                                                                cAATCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGc               <  1:1119012/151‑1 (MQ=255)
                                                                                                                                 aaTCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCt              <  2:825313/151‑1 (MQ=255)
                                                                                                                                 aaTCCTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCt              >  2:96217/1‑151 (MQ=255)
                                                                                                                                    ccTGACTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTc           <  1:483254/151‑1 (MQ=255)
                                                                                                                                        aCTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGAGGCATGTAATGAATTTTGCCGCTGTCCggg       >  2:706693/1‑151 (MQ=255)
                                                                                                                                         cTTGCACCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGt      >  2:814829/1‑151 (MQ=255)
                                                                                                                                              aCCGAACTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTcaac  >  2:1065359/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTTCTACCTGAATATGCGGCATGTAATGAATTTTGCCGCTGTCCGGGTCAAC  >  NC_000913/700832‑701123

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