Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,291,512 G→A 5.7% M1M (GTG→ATG) †‡ nrfE → heme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,5120GA 5.7% 492.7 / 12.0 150M1M (GTG→ATG) †‡nrfEheme lyase (NrfEFG) for insertion of heme into c552, subunit NrfE
Reads supporting (aligned to +/- strand):  ref base G (74/67);  new base A (3/6);  total (77/73)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

Only 100 of 156 total aligned reads displayed.

TGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCG  >  NC_000913/4291250‑4291773
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                     
tggggggggtggtTCTGGTTTGGTGTGGCCGGGTTTGGGCGGTCTTTGCCAATGCGTGCCAAACCGTGGGTCAATGGCAGTTCCGGCATTCCTGCCGGGGTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAGGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCga                                                                                                                                                                                                                                                           <  4:446884/268‑2 (MQ=255)
       ggtggtTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGGTTCTGGACGTATTCCTTCCTGGAGTTGGTTTTCTGGCGTTGTTGTTAAGTATAGGGGCCAACGTGTTGACTCCGTTGACGGc                                                                                                                                                                                                                                                    >  2:516458/1‑275 (MQ=255)
        gtggtTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGCAGTCGGCTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcc                                                                                                                                                                                                                                                               >  2:446884/1‑263 (MQ=255)
                        gtctcagggctggggctgATTGTGCCAATGCTGCTCAAACCGGTGGTCATGCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCgg                                                                                                                                                                                                                                    <  1:144908/268‑1 (MQ=255)
                                              tgCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCCTGCCTGCCatttttc                                                                                                                                                                                                              >  4:103116/1‑269 (MQ=255)
                                                    tgctgcTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCCCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAccc                                                                                                                                                                                                         >  2:243281/1‑271 (MQ=255)
                                                         gctCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTg                                                                                                                                                                                                    >  2:149044/1‑274 (MQ=255)
                                                          ctCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGc                                                                                                                                                                                                   >  2:612098/1‑274 (MQ=255)
                                                                    gggTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTc                                                                                                                                                                                         >  2:291590/1‑274 (MQ=255)
                                                                                gTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCCCTGGCATCGGCATTCTGGCGCAGTTcg                                                                                                                                                                             >  2:331515/1‑274 (MQ=255)
                                                                                 ttCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCCGCGCTGTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGAGGTTGCTGGACGTATTCCTTCCTGAAGCCGGTTTTCCGGCGTTGTTGTTAAGTCTTGGGGTCCACGTGGTGACCCCGTTGATGGCCTTCGCGGGAGTGCGGTTGCCCTTGCCTGCCATGACGCCACTCCCTTGCAGCGGTATTCTGGCGCATCTcgc                                                                                                                                                                            >  4:308960/1‑274 (MQ=255)
                                                                                  tCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgcg                                                                                                                                                                           <  4:149044/274‑1 (MQ=255)
                                                                                           ccTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCGTGCATGATGCGACGCACTTGCATCGGCATTCTGGCGCAGTTCTCGCTCCtgctgc                                                                                                                                                                 >  1:90369/1‑275 (MQ=255)
                                                                                                                   gCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACg                                                                                                                                                                                                                                                                       >  8:734098/1‑149 (MQ=255)
                                                                                                                   gCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACATATTCCTTGCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACg                                                                                                                                                                                                                                                                       >  7:932736/1‑149 (MQ=255)
                                                                                                                       gTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                                                                                                                                   <  7:143623/149‑1 (MQ=255)
                                                                                                                         cTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGa                                                                                                                                                                                                                                                                 <  5:536919/149‑1 (MQ=255)
                                                                                                                           ggTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGAcc                                                                                                                                                                                                                                                               >  5:663942/1‑149 (MQ=255)
                                                                                                                            gTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGAccc                                                                                                                                                                                                                                                              <  8:114188/149‑1 (MQ=255)
                                                                                                                                gCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGtt                                                                                                                                                                                                                                                          <  5:720818/149‑1 (MQ=255)
                                                                                                                                  gtgtTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGa                                                                                                                                                                                                                                                        <  6:102546/149‑1 (MQ=255)
                                                                                                                                   tgtTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGAc                                                                                                                                                                                                                                                       <  6:327963/149‑1 (MQ=255)
                                                                                                                                    gtTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                                                                                                                                           >  5:356043/1‑143 (MQ=255)
                                                                                                                                    gtTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACACGTTGACCCCga                                                                                                                                                                                                                                                           <  6:356043/144‑2 (MQ=255)
                                                                                                                                         tgctgcGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGGTGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCtt                                                                                                                                                                                                                                                 <  6:670177/149‑1 (MQ=255)
                                                                                                                                           ctgcGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCg                                                                                                                  <  2:308960/274‑1 (MQ=255)
                                                                                                                                            tgcGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACgtg                                                                                                                                                                                                                                                                     <  8:432951/126‑1 (MQ=255)
                                                                                                                                              cgcTTTTTCATTCACTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAACGAGGTGGATTCTGGACGTCTTCCTACCCGAAGTCGGTATTCTAGCGTTGTTGTTAAATCTGGCGGTAAACGTGTAGACCCCGTTGAAGGCTTTCgggg                                                                                                                                                                                                                                            <  6:475543/149‑1 (MQ=255)
                                                                                                                                               gcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCggg                                                                                                                                                                                                                                           >  7:785913/1‑149 (MQ=255)
                                                                                                                                               gcTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                                                                                                                                              >  4:199872/1‑162 (MQ=255)
                                                                                                                                                 tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                                                                                                                                                                                         >  6:1008615/1‑149 (MQ=255)
                                                                                                                                                 tttttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAg                                                                                                                                                                                                                                         >  5:188355/1‑149 (MQ=255)
                                                                                                                                                   tttCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc                                                                                                                                                                                                                          <  2:199872/162‑1 (MQ=255)
                                                                                                                                                        ttCTCTACGCCCGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGtt                                                                                                                                                                                                                                  >  5:576942/1‑149 (MQ=255)
                                                                                                                                                           tctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                                                                                                                                               >  5:966601/1‑149 (MQ=255)
                                                                                                                                                           tctACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcg                                                                                                                                                                                                                               >  7:357359/1‑149 (MQ=255)
                                                                                                                                                             tACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                                                                                                                                                                                             >  5:675600/1‑149 (MQ=255)
                                                                                                                                                             tACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCt                                                                                                                                                                                                                             >  5:243628/1‑149 (MQ=255)
                                                                                                                                                             tACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCCTTCTGGCGCAGTTCGCGCCCCTGCTGCTCGCCTTTGGCGTACTCACGGATTGTTTTCTCATCAGCGATTTCTCGGTCCTTTTTgactcacaa                                                                                                >  4:419285/1‑266 (MQ=255)
                                                                                                                                                             tACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCaa                                                                                                >  4:489843/1‑274 (MQ=255)
                                                                                                                                                               ggccgggaggTTGGCGGTGGCGTCATCCATAAAAGGGGTGGTTTCTGGACGTATTCCTTCCTGAGGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACa                                                                                              <  4:243281/265‑1 (MQ=255)
                                                                                                                                                                 ccGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTcg                                                                                                     <  2:419285/265‑1 (MQ=255)
                                                                                                                                                                  cGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCt                                                                                                                                                                                                                        <  7:901485/149‑1 (MQ=255)
                                                                                                                                                                    gACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCta                                                                                        <  2:489843/275‑1 (MQ=255)
                                                                                                                                                                     aCAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                                                                                                                                                                                     <  8:64835/149‑1 (MQ=255)
                                                                                                                                                                     aCAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc                                                                                                                                                                                                                     >  5:756028/1‑149 (MQ=255)
                                                                                                                                                                         ttGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatga                                                                                                                                                                                                                 >  5:694733/1‑149 (MQ=255)
                                                                                                                                                                         ttGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGc                                                                                                                                                                                                   <  4:166344/163‑1 (MQ=255)
                                                                                                                                                                          tGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCct                                                                                  >  2:404562/1‑275 (MQ=255)
                                                                                                                                                                           gACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCatgatg                                                                                                                                                                                                               <  5:1096771/149‑1 (MQ=255)
                                                                                                                                                                               gTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGAc                                                                                                                                                                                                           <  6:586526/149‑1 (MQ=255)
                                                                                                                                                                                    gTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACt                                                                                                                                                                                                      >  7:197348/1‑149 (MQ=255)
                                                                                                                                                                                        gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAATTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGaa                                                                     <  2:22299/274‑1 (MQ=255)
                                                                                                                                                                                        gCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGCGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCACCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGGTGTCGTGGGAAc                                                                    >  2:151963/1‑275 (MQ=255)
                                                                                                                                                                                            gtaaaGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGCCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCgg                                                                                                                                                                                              >  5:836243/3‑149 (MQ=255)
                                                                                                                                                                                                 gagGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTc                                                                                                                                                                                         <  7:731956/149‑1 (MQ=255)
                                                                                                                                                                                                     gggtTTGTGGATGTATTCCTGCTTGAAGTCGGGTGTTGGGTGTGGTTGTTTAGTCGCGGGGTGAACGTGTCGAGCCCGGGGACGGCCTTCGCGGGAGGGCGGTTGCGCTGGCCTGCCAGTATGCGACGCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATGGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAg                                                        <  4:516458/273‑1 (MQ=255)
                                                                                                                                                                                                       gtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGgcgc                                                                                                                                                                                   >  7:25889/1‑149 (MQ=255)
                                                                                                                                                                                                       gtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTcgcg                                                                                                   >  2:317610/1‑229 (MQ=255)
                                                                                                                                                                                                       gtTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTcgcg                                                                                                   <  4:317610/229‑1 (MQ=255)
                                                                                                                                                                                                        tttCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc                                                                                                                                                                                                                              <  4:659635/105‑1 (MQ=255)
                                                                                                                                                                                                           cTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGtt                                                                                                                                                                               >  6:724497/1‑149 (MQ=255)
                                                                                                                                                                                                           cTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCg                                                                                                                  <  2:667393/210‑1 (MQ=255)
                                                                                                                                                                                                            tGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTc                                                                                                                                                                              <  6:629176/149‑1 (MQ=255)
                                                                                                                                                                                                             ggACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcg                                                                                                                                                                             >  8:324244/1‑149 (MQ=255)
                                                                                                                                                                                                              gACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTcgc                                                                                                                                                                            <  7:986134/149‑1 (MQ=255)
                                                                                                                                                                                                                cGTTTTCCTTCCTTAGTTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggc                                             <  4:386846/274‑1 (MQ=255)
                                                                                                                                                                                                                 gTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGggcg                                            >  4:47934/1‑274 (MQ=255)
                                                                                                                                                                                                                   aTTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTCCGCGGGAGTGCGGTTGCCCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTcc                                                                                                                                                                       <  7:669967/149‑1 (MQ=255)
                                                                                                                                                                                                                    ttccttccTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCct                                                                                                                                                                      >  5:450001/1‑149 (MQ=255)
                                                                                                                                                                                                                       cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgc                                                                                                                                                                    >  6:155745/1‑148 (MQ=255)
                                                                                                                                                                                                                       cttcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCctgc                                                                                                                                                                    <  5:155745/148‑1 (MQ=255)
                                                                                                                                                                                                                         tcctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTgg                                                                                                                                                       >  4:73100/1‑159 (MQ=255)
                                                                                                                                                                                                                           ctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                                                                                                                                     <  7:522739/111‑1 (MQ=255)
                                                                                                                                                                                                                           ctGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                                                                                                                                     >  8:522739/1‑111 (MQ=255)
                                                                                                                                                                                                                             gAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTa                                                                                                                                                   <  2:73100/159‑1 (MQ=255)
                                                                                                                                                                                                                                gTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCtt                                                                                                                                                          >  5:13205/1‑149 (MQ=255)
                                                                                                                                                                                                                                  cGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTg                                                                                                                                                        <  6:675600/149‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTg                                                                                                                                                        <  7:486618/149‑1 (MQ=255)
                                                                                                                                                                                                                                   ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGtt                                                                                                                                                                                                                                                                   <  6:696910/41‑1 (MQ=255)
                                                                                                                                                                                                                                   ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTgg                                                                                                                                                       >  6:224126/1‑149 (MQ=255)
                                                                                                                                                                                                                                   ggTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTgg                                                                                                                                                       <  8:217427/149‑1 (MQ=255)
                                                                                                                                                                                                                                       ttCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATg                                                                                                        <  2:485091/192‑1 (MQ=255)
                                                                                                                                                                                                                                        tCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTAc                                                                                                                                                  <  6:796957/149‑1 (MQ=255)
                                                                                                                                                                                                                                            gCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGAc                                                                                                                                              <  5:232943/149‑1 (MQ=255)
                                                                                                                                                                                                                                            gCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGAc                                                                                                                                              <  6:174740/149‑1 (MQ=255)
                                                                                                                                                                                                                                              gttgttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGt                                                                                                                                            <  5:1008615/149‑1 (MQ=255)
                                                                                                                                                                                                                                               ttgttgttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTa                                                                                                                                           <  7:324244/149‑1 (MQ=255)
                                                                                                                                                                                                                                                    gttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGtt                                                                                                                                      >  7:1005548/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCt                                                                                                                                                           >  8:80777/1‑127 (MQ=255)
                                                                                                                                                                                                                                                       aaGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTtc                                                                                                                                   <  6:450001/149‑1 (MQ=255)
                                                                                                                                                                                                                                                           aacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTgg                                                                                                                                                       <  8:505304/123‑1 (MQ=255)
                                                                                                                                                                                                                                                           aacGGGGTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTgg                                                                                                                                                       >  7:505304/3‑125 (MQ=255)
                                                                                                                                                                                                                                                            tcGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCtcatc                                                                                                                              <  8:194002/149‑1 (MQ=255)
                                                                                                                                                                                                                                                               gggTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGc                                                                                                                           >  8:1013074/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                                                                                                                                     <  8:449712/75‑1 (MQ=255)
                                                                                                                                                                                                                                                               gggTCAACACGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACtt                                                                                                                                                                                                     >  7:449712/1‑75 (MQ=255)
                                                                                                                                                                                                                                                                   cAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGAtt                                                                                                                       <  7:861769/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                     aCGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtc                                                                                                                     >  5:4142/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                       gtgtTGACCCCGTTGACGGCCCTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTtctc                                                                                                                   <  7:431677/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                     
TGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCATGATGCGACTCACTTGCATCGGCATTCTGGCGCAGTTCGCGCTCCTGCTGCTCGCCTTTGGCGTACTGACGTATTGTTTTCTCATCAGCGATTTCTCGGTCATTTATGTCGCGCAACATAGCTATAGCCTGCTGTCGTGGGAACTCAAGCTGGCAGCGGTGTGGGGCGGCCATGAAGGTTCGCTGCTGCTTTGGGTGCTGCTGCTTTCCG  >  NC_000913/4291250‑4291773

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