Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,659,314 A→G T416A (ACC→GCC)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,659,3140AG98.5% 247.0 / NA 68T416A (ACC→GCC) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (39/28);  minor base C (1/0);  total (40/28)
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.96e-01

CCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCC  >  NC_000913/1659173‑1659460
                                                                                                                                             |                                                                                                                                                  
cccGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTAtgccattacca                                                                                                                                           <  2:272990/151‑1 (MQ=255)
     ttATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGa                                                                                                                                       >  1:488892/1‑150 (MQ=255)
       atatCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAAGGTCGGCATCAGCGGCGGAAATAGGGGCGCGCGTGAATCGACCTATGCCATTTCCATAGAAc                                                                                                                                     >  2:501277/1‑150 (MQ=255)
             gCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGcc                                                                                                                                  <  2:475736/147‑1 (MQ=255)
             gCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGcc                                                                                                                              <  2:198883/151‑1 (MQ=255)
              ccAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCg                                                                                                                             <  1:537613/151‑1 (MQ=255)
                aGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCgg                                                                                                                            >  2:751888/1‑150 (MQ=255)
                   gCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCt                                                                                                                        <  2:321493/151‑1 (MQ=255)
                    cTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTg                                                                                                                       >  1:307362/1‑151 (MQ=255)
                      gggggCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTgg                                                                                                                      >  1:205206/1‑150 (MQ=255)
                       ggggCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGa                                                                                                                     >  2:652664/1‑150 (MQ=255)
                             acaACGCCAGGCAAACGGCGAACTCACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCg                                                                                                              <  2:783884/151‑1 (MQ=255)
                              caACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCgg                                                                                                             <  1:110584/151‑1 (MQ=255)
                              caACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCCATTTTGACGGGCAATGTCGGCATCAGCTGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTCCCATAGAACGCCTGCCGGTGCTGGATAAt                                                                                                                 <  2:632807/147‑1 (MQ=255)
                                          aaCGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCaa                                                                                                 <  1:848562/151‑1 (MQ=255)
                                            cGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAtc                                                                                               <  2:732274/151‑1 (MQ=255)
                                              gCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAAtctc                                                                                             >  2:701776/1‑151 (MQ=255)
                                                gAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCt                                                                                           <  1:258949/151‑1 (MQ=255)
                                                   ctgactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCt                                                                                        >  1:82256/1‑151 (MQ=255)
                                                    tgactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCtt                                                                                       >  2:573419/1‑151 (MQ=255)
                                                     gactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTc                                                                                      <  1:811962/151‑1 (MQ=255)
                                                     gactgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTc                                                                                      >  2:335215/1‑151 (MQ=255)
                                                      actgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTc                                                                                      >  2:738055/1‑150 (MQ=255)
                                                      actgCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCa                                                                                     <  2:759503/151‑1 (MQ=255)
                                                         gCACGCGCTATTGCTATGTTACCTATTTTGAAGTGCAATGTCGGCATCAGCGTCGGAAATAGTGGCGCGCGTGACTCGACCTATCCCAATCCCATAGAACGCCCGCCGGTGCTGGATATTCCGGTCACAACGACAATTCACTGCTTCAGc                                                                                   >  2:559442/1‑150 (MQ=255)
                                                            cgcgCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTgg                                                                                >  1:272121/1‑150 (MQ=255)
                                                             gcgcTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGAc                                                                              <  1:665750/151‑1 (MQ=255)
                                                             gcgcTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGAc                                                                              <  1:264691/151‑1 (MQ=255)
                                                              cgcTATTGCTATGTTACCTATTTTGCCGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACa                                                                             <  2:661199/151‑1 (MQ=255)
                                                                 tATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGAc                                                                              >  1:765416/1‑147 (MQ=255)
                                                                  aTTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATg                                                                         <  2:569462/151‑1 (MQ=255)
                                                                  aTTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATg                                                                         >  1:170584/1‑151 (MQ=255)
                                                                     gCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTa                                                                      <  2:327950/151‑1 (MQ=255)
                                                                       tATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCa                                                                                                  >  1:5554/1‑121 (MQ=255)
                                                                       tATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCa                                                                                                  <  2:5548/121‑1 (MQ=255)
                                                                        aTGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTAt                                                                     <  1:734970/149‑1 (MQ=255)
                                                                        aTGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCg                                                                   >  1:375811/1‑151 (MQ=255)
                                                                              ccTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATg                                                             >  2:615216/1‑151 (MQ=255)
                                                                                tATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGt                                                           >  2:786956/1‑151 (MQ=255)
                                                                                      gACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCa                                                      <  1:461170/150‑1 (MQ=255)
                                                                                       aCGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCaa                                                     <  2:801651/150‑1 (MQ=255)
                                                                                       aCGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCaa                                                     >  1:657280/1‑150 (MQ=255)
                                                                                       aCGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCaa                                                     >  1:628319/1‑150 (MQ=255)
                                                                                       aCGGGCAATGTAGGCATCAGCGGCGGAAATAGTGGCGCACGTGAATCGACCTATGCCATTACCCAAGAACGCCTGCCGGTGCTGGATAATACGGTCAAAACGTCAATCTCTTGCTTCAGCTGGACGGATTCTATCGAACATGGTACGCaa                                                     >  2:906695/1‑150 (MQ=255)
                                                                                        cGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAAt                                                   >  1:337595/1‑151 (MQ=255)
                                                                                           gCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGAc                                                >  2:60512/1‑151 (MQ=255)
                                                                                            cAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACg                                               >  2:79506/1‑151 (MQ=255)
                                                                                               tGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGc                                             >  1:229067/1‑150 (MQ=255)
                                                                                                  cGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATc                                         <  1:235683/151‑1 (MQ=255)
                                                                                                     cATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCcgc                                      >  1:458562/1‑151 (MQ=255)
                                                                                                     cATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCcgc                                      <  1:514562/151‑1 (MQ=255)
                                                                                                        cAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGa                                    <  2:417099/150‑1 (MQ=255)
                                                                                                        cAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGAc                                   <  2:889923/151‑1 (MQ=255)
                                                                                                        cAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGAc                                   >  2:458134/1‑151 (MQ=255)
                                                                                                         aGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACg                                  >  2:741011/1‑151 (MQ=255)
                                                                                                         aGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGACTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACg                                  >  2:668509/1‑151 (MQ=255)
                                                                                                                gAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcg                           <  1:79528/151‑1 (MQ=255)
                                                                                                                gAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcg                           <  2:661152/151‑1 (MQ=255)
                                                                                                                gAAATAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGCCcg                           >  1:203710/1‑151 (MQ=255)
                                                                                                                 aaaTAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcgc                          >  1:303439/1‑151 (MQ=255)
                                                                                                                  aaTAGTGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCcgcg                         >  2:513671/1‑151 (MQ=255)
                                                                                                                       tGGCGCGCGTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCaa                     >  1:60694/1‑150 (MQ=255)
                                                                                                                           gcgcgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGAt                 <  1:31817/150‑1 (MQ=255)
                                                                                                                            cgcgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATa                >  2:491833/1‑150 (MQ=255)
                                                                                                                              cgTGAATCGACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAAc             >  1:327205/1‑151 (MQ=255)
                                                                                                                                tGAATCGACCTATGCGATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTg           >  2:341842/1‑151 (MQ=255)
                                                                                                                                      gACCTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGAtgtg     >  1:161385/1‑151 (MQ=255)
                                                                                                                                        ccTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGcc   >  1:213397/1‑151 (MQ=255)
                                                                                                                                         cTATGCCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGccc  >  1:655219/1‑151 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
CCCGCTTATATCTGCCAGGGCTGGGGGCCACAACGCCAGGCAAACGGCGAACTGACTGCACGCGCTATTGCTATGTTACCTATTTTGACGGGCAATGTCGGCATCAGCGGCGGAAATAGTGGCGCGCGTGAATCGACCTATACCATTACCATAGAACGCCTGCCGGTGCTGGATAATCCGGTCAAAACGTCAATCTCCTGCTTCAGCTGGACAGATGCTATCGATCATGGTCCGCAAATGACGGCAATCCGCGACGGCGTCCGCGGCAAAGATAAACTGGATGTGCCC  >  NC_000913/1659173‑1659460

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: ‑