Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,419,071 T→G 7.3% pseudogene (30/66 nt) yhdW → pseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,419,0710TG 7.3% 285.0 / 11.7 82pseudogene (30/66 nt)yhdWpseudogene, ABC transporter periplasmic binding protein family;putative transport; Not classified; putative periplasmic binding transport protein; periplasmic‑binding component of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base T (44/32);  new base G (4/2);  total (48/34)
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.98e-01

GAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAGGTA  >  NC_000913/3418928‑3419218
                                                                                                                                               |                                                                                                                                                   
gAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTgccgc                                                                                                                                                >  3:496612/1‑149 (MQ=255)
                    gCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGGCATTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGTAAGCATTATGATAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCCTGCCGTTg                                                                                                                            >  4:556096/1‑149 (MQ=255)
                     cACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTg                                                                                                                            >  2:892187/1‑148 (MQ=255)
                      actactTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCa                                                                                                                          >  4:35430/1‑149 (MQ=255)
                          ctTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAAt                                                                                                                       <  1:314293/148‑1 (MQ=255)
                           ttttGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACGCTGGCGGCAGCCAGtgtg                                                                                                                                         <  3:481368/129‑1 (MQ=255)
                           ttttGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACGCTGGCGGCAGCCAGtgtg                                                                                                                                         >  4:481368/1‑129 (MQ=255)
                             ttGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAgg                                                                                                                   >  4:798563/1‑149 (MQ=255)
                                 ttATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGgcgc                                                                                                                <  2:537182/148‑1 (MQ=255)
                                     aaGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGctgctg                                                                                                           >  4:598117/1‑149 (MQ=255)
                                     aaGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGctgctg                                                                                                           >  4:20195/1‑149 (MQ=255)
                                     aaGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGctgctg                                                                                                           >  3:710847/1‑149 (MQ=255)
                                      agagGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGctgctg                                                                                                           >  2:853565/1‑148 (MQ=255)
                                      agagGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGctgctg                                                                                                           >  1:40559/1‑148 (MQ=255)
                                                 acagacagGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCtt                                                                                                >  2:462078/1‑148 (MQ=255)
                                                  cagacagGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTg                                                                                               <  1:853565/148‑1 (MQ=255)
                                                  cagacagGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTg                                                                                               <  2:40559/148‑1 (MQ=255)
                                                  cagacagGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGa                                                                                              <  4:710847/149‑1 (MQ=255)
                                                          tAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTc                                                                                      <  3:370387/149‑1 (MQ=255)
                                                           aaaaGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTc                                                                                      <  2:480769/148‑1 (MQ=255)
                                                           aaaaGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCa                                                                                     <  3:666862/149‑1 (MQ=255)
                                                             aaGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGgc                                                                                                         >  3:817077/1‑127 (MQ=255)
                                                             aaGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGgc                                                                                                         <  4:817077/127‑1 (MQ=255)
                                                              aGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGtgc                                                                                                                                        <  4:51156/95‑1 (MQ=255)
                                                              aGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGtgc                                                                                                                                        >  3:51156/1‑95 (MQ=255)
                                                               gACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGaa                                                                                  >  1:566415/1‑148 (MQ=255)
                                                                     tCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGt                                                                           <  3:889059/149‑1 (MQ=255)
                                                                        cAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGttt                                                                         <  2:566415/148‑1 (MQ=255)
                                                                        cAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGttt                                                                         >  1:312037/1‑148 (MQ=255)
                                                                                 tataACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAg                                                                             <  3:633733/135‑1 (MQ=255)
                                                                                 tataACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAg                                                                             >  4:633733/1‑135 (MQ=255)
                                                                                 tataACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGc                                                               >  4:286659/1‑149 (MQ=255)
                                                                                      cGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGa                                                           <  1:892187/148‑1 (MQ=255)
                                                                                       gATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATc                                                         >  3:705634/1‑149 (MQ=255)
                                                                                        aTAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATc                                                         >  2:739210/1‑148 (MQ=255)
                                                                                           aTTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGt                                                      <  1:739210/148‑1 (MQ=255)
                                                                                            ttttCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGa                                                    <  3:20195/149‑1 (MQ=255)
                                                                                            ttttCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGa                                                    <  3:286659/149‑1 (MQ=255)
                                                                                              ttCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGAAGCTGGATGCCGTTCAGAAAAAAGGTTTTGTGCAATTCGGGATCAGTGaag                                                  <  4:859688/149‑3 (MQ=255)
                                                                                                 gCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGAt                                               >  4:227375/1‑149 (MQ=255)
                                                                                                   cTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTa                                             >  4:509366/1‑149 (MQ=255)
                                                                                                   cTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTa                                             >  3:678456/1‑149 (MQ=255)
                                                                                                    tCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTAc                                            >  3:372596/1‑149 (MQ=255)
                                                                                                    tCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTAc                                            >  4:394905/1‑149 (MQ=255)
                                                                                                     cacaGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTAcc                                           >  3:163841/1‑149 (MQ=255)
                                                                                                     cacaGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTAcc                                           <  3:598117/149‑1 (MQ=255)
                                                                                                     cacaGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTAc                                            >  1:492473/1‑148 (MQ=255)
                                                                                                       caGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTg                                         >  3:217242/1‑149 (MQ=255)
                                                                                                        aGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTgg                                        <  4:749169/149‑1 (MQ=255)
                                                                                                         ggAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTgg                                        <  2:274936/148‑1 (MQ=255)
                                                                                                          gAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGt                                      >  3:868353/1‑149 (MQ=255)
                                                                                                           aaGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCg                                                              >  1:335721/1‑124 (MQ=255)
                                                                                                           aaGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCg                                                              <  2:335721/124‑1 (MQ=255)
                                                                                                           aaGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGtt                                     <  3:227375/149‑1 (MQ=255)
                                                                                                           aaGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGt                                      >  1:441085/1‑148 (MQ=255)
                                                                                                            aGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCtt                                                                                                <  4:516728/89‑1 (MQ=255)
                                                                                                            aGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCtt                                                                                                >  3:516728/1‑89 (MQ=255)
                                                                                                            aGCATTATGAAAAAGATGAGGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTtc                                    >  3:715893/1‑149 (MQ=255)
                                                                                                             gCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTtc                                    >  2:202252/1‑148 (MQ=255)
                                                                                                                ttATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTctct                                 <  2:312037/148‑1 (MQ=255)
                                                                                                                   tGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATg                             >  4:788505/1‑149 (MQ=255)
                                                                                                                       aaaGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCg                          >  2:509028/1‑148 (MQ=255)
                                                                                                                        aaGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGGCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGAt                        <  4:496612/149‑1 (MQ=255)
                                                                                                                         aGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATg                       >  4:691654/1‑149 (MQ=255)
                                                                                                                          gatgatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGaaaaaa                                                                              <  1:384873/93‑1 (MQ=255)
                                                                                                                          gatgatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGaaaaaa                                                                              >  2:384873/1‑93 (MQ=255)
                                                                                                                          gatgatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGc                      <  3:788505/149‑1 (MQ=255)
                                                                                                                            tgatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTg                    >  3:668544/1‑149 (MQ=255)
                                                                                                                             gatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAg                                                       >  1:856005/1‑113 (MQ=255)
                                                                                                                             gatgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAg                                                       <  2:856005/113‑1 (MQ=255)
                                                                                                                              atgatAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGa                   <  1:509028/148‑1 (MQ=255)
                                                                                                                                 atAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACgg                >  1:29821/1‑148 (MQ=255)
                                                                                                                                   aGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGTTTGATGCCGTTCAGAAAAACGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATCCTGACGGTaa             <  3:556096/149‑1 (MQ=255)
                                                                                                                                      cacnCTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGt           >  1:290140/1‑148 (MQ=255)
                                                                                                                                      cacgCTGGCGGCAGCCAGTGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGtt          >  4:22326/1‑149 (MQ=255)
                                                                                                                                      cacgCTGGCGGCAGCCAGTGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGtt          >  4:647121/1‑149 (MQ=255)
                                                                                                                                      cacgCTGGCGGCAGCCAGTGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGt           >  1:524246/1‑148 (MQ=255)
                                                                                                                                      cacaCTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTa              >  2:828308/1‑145 (MQ=255)
                                                                                                                                      cacaCTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTa              <  1:828308/145‑1 (MQ=255)
                                                                                                                                      cacaCTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGt           >  2:101331/1‑148 (MQ=255)
                                                                                                                                       acacTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGttt         >  3:344853/1‑149 (MQ=255)
                                                                                                                                           tGGCGGCAGCCAGTGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAg     <  3:22326/149‑1 (MQ=255)
                                                                                                                                               tGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAGGTa  >  1:473097/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
GAAATCGTGCGCAAAAGCTGGCACTACTTTTGCTTATAAGAGGATGGCCACAGACAGGTAAAAGACGTTTCCCAAACGTCCTATAACGATAATTTTCGCCTCACAGGAAGCATTATGAAAAAGATGATGATAGCCACACTGGCTGCCGCCAGCGTGCTGCTTGCCGTTGCAAATCAGGCGCTGCTGGCGCGACGCTTGATGCCGTTCAGAAAAAAGGTTTTGTGCAATGCGGGATCAGTGATGGATTACCTGGGTTCTCTTATGCCGATGCTGACGGTAAGTTTTCAGGTA  >  NC_000913/3418928‑3419218

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