Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 832,254 T→A 5.6% intergenic (+18/+214) rhlE → / ← ybiA ATP‑dependent RNA helicase/DUF1768 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913832,2540TA 5.6% 517.6 / 15.6 144intergenic (+18/+214)rhlE/ybiAATP‑dependent RNA helicase/DUF1768 family protein
Reads supporting (aligned to +/- strand):  ref base T (57/79);  new base A (4/4);  total (61/83)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-01

Only 100 of 148 total aligned reads displayed.

CCGGACCCGTCAATCAAAGCCGAACCGATCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCT  >  NC_000913/832003‑832526
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                
ccGGCCCCGTCAATCAAAGCCGACCCGACCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAAc                                                                                                                                                                                                                                                            <  2:170125/274‑1 (MQ=255)
                     gacccGACCCAAAACGGTCGCCGGCAACGTGGCGGGGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATgcg                                                                                                                                                                                                                                      <  4:308799/272‑1 (MQ=255)
                     gAAACGATCCAGAAGGGTGGCAGGCAACGTGGCGGCGGGGGTGGTGGGCAAGGGGGTGGTCGCGGGCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATgcg                                                                                                                                                                                                                                      <  4:314513/275‑1 (MQ=255)
                           aTCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAACCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCAACGc                                                                                                                                                                                                                                >  2:308799/1‑275 (MQ=255)
                                   cGGTCGGCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAACACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                                        <  4:286427/275‑1 (MQ=255)
                                   cGGTCGCCCGGCACGGGGGGGGGGCGGGGNGGGGGAAGGGGGTGGCCGGGGTACCNAGCAACCACGCGGGGGGGAAGGGGGCGGAAAATCTGCAAGCGCGAAACCCGCAGAAAAACCGTCTCGCCGCCTCGGCGATNCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAgg                                                                                                                                                                                                                        <  3:60611/275‑1 (MQ=255)
                                                cGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGaa                                                                                                                                                                                                              <  2:274745/272‑1 (MQ=255)
                                                                    gaaaGGGGGTGGTCGGGGTCAACAGCAACCACGCCGGGGGGAAGGGGGCGAAAAATCTGCAAGCGCGAAACCGGCAGAAAACCCGGCTCGCCGCCTCGGCGAGGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTAGGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCa                                                                                                                                                                                       <  4:322398/273‑1 (MQ=255)
                                                                                  gcgGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGc                                                                                                                                                                          <  2:43025/274‑1 (MQ=255)
                                                                                                ccACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAAc                                                                                                                                                                                                                                                 >  4:263836/1‑189 (MQ=255)
                                                                                                    gCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                                                                                                                                                             <  2:263836/189‑1 (MQ=255)
                                                                                                    gCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTCAACTCTGCAATATATTTAATTTGCATGCTTTTGTAGGCCGGATAAGGCGGTTAcgcc                                                                                                                                                       >  2:41473/1‑275 (MQ=255)
                                                                                                          gggAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGcc                                                                                                                                                                                                                                                                               >  8:493370/1‑149 (MQ=255)
                                                                                                           ggAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGccc                                                                                                                                                                                                                                                                              <  8:1372719/149‑1 (MQ=255)
                                                                                                           ggAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGccc                                                                                                                                                                                                                                                                              <  7:493370/149‑1 (MQ=255)
                                                                                                               ggTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCATAGCCCGGCa                                                                                                                                                                                                                                                                          <  5:14282/149‑1 (MQ=255)
                                                                                                                  ggCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                                                                                                                                                       <  5:1160593/149‑1 (MQ=255)
                                                                                                                  ggCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                                                                                                                                                       <  6:1001497/149‑1 (MQ=255)
                                                                                                                   gcgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAg                                                                                                                                                                                                                                                                      <  6:1319159/149‑1 (MQ=255)
                                                                                                                    cgcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCGTTTATGCCGGGCTATGCCCGGCCCCAGGCTGATGAACAAAATCAAAACTGCCTGATTCCCTACGCTTATCAGGCCTACGTGAACTCTGACATATATTTAATTTCCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATTCTGCATGTca                                                                                                                                        >  4:137274/1‑274 (MQ=38)
                                                                                                                     gcAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCa                                                                                                                                                                                                                                                                       <  5:1285105/146‑1 (MQ=255)
                                                                                                                       aaaaTCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtg                                                                                                                                                                                                                                                                  <  5:197860/149‑1 (MQ=255)
                                                                                                                        aaaTCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                                                                                                                                                 <  6:896455/149‑1 (MQ=255)
                                                                                                                        aaaTCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                                                                                                                                                 >  8:1209467/1‑149 (MQ=255)
                                                                                                                         aaTCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtgat                                                                                                                                                                                                                                                                <  5:1255544/149‑1 (MQ=255)
                                                                                                                            cTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGaa                                                                                                                                                                                                                                                             <  6:1352295/149‑1 (MQ=255)
                                                                                                                            cTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGaa                                                                                                                                                                                                                                                             <  8:945203/149‑1 (MQ=255)
                                                                                                                                 aGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAAc                                                                                                                                                                                                                                                        >  7:993207/1‑149 (MQ=255)
                                                                                                                                      gAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCaaa                                                                                                                                                                                                                                                   >  6:1076585/1‑149 (MQ=255)
                                                                                                                                             gCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCt                                                                                                                                                                                                                                            <  5:306149/149‑1 (MQ=255)
                                                                                                                                              cAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                                                                                                                                                 >  6:608958/1‑127 (MQ=255)
                                                                                                                                              cAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCtga                                                                                                                                                                                                                                                                 <  5:608958/127‑1 (MQ=255)
                                                                                                                                               aGAAAAACCGTCTCGCTGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGa                                                                                                                                                                                                                                          <  7:1244198/149‑1 (MQ=255)
                                                                                                                                                  aaaaCCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGAt                                                                                                                                                                     <  2:318617/215‑1 (MQ=255)
                                                                                                                                                    aaCCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATgcgc                                                                                                                                                                                                                                     >  5:902440/1‑149 (MQ=255)
                                                                                                                                                     aCCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                                                                                                                                                              <  5:718081/139‑1 (MQ=255)
                                                                                                                                                       cGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTAc                                                                                                                                                                                                                                  >  5:199387/1‑149 (MQ=255)
                                                                                                                                                         tctcGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGc                                                                                                                                                                                                                                <  8:302625/149‑1 (MQ=255)
                                                                                                                                                             gccgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTAt                                                                                                                                                                                                                            >  8:657227/1‑149 (MQ=255)
                                                                                                                                                               cgccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          <  7:815517/149‑1 (MQ=255)
                                                                                                                                                                gccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAg                                                                                                                                                                                                                         >  5:947747/1‑149 (MQ=255)
                                                                                                                                                                 ccTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGTTGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTTTGGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTTTTAGCCATTTGAGAAAAAGGGAAAAAATTTAtt                                                                                          >  2:146394/1‑274 (MQ=255)
                                                                                                                                                                    cGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCt                                                                                                                                                                                                                     <  8:993207/149‑1 (MQ=255)
                                                                                                                                                                        gATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGc                                                                                                                                                                                                                                              >  8:85321/1‑120 (MQ=255)
                                                                                                                                                                          tGCCAAACCGGCAGGCGACCAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGTGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCTCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTAcc                                                                                 <  2:103146/275‑1 (MQ=255)
                                                                                                                                                                           gCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          >  4:367746/1‑137 (MQ=255)
                                                                                                                                                                           gCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGaa                                                                                                                                                                                                              <  7:1023650/149‑1 (MQ=255)
                                                                                                                                                                            ccAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAAc                                                                                                                                                                                                             <  7:203677/149‑1 (MQ=255)
                                                                                                                                                                             cAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTAt                                                                                                                                                                                                                                                                                  >  5:972523/1‑79 (MQ=255)
                                                                                                                                                                             cAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTAt                                                                                                                                                                                                                                                                                  <  6:972523/79‑1 (MQ=255)
                                                                                                                                                                              aaaCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACaaa                                                                                                                                                                                                                                                         <  6:1148427/103‑1 (MQ=255)
                                                                                                                                                                              aaaCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACaaa                                                                                                                                                                                                                                                         >  5:1148427/1‑103 (MQ=255)
                                                                                                                                                                              aaaCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAActc                                                                                                                                                                                                           <  7:194980/149‑1 (MQ=255)
                                                                                                                                                                              aaaCCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAActc                                                                                                                                                                                                           >  5:84739/1‑149 (MQ=255)
                                                                                                                                                                                       ggCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAatata                                                                                                                                                                                                  <  8:145068/149‑1 (MQ=255)
                                                                                                                                                                                        gCGAACAACAACGTCGCCGCCGTCCGCGTAAGCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAatatat                                                                                                                                                                                                 >  8:892878/1‑149 (MQ=255)
                                                                                                                                                                                         cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATAtt                                                                                                                                                                                                <  6:1314183/149‑1 (MQ=255)
                                                                                                                                                                                         cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATAtt                                                                                                                                                                                                <  6:1262065/149‑1 (MQ=255)
                                                                                                                                                                                         cGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATAtt                                                                                                                                                                                                <  6:1453273/149‑1 (MQ=255)
                                                                                                                                                                                          gAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTg                                                                                                                                                                                               >  6:226784/1‑149 (MQ=255)
                                                                                                                                                                                              aacaacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAAtt                                                                                                                                                                                           >  8:719749/1‑149 (MQ=255)
                                                                                                                                                                                                caacGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTTCGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTg                                                                                                                                                                                         <  7:892878/149‑1 (MQ=255)
                                                                                                                                                                                                    gTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATg                                                                                                                                                                                     <  7:9046/149‑1 (MQ=255)
                                                                                                                                                                                                         cgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCg                                                                                                                                                                                                                                                                             >  6:952116/1‑56 (MQ=255)
                                                                                                                                                                                                         cgccgTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCtttt                                                                                                                                                                                <  5:849535/149‑1 (MQ=255)
                                                                                                                                                                                                              tCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                                                                                                                                      >  8:574648/1‑106 (MQ=255)
                                                                                                                                                                                                              tCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAgg                                                                                                                                                                           >  7:1386335/1‑149 (MQ=255)
                                                                                                                                                                                                              tCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAgg                                                                                                                                                                           >  8:794376/1‑149 (MQ=255)
                                                                                                                                                                                                                cgcgTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCaa                                                                                                                                                                                                      >  2:121593/1‑120 (MQ=255)
                                                                                                                                                                                                                cgcgTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGAATGCTTTTGTAGGcc                                                                                                                                                                         <  6:1081086/149‑1 (MQ=255)
                                                                                                                                                                                                                   gTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGa                                                                                                                                                                      <  5:800608/149‑1 (MQ=255)
                                                                                                                                                                                                                      aaCCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATaa                                                                                                                                                                   >  7:1093480/1‑149 (MQ=255)
                                                                                                                                                                                                                        ccTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                                                                                                                                                             >  8:778038/1‑73 (MQ=255)
                                                                                                                                                                                                                        ccTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGcc                                                                                                                                                                                                                                             <  7:778038/73‑1 (MQ=255)
                                                                                                                                                                                                                         cTGCCGCTGCGCAGTAATCTTTTATGCCGGGTTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGc                                                                                                                                                                <  7:994757/149‑1 (MQ=255)
                                                                                                                                                                                                                         cTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGc                                                                                                                                                                <  7:343650/149‑1 (MQ=255)
                                                                                                                                                                                                                            ccGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGtt                                                                                                                                                             <  6:345447/149‑1 (MQ=255)
                                                                                                                                                                                                                             cGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGttt                                                                                                                                                            <  5:1076585/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cgcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgcc                                                                                                                                                       <  8:1451224/149‑1 (MQ=255)
                                                                                                                                                                                                                                   gcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                      >  7:493309/1‑149 (MQ=255)
                                                                                                                                                                                                                                   gcAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTAcgccg                                                                                                                                                      >  5:10483/1‑149 (MQ=255)
                                                                                                                                                                                                                                     aGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCa                                                                                                                                                    >  6:321193/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                  >  7:306576/1‑149 (MQ=255)
                                                                                                                                                                                                                                       tAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATc                                                                                                                                                  >  8:1288981/1‑149 (MQ=255)
                                                                                                                                                                                                                                         aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          <  7:1387802/75‑1 (MQ=255)
                                                                                                                                                                                                                                         aTCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCa                                                                                                                                                                                                                          >  8:1387802/1‑75 (MQ=255)
                                                                                                                                                                                                                                                aTGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTc                                                                                                                                         >  8:501497/1‑149 (MQ=255)
                                                                                                                                                                                                                                                 tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAAtt                                                                                                                                                                                           >  1:20518/1‑98 (MQ=255)
                                                                                                                                                                                                                                                 tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAttt                                                                                                                                          <  5:552677/147‑1 (MQ=255)
                                                                                                                                                                                                                                                 tGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAttt                                                                                                                                          >  6:552677/1‑147 (MQ=255)
                                                                                                                                                                                                                                                 tGCCGGGCATAGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTca                                                                                                                                        >  6:782176/1‑149 (MQ=255)
                                                                                                                                                                                                                                                  gCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTcac                                                                                                                                       >  7:908577/1‑149 (MQ=255)
                                                                                                                                                                                                                                                    cGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaa                                                                                                                                     <  7:1290203/149‑1 (MQ=255)
                                                                                                                                                                                                                                                     gggCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCa                                                                                                                                                                                       <  3:20518/98‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCacaaca                                                                                                                                   >  8:153741/1‑149 (MQ=255)
                                                                                                                                                                                                                                                       gCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCAcaacaa                                                                                                                                  <  5:321193/149‑1 (MQ=255)
                                                                                                                                                                                                                                                       gCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCAcaacaa                                                                                                                                  <  8:562975/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAg                                                                                                                                 <  7:657227/149‑1 (MQ=255)
                                                                                                                                                                                                                                                        cTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAg                                                                                                                                 >  7:630192/1‑149 (MQ=255)
                                                                                                                                                                                                                                                         tATGCTAGGCATCAGACTGAAGAACAAACGGAAAACAGCCGTGTGCGTTATGCGTATCAGGCACACGTGTACACTTCAATATAGTTAATTTGCATGCTTTTGTAGTCCGGATAAGGGGTTTACGCTGCATCCGGAATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTAACCACCCAGAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCt  <  4:146394/275‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                
CCGGACCCGTCAATCAAAGCCGAACCGATCCAGAACGGTCGCCAGCAACGTGGCGGCGGCGGTCGTGGGCAAGGTGGTGGTCGCGGTCAACAGCAACCACGCCGTGGGGAAGGTGGCGCAAAATCTGCAAGCGCGAAACCTGCAGAAAAACCGTCTCGCCGCCTCGGCGATGCCAAACCGGCAGGCGAACAACAACGTCGCCGCCGTCCGCGTAAACCTGCCGCTGCGCAGTAATCTTTTATGCCGGGCTATGCCCGGCATCAGGCTGATGAACAAACGCAAAACTGCCTGATGCGCTACGCTTATCAGGCCTACGTGAACTCTGCAATATATTGAATTTGCATGCTTTTGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCATTTCACAACAAGCACTTGTCAGCAATTTGAGAACACGGGAAAAAATTTATCTGTTTTACCGCCCATAAGCCACCGAGACGTTACTTCTCTATAGCCAATTGTTCGCGCAACTCCATTAAAAGGTAGCCCAGTCGATTCT  >  NC_000913/832003‑832526

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