Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,015,192 A→G 6.0% S105G (AGC→GGC)  fliG → flagellar motor switching and energizing component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,015,1920AG 6.0% 273.9 / 5.7 83S105G (AGC→GGC) fliGflagellar motor switching and energizing component
Reads supporting (aligned to +/- strand):  ref base A (38/40);  new base G (3/2);  total (41/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGA  >  NC_000913/2015049‑2015331
                                                                                                                                               |                                                                                                                                           
cGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCgg                                                                                                                                         >  2:515018/1‑148 (MQ=255)
cGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGt                                                                                                                                        <  4:179488/149‑1 (MQ=255)
       tGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCTGGCGGTAtcgcga                                                                                                                                 <  4:273929/149‑6 (MQ=255)
        gAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAAc                                                                                                                                <  4:223062/149‑1 (MQ=255)
         aGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAAc                                                                                                                                >  2:33372/1‑148 (MQ=255)
               aagCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCaa                                                                                                                          >  2:170485/1‑148 (MQ=255)
                   tGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTa                                                                                                                     <  4:57290/149‑1 (MQ=255)
                    gAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGTGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTAt                                                                                                                    >  4:306130/1‑149 (MQ=255)
                    gAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTAt                                                                                                                    >  4:660780/1‑149 (MQ=255)
                        aGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAg                                                                                                                <  4:131743/149‑1 (MQ=255)
                            tGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCAAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCcac                                                                                                            <  4:260848/149‑1 (MQ=255)
                                 cACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACagag                                                                                                        <  1:243607/148‑1 (MQ=255)
                                  aCTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAgcg                                                                                                      >  4:675619/1‑149 (MQ=255)
                                    tGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCTGGCGGTAtcg                                                                                                                                    >  4:317911/1‑115 (MQ=255)
                                    tGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCTGGCGGTAtcg                                                                                                                                    <  3:317911/117‑3 (MQ=255)
                                    tGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGcgcc                                                                                                    >  3:241379/1‑149 (MQ=255)
                                     gAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGc                                                                                                                                           <  1:164046/109‑1 (MQ=255)
                                     gAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGc                                                                                                                                           >  2:164046/1‑109 (MQ=255)
                                     gAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGcgcc                                                                                                    >  2:216491/1‑148 (MQ=255)
                                       atatCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCgccgcc                                                                                                 >  3:876364/1‑149 (MQ=255)
                                         atCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGccgccg                                                                                                <  1:170485/148‑1 (MQ=255)
                                          tCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGAt                                                                                              <  3:704850/149‑1 (MQ=255)
                                          tCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCGCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGAt                                                                                              >  4:523921/1‑149 (MQ=255)
                                           cAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGAt                                                                                              >  2:301864/1‑148 (MQ=255)
                                                ccAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGAtt                                                                                        >  4:695535/1‑149 (MQ=255)
                                                  aaCGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGt                                                                                                                                        >  1:490353/1‑99 (MQ=255)
                                                  aaCGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAGCGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTc                                                                                       >  2:96315/1‑148 (MQ=255)
                                                  aaCGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCAAGCGGt                                                                                                                                        <  2:490353/99‑1 (MQ=255)
                                                   aCGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTcgc                                                                                     <  4:540956/149‑1 (MQ=255)
                                                     gATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGa                                                                                   <  3:479732/149‑1 (MQ=255)
                                                      aTTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGa                                                                                   <  1:69706/148‑1 (MQ=255)
                                                        tATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCaa                                                                                                                          >  1:690648/1‑107 (MQ=255)
                                                        tATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCaa                                                                                                                          <  2:690648/107‑1 (MQ=255)
                                                        tATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGa                                                                                <  3:801266/149‑1 (MQ=255)
                                                        tATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGa                                                                                <  3:745328/149‑1 (MQ=255)
                                                        tATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGa                                                                                <  3:695535/149‑1 (MQ=255)
                                                                       ttGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTa                                                                  <  1:301864/148‑1 (MQ=255)
                                                                        tGGTCAAAGCTCTGTGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATc                                                                <  3:306130/149‑1 (MQ=255)
                                                                            cAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGcc                                                             >  1:509882/1‑148 (MQ=255)
                                                                             aaaGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGccac                                                           >  3:123034/1‑149 (MQ=255)
                                                                              aaGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGcgccgc                                                                                                  >  2:327764/1‑109 (MQ=255)
                                                                              aaGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGcgccgc                                                                                                  <  1:327764/109‑1 (MQ=255)
                                                                                  tctGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTc                                                      >  4:500059/1‑149 (MQ=255)
                                                                                        tGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTg                                                 >  2:473942/1‑148 (MQ=255)
                                                                                        tGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCAATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGc                                                <  4:666520/149‑1 (MQ=255)
                                                                                         gaagaaCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGc                                                >  1:74191/1‑148 (MQ=255)
                                                                                           agaaCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCAt                                              >  1:692542/1‑148 (MQ=255)
                                                                                           agaaCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATc                                             >  4:349456/1‑149 (MQ=255)
                                                                                            gaaCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCt                                            <  3:349342/149‑1 (MQ=255)
                                                                                            gaaCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCt                                            >  4:422074/1‑149 (MQ=255)
                                                                                              aCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGcc                                                             <  1:871080/130‑1 (MQ=255)
                                                                                              aCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGcc                                                             >  2:871080/1‑130 (MQ=255)
                                                                                              aCGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGa                                          <  3:660780/149‑1 (MQ=255)
                                                                                               cGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGa                                          <  2:100711/148‑1 (MQ=255)
                                                                                                     gccAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAgcgcgc                                   >  4:285278/1‑149 (MQ=255)
                                                                                                     gccAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAgcgcg                                    <  2:74191/148‑1 (MQ=255)
                                                                                                          ccTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCaa                               <  1:473942/148‑1 (MQ=255)
                                                                                                                ggAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATc                                                                >  4:666999/1‑109 (MQ=255)
                                                                                                                ggAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATc                                                                <  3:666999/109‑1 (MQ=255)
                                                                                                                ggAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCTCCCAAGccgccg                        <  4:123034/149‑1 (MQ=255)
                                                                                                                ggAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAgccgcc                         >  1:97433/1‑148 (MQ=255)
                                                                                                                 gAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGa                       >  3:329771/1‑149 (MQ=255)
                                                                                                                      tatTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTc                  >  3:345806/1‑149 (MQ=255)
                                                                                                                       atTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCt                 >  4:716372/1‑149 (MQ=255)
                                                                                                                        ttCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTTAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTg                >  4:524761/1‑149 (MQ=255)
                                                                                                                         tctcGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTg                >  2:193100/1‑148 (MQ=255)
                                                                                                                           tcGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCg             >  4:783029/1‑149 (MQ=255)
                                                                                                                            cGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCg             >  1:312467/1‑148 (MQ=255)
                                                                                                                            cGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCg             >  2:544008/1‑148 (MQ=255)
                                                                                                                             gAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCgtt           <  4:345806/149‑1 (MQ=255)
                                                                                                                             gAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCgtt           <  3:349456/149‑1 (MQ=255)
                                                                                                                                aCTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCgttgt         <  1:515018/148‑1 (MQ=255)
                                                                                                                                aCTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCgttgt         <  2:97433/148‑1 (MQ=255)
                                                                                                                                 cTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCgttgtt        <  2:321156/148‑1 (MQ=255)
                                                                                                                                  tCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCg      >  4:391576/1‑149 (MQ=255)
                                                                                                                                   ttcaaTACCGCTGGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCg      >  2:58646/5‑148 (MQ=255)
                                                                                                                                   cgcgATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGa     <  3:500059/149‑1 (MQ=255)
                                                                                                                                   cccgATACCGCCAGCGGTATTGAAACGCTCAACTGTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGa     <  3:739879/147‑1 (MQ=255)
                                                                                                                                    gcgATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGa     <  1:721600/148‑1 (MQ=255)
                                                                                                                                    gcgATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGAt    <  4:196409/149‑1 (MQ=255)
                                                                                                                                      aaTACCGCTGGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGa  >  4:614399/2‑149 (MQ=255)
                                                                                                                                       aTACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGa  <  1:216491/148‑1 (MQ=255)
                                                                                                                                       aTACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGa  <  2:312467/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
CGGAGTTTGAGCAAGAAGCTGAACAGTTTGCCGCACTGAATATCAACGCCAACGATTATCTGCGCTCGGTATTGGTCAAAGCTCTGGGTGAAGAACGTGCCGCCAGCCTGCTGGAAGATATTCTCGAAACTCGCGATACCGCCAGCGGTATTGAAACGCTCAACTTTATGGAGCCACAGAGCGCCGCCGATCTGATTCGCGATGAGCATCCGCAAATTATCGCCACCATTCTGGTGCATCTGAAGCGCGCCCAAGCCGCCGATATTCTGGCGTTGTTCGATGA  >  NC_000913/2015049‑2015331

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