Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,407,641 A→G V52A (GTC→GCC)  ydaM ← diguanylate cyclase, csgD regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,407,6410AG100.0% 266.5 / NA 72V52A (GTC→GCC) ydaMdiguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (40/32);  total (40/32)

CACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGC  >  NC_000913/1407493‑1407787
                                                                                                                                                    |                                                                                                                                                  
cACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCg                                                                                                                                                  <  2:472055/151‑1 (MQ=255)
 aCCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCgg                                                                                                                                                 >  2:505317/1‑151 (MQ=255)
     gTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAAt                                                                                                                                              >  1:1134187/1‑150 (MQ=255)
     gTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATa                                                                                                                                             <  2:328500/151‑1 (MQ=255)
       tctgtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAgg                                                                                                                                           >  1:486547/1‑151 (MQ=255)
        ctgtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAggg                                                                                                                                          >  1:1176310/1‑151 (MQ=255)
           tttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTcc                                                                                                                                       >  1:788626/1‑151 (MQ=255)
            ttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCg                                                                                                                                      >  2:91910/1‑151 (MQ=255)
              ctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGt                                                                                                                                    >  1:304587/1‑151 (MQ=255)
                gtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGtt                                                                                                                                   <  2:1194563/150‑1 (MQ=255)
                gtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTc                                                                                                                                  <  2:621712/151‑1 (MQ=255)
                  cAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCg                                                                                                                                 <  1:614225/150‑1 (MQ=255)
                    gCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAg                                                                                                                              <  2:613075/151‑1 (MQ=255)
                     cTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGc                                                                                                                             >  1:35201/1‑151 (MQ=255)
                         cGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGAc                                                                                                                         <  1:325052/151‑1 (MQ=255)
                          gCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACt                                                                                                                        >  1:284923/1‑151 (MQ=255)
                          gCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACt                                                                                                                        >  1:1153500/1‑151 (MQ=255)
                           cAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACt                                                                                                                        >  2:859237/1‑150 (MQ=255)
                            aaaaCAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTaa                                                                                                                      <  2:724552/151‑1 (MQ=255)
                                    ccGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAgg                                                                                                              <  2:390975/151‑1 (MQ=255)
                                    ccGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAgg                                                                                                              <  1:812210/151‑1 (MQ=255)
                                     cGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGt                                                                                                             >  2:11639/1‑151 (MQ=255)
                                      gACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTc                                                                                                            >  1:31778/1‑151 (MQ=255)
                                      gACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTc                                                                                                            >  1:399831/1‑151 (MQ=255)
                                      gACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTc                                                                                                            >  2:165671/1‑151 (MQ=255)
                                         aGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTg                                                                                                         <  1:1125593/151‑1 (MQ=255)
                                           cTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGc                                                                                                        >  1:489296/1‑150 (MQ=255)
                                            tCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCAt                                                                                                      >  1:1184207/1‑151 (MQ=255)
                                             cAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCAt                                                                                                      <  2:1015293/150‑1 (MQ=255)
                                                 gCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAAc                                                                                                  >  1:977309/1‑150 (MQ=255)
                                                  cTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAAc                                                                                                  >  2:1156929/1‑149 (MQ=255)
                                                   tGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGt                                                                                               >  1:455767/1‑151 (MQ=255)
                                                       tCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcg                                                                                           >  1:566716/1‑151 (MQ=255)
                                                       tCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTc                                                                                            <  1:1205410/150‑1 (MQ=255)
                                                        ccTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcg                                                                                           <  2:164626/150‑1 (MQ=255)
                                                        ccTCTTCTTTACGCTGAACAGTGAGGATTTCGACAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcg                                                                                           <  1:678825/150‑1 (MQ=255)
                                                          tcttctTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcgccg                                                                                        >  1:1083923/1‑151 (MQ=255)
                                                           cttcttTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcgccgc                                                                                       >  1:1170361/1‑151 (MQ=255)
                                                            ttcttTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCcgcg                                                                                      >  2:419292/1‑151 (MQ=255)
                                                               tttACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCt                                                                                    >  1:734748/1‑150 (MQ=255)
                                                                         aCAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATaaa                                                                         <  2:948609/151‑1 (MQ=255)
                                                                             tGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTa                                                                      >  2:1039449/1‑150 (MQ=255)
                                                                               aGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACt                                                                   >  2:372915/1‑151 (MQ=255)
                                                                                ggATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTg                                                                  <  2:621259/151‑1 (MQ=255)
                                                                                ggATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTg                                                                  <  2:215432/151‑1 (MQ=255)
                                                                                         taataTCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGGGGTGGCGGTATAGGGGCCGGTTCGTAGCTGGCTAAACGTGAGGTCTTGCATCAACAGGTACGCCTTGCTATTGGCATAAATTAACTGTTCCTCa                                                           >  2:263507/1‑149 (MQ=255)
                                                                                           atatCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAgg                                                       <  1:508888/151‑1 (MQ=255)
                                                                                             atCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAggg                                                      <  1:598799/150‑1 (MQ=255)
                                                                                               cGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCg                                                    >  1:268345/1‑150 (MQ=255)
                                                                                                gtgtTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGaa                                                  <  2:601828/151‑1 (MQ=255)
                                                                                                 tgtTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGaa                                                  >  2:488390/1‑150 (MQ=255)
                                                                                                  gtTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGaaa                                                 <  1:175413/150‑1 (MQ=255)
                                                                                                     ggTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGAt                                             >  1:1220277/1‑151 (MQ=255)
                                                                                                          tGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAg                                         <  1:121414/150‑1 (MQ=255)
                                                                                                            cAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGAc                                       >  1:463077/1‑150 (MQ=255)
                                                                                                            cAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGAc                                       >  1:838450/1‑150 (MQ=255)
                                                                                                               aTCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAgg                                    >  2:842436/1‑150 (MQ=255)
                                                                                                               aTCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGa                                   <  1:505338/151‑1 (MQ=255)
                                                                                                                 cGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACt                                 >  1:43173/1‑151 (MQ=255)
                                                                                                                    agagGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGt                              <  2:802745/151‑1 (MQ=255)
                                                                                                                     gagGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTg                             <  2:216237/151‑1 (MQ=255)
                                                                                                                       ggTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAg                           <  2:1134209/151‑1 (MQ=255)
                                                                                                                       ggTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAg                           <  1:11641/151‑1 (MQ=255)
                                                                                                                          aTTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTaa                        <  1:1166224/151‑1 (MQ=255)
                                                                                                                          aTTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTa                         >  2:617564/1‑150 (MQ=255)
                                                                                                                             tCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTc                     >  2:932959/1‑151 (MQ=255)
                                                                                                                              cGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTcc                    >  2:570736/1‑151 (MQ=255)
                                                                                                                              cGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTc                     <  2:563683/150‑1 (MQ=255)
                                                                                                                                        ttttGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTg           >  2:1122685/1‑150 (MQ=255)
                                                                                                                                          ttGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGaa         <  2:1184221/150‑1 (MQ=255)
                                                                                                                                             tgAGGAGGCGGAATATTTTCCGGTTCGTAGCTGACTAAACGTGAGTCCTTGCATCAACAGTTTCGCCGCGCTATTGGCATACATTAACTGTTCCTCGAAGGGCGAAACGATCCAGACAGGAGTGGAGAGAAAGTCCGGGGTATTGAAGtt      >  2:49780/1‑150 (MQ=255)
                                                                                                                                               aggaggCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtgtg   >  1:1181232/1‑151 (MQ=255)
                                                                                                                                                ggagGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGAGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc  <  2:284919/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGC  >  NC_000913/1407493‑1407787

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