Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,639,557 A→T 7.7% intergenic (+27/+33) creD → / ← arcA inner membrane protein/response regulator in two‑component regulatory system with ArcB or CpxA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,639,5570AT 7.7% 284.2 / 5.7 86intergenic (+27/+33)creD/arcAinner membrane protein/response regulator in two‑component regulatory system with ArcB or CpxA
Reads supporting (aligned to +/- strand):  ref base A (46/33);  new base T (4/3);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTC  >  NC_000913/4639410‑4639703
                                                                                                                                                   |                                                                                                                                                   
ggtggCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTaaaaa                                                                                                                                                    <  2:868211/149‑1 (MQ=255)
    gcgcTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAgcgc                                                                                                                                                <  1:249564/149‑1 (MQ=255)
        tGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTTTTTAGCGCCG‑tt                                                                                                                                           <  2:109760/149‑1 (MQ=255)
        tGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTTTTTAGCGCCG‑tt                                                                                                                                           <  1:473892/149‑1 (MQ=255)
          gCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttt                                                                                                                                         >  4:667808/1‑149 (MQ=255)
            cGGCAGGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAGAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAGGGTTGAAAAATAAAAACGGCGCTAAATAGCACCG‑tttttt                                                                                                                                       >  4:608603/1‑149 (MQ=255)
             ggCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑ttttttt                                                                                                                                      <  3:442417/149‑1 (MQ=255)
             ggCATGATGTTTGTGACCAGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑ttttttt                                                                                                                                      <  3:295598/149‑1 (MQ=255)
              gCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttttttt                                                                                                                                     <  4:107570/149‑1 (MQ=255)
              gCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑tttttttt                                                                                                                                     <  2:87007/149‑1 (MQ=255)
                 tgatgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAA‑AAAAACGGCGCTTTTTAGCGCCG‑tttttatt                                                                                                                                     <  2:655786/145‑1 (MQ=255)
                 tgatgTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAA‑AAAAACGGCGCTTTTTAGCGCCG‑tttttatt                                                                                                                                     >  1:655786/1‑145 (MQ=255)
                      tttGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACggtgg                                                                                                                             >  4:183287/1‑149 (MQ=255)
                         gtgACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTaa                                                                                                                          <  4:450559/149‑1 (MQ=255)
                              ccGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCg                                                                                                                     >  4:302703/1‑149 (MQ=255)
                               cGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGCAAAATAAAAACGGCGCTAAAAGGCGCCT‑TTTTTTTTGACGGTGGTAAAGCCGa                                                                                                                    <  2:805154/149‑1 (MQ=255)
                               cGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGa                                                                                                                    >  2:329042/1‑149 (MQ=255)
                                gTAATATCGACTGATATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAAGGGAGGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGAt                                                                                                                   <  2:420785/149‑1 (MQ=255)
                                    tatCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAAt                                                                                                               <  1:644748/149‑1 (MQ=255)
                                        gACTGGTATGCGTTTTCACTGCCAAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTc                                                                                                           >  4:56804/1‑149 (MQ=255)
                                         aCTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTcc                                                                                                          <  4:259462/149‑1 (MQ=255)
                                          cTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCa                                                                                                         <  2:124194/149‑1 (MQ=255)
                                              tATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATc                                                                                                     <  4:878080/149‑1 (MQ=255)
                                              tATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑‑TTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCa                                                                                                    >  3:48761/1‑149 (MQ=255)
                                                tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑ttttttt                                                                                                                                      >  3:664619/1‑114 (MQ=255)
                                                tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑ttttttt                                                                                                                                      <  4:664619/114‑1 (MQ=255)
                                                tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTGAAGCCGATTAATCTTCCAGATCAc                                                                                                   >  1:304841/1‑149 (MQ=255)
                                                tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCAc                                                                                                   >  1:384400/1‑149 (MQ=255)
                                                tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCAc                                                                                                   >  2:718095/1‑149 (MQ=255)
                                                   gTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                                >  4:149643/1‑149 (MQ=255)
                                                   gTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                                >  3:156124/1‑149 (MQ=255)
                                                   gTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                                >  2:789096/1‑149 (MQ=255)
                                                    ttttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCa                                                                                                         <  1:338516/140‑1 (MQ=255)
                                                    ttttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCa                                                                                                         >  2:338516/1‑140 (MQ=255)
                                                    ttttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCa                                                                                               <  3:667808/149‑1 (MQ=255)
                                                          cTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCg                                                                                         >  1:533211/1‑149 (MQ=255)
                                                           tGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGa                                                                                        >  4:204100/1‑149 (MQ=255)
                                                              cGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATaa                                                                                     <  2:285901/149‑1 (MQ=255)
                                                               gAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAAc                                                                                    >  2:334447/1‑149 (MQ=255)
                                                                aaaaTGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAAcc                                                                                   >  3:570419/1‑149 (MQ=255)
                                                                    tGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCa                                                                               >  1:168190/1‑149 (MQ=255)
                                                                           cAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGaa                                                                        >  1:292018/1‑149 (MQ=255)
                                                                           cAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGaa                                                                        >  3:432326/1‑149 (MQ=255)
                                                                                aagaagTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAtggtg                                                                   >  1:360701/1‑149 (MQ=255)
                                                                                aagaagTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑‑TTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAAtggtgg                                                                  <  4:48761/149‑1 (MQ=255)
                                                                                     gTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgat                                                              >  1:368789/1‑149 (MQ=255)
                                                                                     gTTACAACGGACAATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgat                                                              <  4:305047/149‑1 (MQ=255)
                                                                                         cAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAtt                                                          >  1:24441/1‑149 (MQ=255)
                                                                                          aaCGGACGCTGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatga                                                            >  3:620787/1‑146 (MQ=255)
                                                                                          aaCGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCgatga                                                            <  4:620787/146‑1 (MQ=255)
                                                                                          aaCGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAttt                                                         >  2:889152/1‑149 (MQ=255)
                                                                                          aaCGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGAttt                                                         >  3:367480/1‑149 (MQ=255)
                                                                                             ggACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccg                                                      <  4:321705/149‑1 (MQ=255)
                                                                                              gACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccgg                                                     >  1:759364/1‑149 (MQ=255)
                                                                                               aCGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggc                                                    <  1:329042/149‑1 (MQ=255)
                                                                                                  aTGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgta                                                 >  3:880354/1‑149 (MQ=255)
                                                                                                  aTGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgta                                                 >  1:516609/1‑149 (MQ=255)
                                                                                                    gAGTTACGTATCTGGAAATCAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATAGTGGCGATGATGtccggcgtatc                                               >  2:823805/1‑149 (MQ=255)
                                                                                                       ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                            >  4:876850/1‑149 (MQ=255)
                                                                                                       ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                            >  1:826881/1‑149 (MQ=255)
                                                                                                       ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                            >  1:70232/1‑149 (MQ=255)
                                                                                                           gTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccggcgta                                        <  1:895257/149‑1 (MQ=255)
                                                                                                             aTCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGa                                      >  3:344062/1‑149 (MQ=255)
                                                                                                              tCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAt                                     >  3:367297/1‑149 (MQ=255)
                                                                                                              tCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGAt                                     <  4:367480/149‑1 (MQ=255)
                                                                                                                tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGt                                                                           <  4:848140/109‑1 (MQ=255)
                                                                                                                tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGt                                                                           >  3:848140/1‑109 (MQ=255)
                                                                                                                tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTc                                   >  1:823638/1‑149 (MQ=255)
                                                                                                                tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTc                                   >  1:764484/1‑149 (MQ=255)
                                                                                                                tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTc                                   >  2:215349/1‑149 (MQ=255)
                                                                                                                  gAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGa                                 <  1:689402/149‑1 (MQ=255)
                                                                                                                    aaTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGaaa                               <  1:365729/149‑1 (MQ=255)
                                                                                                                     aTAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGGCGGTGGTAAAGCCGATTAAACTTTCAGATCACCGCCGAAGGGATAACCTTCACCCTGAATGGTGGCGATGATTTCCGGCGTATCCGGCCTAGATTCGAAAt                              >  3:698279/1‑149 (MQ=255)
                                                                                                                        aGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGtt                           >  3:372431/1‑149 (MQ=255)
                                                                                                                        aGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGtt                           <  1:334447/149‑1 (MQ=255)
                                                                                                                          gTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTa                         <  1:718095/149‑1 (MQ=255)
                                                                                                                          gTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTa                         <  2:668972/149‑1 (MQ=255)
                                                                                                                           ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                        >  3:4096/1‑149 (MQ=255)
                                                                                                                           ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                        >  3:438574/1‑149 (MQ=255)
                                                                                                                           ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                        <  4:704183/149‑1 (MQ=255)
                                                                                                                           ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                        >  3:771769/1‑149 (MQ=255)
                                                                                                                           ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                        <  2:516609/149‑1 (MQ=255)
                                                                                                                                 aaaaaaaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAc                  >  4:288164/4‑149 (MQ=255)
                                                                                                                                     aaaaCGGCGCTTTTTAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATAcggcg              >  2:515353/1‑149 (MQ=255)
                                                                                                                                       aaCGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGa            >  4:510936/1‑149 (MQ=255)
                                                                                                                                        aCGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGAt           <  4:372431/149‑1 (MQ=255)
                                                                                                                                            cgcTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTc       <  1:215349/149‑1 (MQ=255)
                                                                                                                                                 aaaaGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTc  >  2:48513/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTC  >  NC_000913/4639410‑4639703

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