Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 887,392 T→A 7.1% intergenic (‑53/‑31) ybjJ ← / → rcdA putative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913887,3920TA 7.1% 632.2 / 27.4 181intergenic (‑53/‑31)ybjJ/rcdAputative drug efflux MFS transporter, inner membrane protein/transcriptional regulator of csgD and ybiJI; autoregulator
Reads supporting (aligned to +/- strand):  ref base T (92/76);  new base A (6/7);  total (98/83)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.01e-01

Only 100 of 191 total aligned reads displayed.

AAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAG  >  NC_000913/887118‑887658
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                           
     ttttCCCTAACCACCCCTAACCTAGAATACCGCTCATCACACCGCTAGACAGACCAAAGAGAACCCCGGCCATGTCAGCGCTCGAGACAGAGAGGATATCGCGGAGAGCGGGCGTACAGGTGGCCCAGGACGCCATTAACAGGCAGAGAAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGAAAGGGTGTGTACATTTGTa                                                                                                                                                                                                                                                                        <  4:441651/274‑1 (MQ=255)
      tttCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGAAAGGTTGTGTACATTTGTac                                                                                                                                                                                                                                                                       >  2:114045/1‑275 (MQ=255)
                  cGCCCAGCAGAGATTGCAGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCaga                                                                                                                                                                                                                                                            <  4:178051/274‑1 (MQ=255)
                   gCCGAGCAGCGAATACCGCTCATCGAACCGATTGACACATTAACGAGAACACCGCCCACTTTAGCGACCGCGAGCGAGAGAATATCGCGGATAGCAGGCGAATGGGTCCCCCAGGCCGTCATTAACACGCCCGGCACAAAGAAGAACATTAACAGCTCACAGGTCTGGAGTTTCAATGCATTACGGCACGGATTTTCGTGCATAGATTACGTCAAAATAAGAAGTGGGCAGACAACATTAGACCGGCTGTTGACCTTTGTAAAAAATTGcgagaa                                                                                                                                                                                                                                                          >  4:314046/1‑270 (MQ=255)
                                   cGCCCATCGAACCGATCGACAGACCAAAGAGAACACAGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAAGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACAAAATTGCAGAGAGAGGAAATGACATGCg                                                                                                                                                                                                                                          <  2:517233/275‑1 (MQ=255)
                                          cGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                                                                                                                                        >  2:88217/1‑238 (MQ=255)
                                           gAACCGATCGACAGACCAAAGGGCACACCGCCCCTTTCAGCGATCGAGACAGAGAAAATATCGCGGATAACAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCACAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcg                                                                                                                                                                                                                                       <  4:114044/270‑1 (MQ=255)
                                              ccGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTAcaca                                                                                                                                                                                                                                                                    <  4:88216/238‑1 (MQ=255)
                                                              cgataaCACCGCCCATTTCAGCGATCGAGACAAAGAGAATATTGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACATGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgaga                                                                                                                                                                                                               <  2:327324/271‑1 (MQ=255)
                                                              agagAAAACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCgaga                                                                                                                                                                                                               <  2:411104/275‑1 (MQ=255)
                                                                       cGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGCGAAAAAATTAt                                                                                                                                                                                                      >  1:136533/1‑275 (MQ=255)
                                                                                                                 cAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTGCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGCTGTGTACATTTGTACACAATTGCAGCGAGAGGAAATGAAATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAACTATCCAGGCAACGCGGGAGGCGGGGAAACTTTCCAGAATACAtcc                                                                                                                                                            >  4:154511/1‑273 (MQ=255)
                                                                                                                            gTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGACGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGGGAGAGTAAATGACATGCGTCGCGCTACCTATCAGCAACGGCGAGAAAAAATTACCCAGTCCACGCTGGAGGCGGTGAAACTTACCGAAATACAATCTGTTTCGcct                                                                                                                                                  >  4:537637/1‑272 (MQ=255)
                                                                                                                                 ccAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                                                                                                                                          >  8:715116/1‑149 (MQ=255)
                                                                                                                                  cAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTg                                                                                                                                                                                                                                                                          >  6:275361/1‑148 (MQ=255)
                                                                                                                                        gCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTAcaca                                                                                                                                                                                                                                                                    <  6:69864/148‑1 (MQ=255)
                                                                                                                                        gCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaa                                                                                                                                                                                                                                                                   <  7:830505/149‑1 (MQ=255)
                                                                                                                                            ttAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTAcaca                                                                                                                                                                                                                                                                    <  6:893156/144‑1 (MQ=255)
                                                                                                                                                  ggCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                                                                                                                                                                                                         >  7:907007/1‑149 (MQ=255)
                                                                                                                                                  ggCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                                                                                                                                                                                                         <  7:247973/149‑1 (MQ=255)
                                                                                                                                                   gCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagag                                                                                                                                                                                                                                                         <  6:537286/148‑1 (MQ=255)
                                                                                                                                                    ccTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGACGAATTTACGGTCATAGATCACGTCAAAATAAAAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCCCGCCGCAGGCGGTGAACCTTTACGGAATACATTCTTTTAGCCCCACACAAATTGCTTCCCTTTCCgggg                                                                                                                        >  4:387550/1‑275 (MQ=255)
                                                                                                                                                    ccTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                                                                                                                                                                                                        <  5:884214/148‑1 (MQ=255)
                                                                                                                                                     cTGGCAAAAAGAAGAACATAAACAGAGCCCAGGTTCGGAGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGt                                                                                                                       <  2:298069/275‑1 (MQ=255)
                                                                                                                                                      tGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAgg                                                                                                                                                                                                                                                      >  6:589315/1‑148 (MQ=255)
                                                                                                                                                       ggCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                                                                                                                                        <  2:592210/129‑1 (MQ=255)
                                                                                                                                                            aaaGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTACACaa                                                                                                                                                                                                                                                                   >  5:585945/1‑129 (MQ=255)
                                                                                                                                                              agaagaACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACAt                                                                                                                                                                                                                                             <  7:817717/149‑1 (MQ=255)
                                                                                                                                                                aagaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGc                                                                                                                                                                                                                                           >  8:548820/1‑149 (MQ=255)
                                                                                                                                                                  gaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAttt                                                                                                                                                                                                                                                                           <  8:264794/115‑1 (MQ=255)
                                                                                                                                                                  gaaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                                                                                                                                        <  8:824150/118‑1 (MQ=255)
                                                                                                                                                                   aaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaaa                                                                                                                                                                             >  2:358345/1‑208 (MQ=255)
                                                                                                                                                                   aaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTAcaca                                                                                                                                                                                                                                                                    <  5:675472/121‑1 (MQ=255)
                                                                                                                                                                   aaCATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAAATGTAcaca                                                                                                                                                                                                                                                                    >  6:675472/1‑121 (MQ=255)
                                                                                                                                                                      aaaaaCAGCGACCAGGGTCGGGGTTTCGATACATTACGTGAAGAATTTACGGACATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCGAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATg                                                                                                       <  4:213250/272‑1 (MQ=255)
                                                                                                                                                                       tAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTcgcg                                                                                                                                                                                                                                     >  6:872837/1‑148 (MQ=255)
                                                                                                                                                                       tAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCgcnc                                                                                                                                                                                                                                    <  7:329186/149‑1 (MQ=255)
                                                                                                                                                                       tAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACttt                                                                                                                                                                         <  4:358343/208‑1 (MQ=255)
                                                                                                                                                                          aCAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                                                                                                                                                                                                        <  7:181307/126‑1 (MQ=255)
                                                                                                                                                                          aCAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCagagaga                                                                                                                                                                                                                                                        >  8:181307/1‑126 (MQ=255)
                                                                                                                                                                           cAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                                                                                                                                                                                                 <  6:535591/148‑1 (MQ=255)
                                                                                                                                                                           cAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                                                                                                                                                                                                 >  5:426099/1‑148 (MQ=255)
                                                                                                                                                                           cAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAAAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTaa                                                                                                                                                                                                                                 <  5:872461/148‑1 (MQ=255)
                                                                                                                                                                             gcgcCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACctacta                                                                                               >  1:143600/1‑275 (MQ=255)
                                                                                                                                                                                 ccAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                                                                                                                                                                                          >  7:699899/1‑149 (MQ=255)
                                                                                                                                                                                  cAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCg                                                                                                                                                                                                                          <  5:765361/148‑1 (MQ=255)
                                                                                                                                                                                          gCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCg                                                                                                                                                                                                                  <  5:720959/148‑1 (MQ=255)
                                                                                                                                                                                               tCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTGTTTTTCAGGAATtgatga                                                                             >  2:379561/1‑275 (MQ=255)
                                                                                                                                                                                                  aTGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGaaaaaa                                                                                                                                                                                                          >  6:631000/1‑148 (MQ=255)
                                                                                                                                                                                                        tACGTGAAGAATTTACGGCCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGCGAGGAAATGACATGCGTCGCGCTAGCGAACCGCAACGGCGAGAAAAAATTATCCAGGCCACGCGGGAGGCGGTGAAACTTTACGGAATACATGATGTTACGCACCGCAAAATTGATACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgttgattag                                                                    >  4:281081/1‑270 (MQ=255)
                                                                                                                                                                                                         aCGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                                                                                                                                                                                                  <  8:699899/149‑1 (MQ=255)
                                                                                                                                                                                                          cGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCtt                                                                                                                               >  2:770971/1‑214 (MQ=255)
                                                                                                                                                                                                              aagaaNTTACNGTCATAGATANTGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCg                                                                                                                           <  4:770969/214‑1 (MQ=255)
                                                                                                                                                                                                                    ttACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTgg                                                                                                                                                                                        >  5:502857/1‑148 (MQ=255)
                                                                                                                                                                                                                     tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCCAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATc                                                                                                                                                                                                     <  8:442883/134‑1 (MQ=255)
                                                                                                                                                                                                                     tACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATc                                                                                                                                                                                                     >  7:442883/1‑134 (MQ=255)
                                                                                                                                                                                                                           cATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACAtt                                                                                                                                                                                                                                                                            >  1:203577/1‑57 (MQ=255)
                                                                                                                                                                                                                             tAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGa                                                                                                                                                                               >  5:908217/1‑148 (MQ=255)
                                                                                                                                                                                                                               gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTa                                                                                                                                                                                                                                                                        <  3:203576/57‑1 (MQ=255)
                                                                                                                                                                                                                               gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCgg                                                                                                                                                                                  <  5:324229/143‑1 (MQ=255)
                                                                                                                                                                                                                               gATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGaaa                                                                                                                                                                             <  5:275361/148‑1 (MQ=255)
                                                                                                                                                                                                                                aTCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACt                                                                                                                                                                           >  7:596719/1‑149 (MQ=255)
                                                                                                                                                                                                                                 tCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGc                                                                                                                                                                                                <  1:46293/127‑1 (MQ=255)
                                                                                                                                                                                                                                   aCGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACttt                                                                                                                                                                         <  6:578242/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                                                                                                                                                        <  6:42812/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTa                                                                                                                                                                        <  5:467879/148‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTgagaac                                        >  2:541415/1‑273 (MQ=255)
                                                                                                                                                                                                                                     gTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACg                                                                                                                                                                      >  8:78717/1‑149 (MQ=255)
                                                                                                                                                                                                                                              agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACAt                                                                                                                                                              <  6:502857/148‑1 (MQ=255)
                                                                                                                                                                                                                                              agaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATg                                                                                                                                                             <  8:23989/149‑1 (MQ=255)
                                                                                                                                                                                                                                               gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                                                                                                                                                                                                  <  8:838494/111‑1 (MQ=255)
                                                                                                                                                                                                                                               gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAg                                                                                                                                                                                                  >  7:838494/1‑111 (MQ=255)
                                                                                                                                                                                                                                               gaagaGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGc                                                                                                                                                            <  8:907007/149‑1 (MQ=255)
                                                                                                                                                                                                                                                 agagGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCt                                                                                                                                                           >  6:387983/1‑148 (MQ=255)
                                                                                                                                                                                                                                                   agGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGtt                                                                                                                                                        <  8:150007/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       aaGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATcc                                                                                                                                                                                                    >  4:607970/1‑101 (MQ=255)
                                                                                                                                                                                                                                                        aGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGc                                                                                                                                                                                                <  6:86060/104‑1 (MQ=255)
                                                                                                                                                                                                                                                        aGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCt                                                                                                                                                                                          <  5:566790/110‑1 (MQ=255)
                                                                                                                                                                                                                                                         gACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCAc                                                                                                                                                  >  7:6903/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          acaacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCAc                                                                                                                                                  >  5:706337/1‑148 (MQ=255)
                                                                                                                                                                                                                                                            aacaTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGc                                                                                                                                                                                                                   <  7:229300/81‑1 (MQ=255)
                                                                                                                                                                                                                                                               aTTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATtgatg                                                                              >  2:206560/1‑210 (MQ=255)
                                                                                                                                                                                                                                                                   gCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAgtt                                                                          <  4:206559/210‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGt                                                                                                                                                         <  6:815052/131‑1 (MQ=255)
                                                                                                                                                                                                                                                                    cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTg                                                                                                                                       >  6:506244/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTg                                                                                                                                       >  6:171411/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                    cAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGc                                                                                                                                      <  8:803065/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aaGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGc                                                                                                                                      >  5:615746/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     aaGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCt                                                                                                                                     >  7:748318/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCg                                                                                                                           >  2:319152/1‑156 (MQ=255)
                                                                                                                                                                                                                                                                        gTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCCATATCAGGCAGTTTTTAGCGATGtt    >  2:436240/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                         ttGTGTACAAATGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTAccc                                                                                                                                 >  8:347418/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                         ttGTGTACAAATGTACACAATCGCAGAGAGAGGAAATGACATGCGTCGCGCTACCGATCCGCAACGGCGGGAAAAAAATATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATTCTGTTACGCACCGCAAAATTGCTAccc                                                                                                                                 >  8:293923/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGATAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTg                                                                                                                              >  7:512430/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                            tgtACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTg                                                                                                                              >  7:387425/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGAt                                                                                                                                                                                                                             <  7:267545/53‑1 (MQ=255)
                                                                                                                                                                                                                                                                              tACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGcc                                                                                                                            >  8:285972/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                cATTTGTACACAATTGCAGAGCGAGGAATTGACATGCGTCGCGCTAACGATCCGCAACGGCgaga                                                                                                                                                                                                               <  2:163811/65‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  tttGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGaa                                                                                    >  2:620085/1‑185 (MQ=255)
                                                                                                                                                                                                                                                                                  tttGTACACAACTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCggg                                                                                                                         >  5:83360/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                           
AAAGCGTTTCACTAACCACGCCGAGCAGAGAATACCGCTCATCGAACCGATCGACAGACCAAAGAGAACACCGCCCATTTCAGCGATCGAGACAGAGAGAATATCGCGGATAGCAGGCGTACGGGTTGCCCAGGACGCCATTAACAGGCCTGGCAAAAAGAAGAACATAAACAGCGCCCAGGTTCGGCGTTTCAATGCATTACGTGAAGAATTTACGGTCATAGATCACGTCAAAATAAGAAGAGGGAAGACAACATTAGCAAGGTTGTGTACATTTGTACACAATTGCAGAGAGAGGAAATGACATGCGTCGCGCTAACGATCCGCAACGGCGAGAAAAAATTATCCAGGCCACGCTGGAGGCGGTGAAACTTTACGGAATACATGCTGTTACGCACCGCAAAATTGCTACCCTTGCCGGGGTACCGTTGGGATCGATGACCTACTATTTTTCAGGAATTGATGAGTTGTTACTGGAGGCGTTCAGCAGTTTTACTGAGATCATGTCCCGGCAATATCAGGCATTTTTTAGCGATGTTAG  >  NC_000913/887118‑887658

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