Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,364,995 T→G 10.0% intergenic (+16/+23) nudI → / ← ais nucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,364,9950TG10.0% 426.7 / 35.7 137intergenic (+16/+23)nudI/aisnucleoside triphosphatase/putative LPS core heptose(II)‑phosphate phosphatase
Reads supporting (aligned to +/- strand):  ref base T (75/48);  new base G (8/6);  total (83/54)
Fisher's exact test for biased strand distribution p-value = 7.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.39e-01

Only 100 of 138 total aligned reads displayed.

TGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATG  >  NC_000913/2364853‑2365143
                                                                                                                                              |                                                                                                                                                    
 gCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGc                                                                                                                                               >  3:820450/1‑149 (MQ=255)
  ccAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCt                                                                                                                                              >  4:592547/1‑149 (MQ=255)
    aaCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGAAAGCCATTACTGGCTGt                                                                                                                                            >  1:48923/1‑149 (MQ=255)
     aCCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTAGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGTc                                                                                                                                           <  3:168915/149‑1 (MQ=255)
      ccGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCa                                                                                                                                          <  3:592547/149‑1 (MQ=255)
       cGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAg                                                                                                                                         >  3:94182/1‑149 (MQ=255)
         agaAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGAt                                                                                                                                       <  1:752726/149‑1 (MQ=255)
             gTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaa                                                                                                                                   >  1:282680/1‑149 (MQ=255)
                 aaaTAAACGAAGAGTTTCAGGACTACGCGCGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAAtt                                                                                                                               >  4:61744/1‑149 (MQ=255)
                   aTAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTa                                                                                                                             <  1:995670/149‑1 (MQ=255)
                   aTAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTa                                                                                                                             <  2:608401/149‑1 (MQ=255)
                     aaaCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGt                                                                                                                           <  3:211486/149‑1 (MQ=255)
                      aaCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTg                                                                                                                          >  3:40824/1‑149 (MQ=255)
                       aCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                 >  4:785504/1‑141 (MQ=255)
                       aCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCAGTAATGGCTGt                                                                                                                                            >  4:540365/1‑130 (MQ=255)
                       aCGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATCATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGaaaa                                                                                                                                 <  3:785504/141‑1 (MQ=255)
                          aaGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTa                                                                                                                      >  1:829325/1‑149 (MQ=255)
                              gTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGa                                                                                                                  >  1:935378/1‑149 (MQ=255)
                                 tCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAAtt                                                                                                               >  3:699106/1‑149 (MQ=255)
                                           gcgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAg                                                                                                     >  4:558380/1‑149 (MQ=255)
                                            cgTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGa                                                                                                    <  2:282680/149‑1 (MQ=255)
                                                ggTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATaaa                                                                                                <  2:476040/149‑1 (MQ=255)
                                                 gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                               >  3:84947/1‑149 (MQ=255)
                                                 gTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAAc                                                                                               >  4:226283/1‑149 (MQ=255)
                                                  tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACt                                                                                              >  1:923321/1‑149 (MQ=255)
                                                  tAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACt                                                                                              >  2:1076728/1‑149 (MQ=255)
                                                       ccTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGcc                                                                                         >  1:352605/1‑149 (MQ=255)
                                                        cTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCt                                                                                        >  1:485357/1‑149 (MQ=255)
                                                         tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                                                                       <  3:416461/149‑1 (MQ=255)
                                                         tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                                                                       >  1:496482/1‑149 (MQ=255)
                                                         tGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtt                                                                                       >  2:1019841/1‑149 (MQ=255)
                                                          gAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCttt                                                                                      >  3:347761/1‑149 (MQ=255)
                                                           aaGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCtttt                                                                                     <  1:457104/149‑1 (MQ=255)
                                                              aTCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                                                                  >  1:1035988/1‑149 (MQ=255)
                                                                 tGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCcaca                                                                               <  1:1019841/149‑1 (MQ=255)
                                                                      cATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCAcatgcat                                                                          <  4:74996/149‑1 (MQ=255)
                                                                        ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGa                                                                        >  1:136071/1‑149 (MQ=255)
                                                                        ttATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGa                                                                        <  1:113192/149‑1 (MQ=255)
                                                                            gATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTaa                                                                    >  1:529456/1‑149 (MQ=255)
                                                                              tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                  >  3:394994/1‑149 (MQ=255)
                                                                              tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                  >  3:445328/1‑149 (MQ=255)
                                                                              tttGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAAc                                                                  >  3:54808/1‑149 (MQ=255)
                                                                                        gCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGATTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATa                                                        >  4:611979/1‑149 (MQ=255)
                                                                                         ccACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa                                                       >  4:670720/1‑149 (MQ=255)
                                                                                         ccACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa                                                       >  3:464626/1‑149 (MQ=255)
                                                                                         ccACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATaa                                                       <  3:313448/149‑1 (MQ=255)
                                                                                          cACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt                                                      <  1:371577/149‑1 (MQ=255)
                                                                                               gAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGt                                                 >  3:839331/1‑149 (MQ=255)
                                                                                                aaaaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt                                                >  3:691675/1‑149 (MQ=255)
                                                                                                aaaaaCGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGtt                                                <  2:388568/149‑1 (MQ=255)
                                                                                                    acgttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTaaa                                            <  2:521672/149‑1 (MQ=255)
                                                                                                    acgttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTaaa                                            <  4:324551/149‑1 (MQ=255)
                                                                                                     cgttacgttTGAAAGGTCTTTTGTATCAATGATGACAGCCAGTACTGGCTGTCAGATTGAAAATTTAGTGGTTGACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTGTAAAt                                           <  4:1018044/149‑1 (MQ=255)
                                                                                                       ttacgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGt                                         >  3:407078/1‑149 (MQ=255)
                                                                                                         acgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcg                                       <  3:670720/149‑1 (MQ=255)
                                                                                                         acgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcg                                       >  4:253873/1‑149 (MQ=255)
                                                                                                          cgttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgc                                      >  1:198391/1‑149 (MQ=255)
                                                                                                           gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                                                                 <  2:150122/41‑1 (MQ=255)
                                                                                                           gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTg                                                                                                                                                 >  1:150122/1‑41 (MQ=255)
                                                                                                           gttTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgcg                                     >  3:1018673/1‑149 (MQ=255)
                                                                                                           gttTGAAAGGCCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTcgcg                                     >  4:244346/1‑149 (MQ=255)
                                                                                                                aaaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGGCTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACg                                >  4:825459/1‑149 (MQ=255)
                                                                                                                aaaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACg                                >  3:435818/1‑149 (MQ=255)
                                                                                                                 aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGc                               >  2:159745/1‑149 (MQ=255)
                                                                                                                 aaGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGc                               >  4:781450/1‑149 (MQ=255)
                                                                                                                  aGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCt                              >  3:18811/1‑149 (MQ=255)
                                                                                                                    gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTc                                                                                                              >  3:513213/1‑67 (MQ=255)
                                                                                                                    gTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                           >  4:402972/1‑134 (MQ=255)
                                                                                                                       ttctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCt                         >  1:467285/1‑149 (MQ=255)
                                                                                                                         ctGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCttt                       <  3:253873/149‑1 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                                                                  >  2:505524/1‑89 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTcc                                                                                  <  1:505524/89‑1 (MQ=255)
                                                                                                                          tGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTa                      >  3:174524/1‑149 (MQ=255)
                                                                                                                           gTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTTCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAg                     >  1:332328/1‑149 (MQ=255)
                                                                                                                              aCAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCaa                  >  1:893285/1‑149 (MQ=255)
                                                                                                                               cAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAat                 <  4:445328/149‑1 (MQ=255)
                                                                                                                                aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                                                                          <  3:494943/59‑1 (MQ=255)
                                                                                                                                aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCa                                                                                                          >  4:494943/1‑59 (MQ=255)
                                                                                                                                aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAata                <  4:839331/149‑1 (MQ=255)
                                                                                                                                aaTGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAata                >  1:850765/1‑149 (MQ=255)
                                                                                                                                  tgatgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatata              <  2:1263/149‑1 (MQ=255)
                                                                                                                                  atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatata              >  2:1072682/4‑149 (MQ=255)
                                                                                                                                  atctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatata              >  4:566002/4‑149 (MQ=255)
                                                                                                                                   tctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat             <  3:566002/147‑1 (MQ=255)
                                                                                                                                   tctgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat             >  3:1065146/3‑149 (MQ=255)
                                                                                                                                     tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAatatat             >  1:806436/1‑147 (MQ=255)
                                                                                                                                     tgaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGt           >  2:1095881/1‑149 (MQ=255)
                                                                                                                                     tgaCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGt           >  4:71815/1‑149 (MQ=255)
                                                                                                                                      gaCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTc          >  1:288657/1‑149 (MQ=255)
                                                                                                                                       aCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTTCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCa         <  2:332328/149‑1 (MQ=255)
                                                                                                                                       aCAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGCCGCGTCACGCTTATCTTTAGCAATATATGTCa         >  2:800936/1‑149 (MQ=255)
                                                                                                                                       aCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt                                                      <  1:728261/104‑1 (MQ=255)
                                                                                                                                       aCAGCCAGTAATGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAAt                                                      >  2:728261/1‑104 (MQ=255)
                                                                                                                                         aGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAgg       <  4:401297/149‑1 (MQ=255)
                                                                                                                                          gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                           >  2:435869/1‑112 (MQ=255)
                                                                                                                                          gCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAAt                                           <  1:435869/112‑1 (MQ=255)
                                                                                                                                           ccATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCa     >  3:563913/1‑149 (MQ=255)
                                                                                                                                            cATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCaa    <  2:829325/149‑1 (MQ=255)
                                                                                                                                             aTTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAAt   >  2:953226/1‑149 (MQ=255)
                                                                                                                                              ttACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATg  >  2:816575/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
TGCCAACCGAGAAGTGAAAATAAACGAAGAGTTTCAGGACTACGCGTGGGTAAAACCTGAAGATCTGGTGCATTATGATTTGAATGTCGCCACCCGAAAAACGTTACGTTTGAAAGGTCTTCTGTAACAATGATGACAGCCATTACTGGCTGTCAGATTGAAAATTTAGTGGTTAACGAATTCCCCATCCAGATAAACTTTGCCTTTTTCCACATGCATGACTAAACCATCCAGATAATCAGGTTTAAATGTCGCGTCACGCTTATCTTTAGCAATATATGTCAGGCAATG  >  NC_000913/2364853‑2365143

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