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,6410AG99.1% 411.2 / ‑4.8 112V52A (GTC→GCC) ydaMdiguanylate cyclase, csgD regulator
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (55/56);  total (55/57)
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.87e-01

Only 100 of 112 total aligned reads displayed.

CACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGC  >  NC_000913/1407493‑1407787
                                                                                                                                                    |                                                                                                                                                  
cACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCg                                                                                                                                                  <  1:208890/151‑1 (MQ=255)
cACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCg                                                                                                                                                  >  2:266297/1‑151 (MQ=255)
   cGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGa                                                                                                                                                >  4:187183/1‑150 (MQ=255)
    ggTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAAt                                                                                                                                              >  2:913810/1‑151 (MQ=255)
    ggTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAAt                                                                                                                                              <  3:1209069/151‑1 (MQ=255)
     ggtctttttctgtCAGCTTTCTCAAAACAAGCCGACAGCTCAATGCGGGTTCCTCTGCTTTACGATGAACAGTGAGTATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTGCGGCAGTTCTTTTTGTGAGGAGACGGaa                                                                                                                                               <  4:417619/147‑1 (MQ=255)
     gTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAAt                                                                                                                                              >  3:508656/1‑150 (MQ=255)
      ttctgtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAAt                                                                                                                                              >  3:1358619/1‑149 (MQ=255)
      ttctgtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAAt                                                                                                                                              >  2:744054/1‑149 (MQ=255)
       tctgtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAgg                                                                                                                                           >  3:1259951/1‑151 (MQ=255)
       tctgtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAg                                                                                                                                            <  2:657992/150‑1 (MQ=255)
          gtttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTc                                                                                                                                        >  3:353829/1‑151 (MQ=255)
            ttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTcc                                                                                                                                       <  3:1234792/150‑1 (MQ=255)
            ttctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCg                                                                                                                                      >  3:878156/1‑151 (MQ=255)
              ctgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGt                                                                                                                                    >  2:1045512/1‑151 (MQ=255)
               tgtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGtt                                                                                                                                   >  3:128461/1‑151 (MQ=255)
                gtCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGt                                                                                                                                    <  3:937515/149‑1 (MQ=255)
                 tCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTc                                                                                                                                  >  2:769229/1‑150 (MQ=255)
                 tCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCg                                                                                                                                 >  2:1109672/1‑151 (MQ=255)
                             aaaCAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTaa                                                                                                                      >  3:974911/1‑150 (MQ=255)
                               aCAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACg                                                                                                                   >  2:1398337/1‑151 (MQ=255)
                               aCAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATCTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACg                                                                                                                   >  2:724278/1‑151 (MQ=255)
                                cAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGt                                                                                                                  <  2:519994/151‑1 (MQ=255)
                                   gCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTg                                                                                                                 >  2:492477/1‑149 (MQ=255)
                                     cGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGt                                                                                                             >  2:1060117/1‑151 (MQ=255)
                                      gACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTc                                                                                                            <  4:109402/151‑1 (MQ=255)
                                      gACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTc                                                                                                            >  2:281931/1‑151 (MQ=255)
                                        cAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCtt                                                                                                          <  3:266323/151‑1 (MQ=255)
                                        cAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCtt                                                                                                          >  3:1140923/1‑151 (MQ=255)
                                        cAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCt                                                                                                           <  3:1398408/150‑1 (MQ=255)
                                         aGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCtt                                                                                                          <  3:799534/150‑1 (MQ=255)
                                         aGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTg                                                                                                         >  2:998533/1‑151 (MQ=255)
                                         aGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTg                                                                                                         <  4:153977/151‑1 (MQ=255)
                                              aaTGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATc                                                                                                     >  3:757131/1‑150 (MQ=255)
                                              aaTGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCa                                                                                                    <  3:1238091/151‑1 (MQ=255)
                                               aTGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCa                                                                                                    >  3:805510/1‑150 (MQ=255)
                                               aTGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCa                                                                                                    >  2:553854/1‑150 (MQ=255)
                                                 gCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAAc                                                                                                  >  3:542525/1‑150 (MQ=255)
                                                  cTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAAc                                                                                                  >  2:876889/1‑149 (MQ=255)
                                                    gTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAg                                                                                                <  2:708773/149‑1 (MQ=255)
                                                    gTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGtt                                                                                              <  1:20757/151‑1 (MQ=255)
                                                    gTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGtt                                                                                              <  4:52386/151‑1 (MQ=255)
                                                      ttCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGttt                                                                                             >  3:951837/1‑150 (MQ=255)
                                                       tCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTc                                                                                            <  3:968342/150‑1 (MQ=255)
                                                        ccTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcgc                                                                                          >  3:385372/1‑151 (MQ=255)
                                                         ctctTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcgcc                                                                                         >  3:211816/1‑151 (MQ=255)
                                                           cttcttTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTcgccgc                                                                                       >  3:108403/1‑151 (MQ=255)
                                                             tcttTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCcgcg                                                                                      <  4:159159/150‑1 (MQ=255)
                                                                ttACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTa                                                                                   <  2:974870/150‑1 (MQ=255)
                                                                ttACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTa                                                                                   <  4:364914/150‑1 (MQ=255)
                                                                 tACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTAtt                                                                                 >  1:265244/1‑151 (MQ=255)
                                                                  aCGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTAtt                                                                                 <  3:744086/150‑1 (MQ=255)
                                                                  aCGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTg                                                                                >  3:936517/1‑151 (MQ=255)
                                                                      tGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCAt                                                                            >  2:1090348/1‑151 (MQ=255)
                                                                          cAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAAt                                                                        <  2:951800/151‑1 (MQ=255)
                                                                            gTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTa                                                                      <  1:47380/151‑1 (MQ=255)
                                                                              gAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTaa                                                                     <  3:307518/150‑1 (MQ=255)
                                                                              gAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAAc                                                                    >  2:496474/1‑151 (MQ=255)
                                                                                ggATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTg                                                                  <  2:353800/151‑1 (MQ=255)
                                                                                ggATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTg                                                                  <  3:1109725/151‑1 (MQ=255)
                                                                                 gATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGt                                                                 <  2:128445/151‑1 (MQ=255)
                                                                                      cGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCt                                                             <  2:622171/150‑1 (MQ=255)
                                                                                      cGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTc                                                            >  3:477128/1‑151 (MQ=255)
                                                                                      cGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTc                                                            >  4:352478/1‑151 (MQ=255)
                                                                                      cGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTc                                                            >  1:35597/1‑151 (MQ=255)
                                                                                       gATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCa                                                           <  1:8857/151‑1 (MQ=255)
                                                                                       gATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCa                                                           <  4:407339/151‑1 (MQ=255)
                                                                                          aataTCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCaaa                                                         <  2:1391585/150‑1 (MQ=255)
                                                                                          aataTCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCaa                                                          >  2:414360/1‑149 (MQ=255)
                                                                                           atatCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAAc                                                                                                                    <  2:396237/90‑1 (MQ=255)
                                                                                           atatCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAAc                                                                                                                    >  3:396269/1‑90 (MQ=255)
                                                                                             atCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGc                                                     <  3:1203538/151‑1 (MQ=255)
                                                                                               cGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCg                                                    <  4:240900/150‑1 (MQ=255)
                                                                                               cGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGa                                                   <  3:496511/151‑1 (MQ=255)
                                                                                                  gtTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAAc                                                <  2:284365/151‑1 (MQ=255)
                                                                                                   ttGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACg                                               >  1:44489/1‑151 (MQ=255)
                                                                                                      gTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATc                                            <  2:385336/151‑1 (MQ=255)
                                                                                                            cAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACa                                      <  2:757097/151‑1 (MQ=255)
                                                                                                               aTCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAgg                                    >  4:435146/1‑150 (MQ=255)
                                                                                                                 cGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGAc                                  <  2:878121/150‑1 (MQ=255)
                                                                                                                 cGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGAc                                  <  3:377385/150‑1 (MQ=255)
                                                                                                                  ggAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTg                                <  3:1019396/151‑1 (MQ=255)
                                                                                                                      agGTATTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGa                            <  2:936479/151‑1 (MQ=255)
                                                                                                                          aTTTCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTa                         >  2:1028350/1‑150 (MQ=255)
                                                                                                                            ttCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAg                       <  2:1358549/150‑1 (MQ=255)
                                                                                                                            ttCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGt                      >  1:407835/1‑151 (MQ=255)
                                                                                                                             tCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGt                      <  3:1090396/150‑1 (MQ=255)
                                                                                                                             tCGGCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTc                     >  1:96833/1‑151 (MQ=255)
                                                                                                                                gCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTcc                    >  3:896102/1‑149 (MQ=255)
                                                                                                                                gCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTcc                    >  4:299681/1‑149 (MQ=255)
                                                                                                                                gCAGTTCTTTTTGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAg                  <  1:329618/151‑1 (MQ=255)
                                                                                                                                       tttttGTGAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTg           >  1:232310/1‑151 (MQ=255)
                                                                                                                                         tttGTGAGGAGGCGGAATAGTGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCCAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGAGTAt             <  4:153681/147‑1 (MQ=255)
                                                                                                                                            gtgAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGtt      <  3:1045558/151‑1 (MQ=255)
                                                                                                                                            gtgAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGt       <  3:281958/150‑1 (MQ=255)
                                                                                                                                             tgAGGAGGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGtt      <  2:692709/150‑1 (MQ=255)
                                                                                                                                              gaggagGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtgt    <  2:1259886/151‑1 (MQ=255)
                                                                                                                                              gaggagGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTtgt    <  2:108387/151‑1 (MQ=255)
                                                                                                                                                ggagGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc  >  2:489113/1‑151 (MQ=255)
                                                                                                                                                ggagGCGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGc  >  3:734674/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CACCGGTTCTGTTTCTGTCAGCTTTCGCAAAACAAGCCGACAGCTCAATGCTGTTTCCTCTTCTTTACGCTGAACAGTGAGGATTTCGATAATATCGTGTTGGTTTTGCAGATCGGAGAGGTATTTCGGCAGTTCTTTTTGTGAGGAGACGGAATAGGGTCCGGTTCGTAGCTGACTAAACGTGAGGTCTTGCATCAACAGTTTCGCCGCGCTATTGGCATAAATTAACTGTTCCTCAAAGGGCGAAACGATCCAGACAGGACTGGTGAGTAAGTCCAGGGTATTGAAGTTGTGC  >  NC_000913/1407493‑1407787

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