Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,639,556 A→T 6.7% intergenic (+26/+34) 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,5560AT 6.7% 288.1 / 7.8 86intergenic (+26/+34)creD/arcAinner membrane protein/response regulator in two‑component regulatory system with ArcB or CpxA
Reads supporting (aligned to +/- strand):  ref base A (47/33);  new base T (3/3);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 6.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTC  >  NC_000913/4639408‑4639703
                                                                                                                                                    |                                                                                                                                                    
gtggtggCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTaaa                                                                                                                                                      >  1:158904/1‑149 (MQ=255)
  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)
                                                                                                                                                    |                                                                                                                                                    
GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTCACGTC  >  NC_000913/4639408‑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: ‑