Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,639,554 A→T 6.8% intergenic (+24/+36) 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,5540AT 6.8% 286.9 / 7.8 86intergenic (+24/+36)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 = 9.93e-01

GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTC  >  NC_000913/4639408‑4639698
                                                                                                                                                  |                                                                                                                                                 
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                                                                                                                                 <  4:664619/114‑1 (MQ=255)
                                                  tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑ttttttt                                                                                                                                 >  3:664619/1‑114 (MQ=255)
                                                  tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTGAAGCCGATTAATCTTCCAGATCAc                                                                                              >  1:304841/1‑149 (MQ=255)
                                                  tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCAc                                                                                              >  2:718095/1‑149 (MQ=255)
                                                  tGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCAc                                                                                              >  1:384400/1‑149 (MQ=255)
                                                     gTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                           >  2:789096/1‑149 (MQ=255)
                                                     gTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                           >  4:149643/1‑149 (MQ=255)
                                                     gTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGc                                                                                           >  3:156124/1‑149 (MQ=255)
                                                      ttttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCa                                                                                                    >  2:338516/1‑140 (MQ=255)
                                                      ttttCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCGTTTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCa                                                                                                    <  1:338516/140‑1 (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                                            >  1:516609/1‑149 (MQ=255)
                                                                                                    aTGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgta                                            >  3:880354/1‑149 (MQ=255)
                                                                                                      gAGTTACGTATCTGGAAATCAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATAGTGGCGATGATGtccggcgtatc                                          >  2:823805/1‑149 (MQ=255)
                                                                                                         ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                       >  1:70232/1‑149 (MQ=255)
                                                                                                         ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                       >  4:876850/1‑149 (MQ=255)
                                                                                                         ttACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTtccggcgtatccgg                                       >  1:826881/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                              >  2:215349/1‑149 (MQ=255)
                                                                                                                  tGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTc                              >  1:764484/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: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)
                                                                                                                             ttGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTAc                   >  3:4096/1‑149 (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)
                                                                                                                                                  |                                                                                                                                                 
GTGGTGGCGCTGGCCGGCATGATGTTTGTGACCCGTAATATCGACTGGTATGCGTTTTCACTGCCGAAAATGAAAGCCAGTAAAGAAGTTACAACGGACGATGAGTTACGTATCTGGAAATAAGGTTGAAAAATAAAAACGGCGCTAAAAAGCGCCG‑TTTTTTTTGACGGTGGTAAAGCCGATTAATCTTCCAGATCACCGCAGAAGCGATAACCTTCACCGTGAATGGTGGCGATGATTTCCGGCGTATCCGGCGTAGATTCGAAATGTTTACGAATACGGCGGATCGTC  >  NC_000913/4639408‑4639698

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: ‑