Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,293 T→C H41R (CAC→CGC)  fecD ← ferric citrate ABC transporter membrane subunit FecD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,2930TC97.6% 239.2 / NA 83H41R (CAC→CGC) fecDferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (41/40);  minor base A (2/0);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 4.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

CCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGT  >  NC_000913/4512054‑4512536
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
ccATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCACCTGTATCAGCACGACCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCCATACATAATAATGCTCGCGTCCGGCCTGc                                                                                                                                                                                                                                            >  1:110758/1‑251 (MQ=255)
              aaGGCCAGCAGCGGCAGCACCAACACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCCCAGAGGCCAGGATGGCGGCATGGGTAACACCAAGAATAGCCAGTGATGCCAGCGGGATGCGCACACTCCCCAGGCTACGCACGCCCGCCACCGCGAGGGCAGCCCCACCAAACAGAGCCAGCATCAAGAGCAGCAGTCGGGACACCATCAACACATCATAATGCTCGCCGCCGGCCTGcccgacgagcaagcg                                                                                                                                                                                                                              >  2:399487/1‑238 (MQ=255)
               aGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCagcagc                                                                                                                                                                                                                             <  1:485270/251‑1 (MQ=255)
                  ccAGCAGCGGCAGCACCATCACGGGCAGTGACGGCCTAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGACTGCACCGACAAACAGAGCCAGCAACAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGATCGCGTCCGGCCTACCAGTCGGTcagcaacacg                                                                                                                                                                                                                          >  1:177850/1‑246 (MQ=255)
                        gcGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcac                                                                                                                                                                                                                   <  1:191373/251‑1 (MQ=255)
                            cagcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggcc                                                                                                                                                                                                               >  1:216126/1‑250 (MQ=255)
                             agcaCCATCCCGGGCGGTGACGGCATAGCAATGAGGGCCCCCACAGAGGCCAGGCGGGGGGCATGTTAAACCCCTAGAAAATCCGGTGAGGCAAGGGGGTGGGGCACAACCCCTGTTATCCGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggcac                                                                                                                                                                                                              <  2:655321/251‑1 (MQ=255)
                              gcaCCATCACGGGCAGGGACGGCATAAGAAGAAGAGCCCCCACAGAGGCCAGGCTGGCGTCATGTTTAACACCGAAAATATCCGGTGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCTCGCcacggcacg                                                                                                                                                                                                             <  1:483997/251‑1 (MQ=255)
                               ccccATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACCCCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggcacgg                                                                                                                                                                                                            <  1:626012/249‑1 (MQ=255)
                                     cACGGGCAGTGAAGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCAGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCAGCCAGACGGTCAGGAGCGCGCGCCACGGCCCGGGGaac                                                                                                                                                                                                       >  2:680173/1‑248 (MQ=255)
                                               accGGCAAAAGAATCAGAGCCCCCCCATGGGCCAGGCTGGCGGCGTGTTTAACACCGATAATACCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCTGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATg                                                                                                                                                                                            <  1:189600/249‑1 (MQ=255)
                                                  ggCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATATAATTCTCGCGTCCGGCCTGCCCGTCGGGTAACCGCGCGCGCCCACGCCCGGGGGTCCCTCCCATATGGAgg                                                                                                                                                                                        >  1:66962/1‑250 (MQ=255)
                                                         gaagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACCAACCGTGCCAGCAGCAAGCGCGGCAGTCGGTCCTCCCTCACTACATAATCATGCTCGCGTCCGGCCTGCCAGTCGGTCCGCAGCGCGCGCCACGGCACGGGGACCACTCCCATATGGAGTggtacc                                                                                                                                                                                  >  2:344200/1‑246 (MQ=255)
                                                          aagTAGCGCCCCACCAGAGGCCAGGTGGGCGGCATGGTTAACCCCTAGAACACCCGCGGATGCCAGGGGGTTGCCCCCAATCCCCGGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAACAGCGCGAGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACa                                                                                                                                                                                 <  2:141710/251‑1 (MQ=255)
                                                          aagTAGAGCCCCCACAGCGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACAGGGATCACTCCCATATGGAGTGATAAc                                                                                                                                                                                  >  1:299562/1‑250 (MQ=255)
                                                              agagCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAgcgc                                                                                                                                                                             <  1:279603/251‑1 (MQ=255)
                                                               gagCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAgcgc                                                                                                                                                                             >  1:510604/1‑250 (MQ=255)
                                                                 gCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGcac                                                                                                                                                                           >  1:390603/1‑250 (MQ=255)
                                                                   ccccACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGc                                                                                                                                                                        <  1:364822/251‑1 (MQ=255)
                                                                    cccACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGAAAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCCGCGCGCGCCACGGCACGGGGCTCACTCCCCTCTGGAGTGATAACAGCGcaccgat                                                                                                                                                                       >  2:489203/1‑247 (MQ=255)
                                                                    cccAAAGAGGCCAGGCTGGCGGCAGGGTTAACACCGAGAAGATCCGGGGATGCAAGAGGGTGGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGcac                                                                                                                                                                           <  1:489203/247‑1 (MQ=255)
                                                                      cacaGAGGCCTGGCGGGCGGCATGGTTAACCCCGAGAATATCCGGGGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCg                                                                                                                                                                     <  1:627318/251‑1 (MQ=255)
                                                                      cacaGAGGCCAGGCTGGGGGCATGGTTAACACCGAGAAAACCCGGTGATGCCAGGGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCg                                                                                                                                                                     <  1:481111/251‑1 (MQ=255)
                                                                       acaGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGc                                                                                                                                                                    >  1:177174/1‑251 (MQ=255)
                                                                        caGAGTCCAGGCTGGCGGCATGGCTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCGCAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCGGCACCGACAAACAGGGCCAGCAGCAAGCGCGGCAGACGGTACTCCCTCAATACATAATAATGCTCGCGGCGGGCCTGCAAGTCGGGCAGACGCGAGCGCCACGGCACGGGGATCCCTCCCATATGGCGGGATAAAAGCGCACCGccatcc                                                                                                                                                                   >  2:487965/1‑247 (MQ=255)
                                                                           gggCCAGGCTGGCGGCAGGGTAAACACCGAGAATACCCGGGGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAgg                                                                                                                                                                <  1:344200/250‑1 (MQ=255)
                                                                            ggCCAGGCTGGCGGCAGGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAggg                                                                                                                                                               <  1:216166/251‑1 (MQ=255)
                                                                               cAGGCTGGCGGCATGGTTAACCCCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCaa                                                                                                                                                           <  1:240902/251‑1 (MQ=255)
                                                                               cAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATACTGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCCCGGCACGGGGATCACTCCCATCTGGCGGGATACCAGCGCCCAGCCCACCAGGGCa                                                                                                                                                             >  2:594030/1‑250 (MQ=255)
                                                                                   ctagcagcATGGATAACACCGAGAATATCCGGTGCTGCCAGCGGGTTGCGCACCATCCCATGCATCAGCACGCCCGCCACGGCGAGGGCTGCAACGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCACGCATCCGGCCTCCCAGTCGGTAAGCAGCGCGCGCCACGGCCCGGGGGTCACTCCCCTATGGAGGGACACCAGCGCACCGACAGCCAGGGCAAggcg                                                                                                                                                       >  2:618863/7‑249 (MQ=255)
                                                                                       gggcATGGTTAACACCGAGAATATCCGGGGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatg                                                                                                                                                   <  1:307534/250‑1 (MQ=255)
                                                                                       cggcATGGTTACCCCCGAGAATATCCGGTGAGGCCAGAGGGTTGCCCACAACCCCCTGTATCAGCACGCCCGCCACGTCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACGCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatg                                                                                                                                                   <  1:519821/251‑1 (MQ=255)
                                                                                        ggcATGGTTAACACCGAGAATATCCGGTGATGCCAGCGGGTTGCGCACAATCACCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCACATGTGGAGTGATAACAGCGCACAGCCCGCCCGGGCACGGGtgatgg                                                                                                                                                  >  2:501706/1‑250 (MQ=255)
                                                                                          cATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTCTCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATAAATACATAATAATGCTCGCGGCCGGCCCGCCAGTCGGTCATCAGCGCGCACCACGGCACGGGGCTCACTCCCCTATGGCGTGATGACAGAGATCAGCCACCCAGGGCAAgggagaagaa                                                                                                                                                 >  2:462861/1‑243 (MQ=255)
                                                                                             ggtaaaccccGAGAATATCCGGTGATGCCAGAGGGGTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATa                                                                                                                                             <  1:461439/243‑1 (MQ=255)
                                                                                               ttAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCAAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCGCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGGGATAAAAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                                                           >  1:554159/1‑251 (MQ=255)
                                                                                                tAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAcc                                                                                                                                          >  1:168841/1‑251 (MQ=255)
                                                                                                aaaccccGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCGGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAc                                                                                                                                           <  1:490120/245‑1 (MQ=255)
                                                                                                 aaccccGAGAATACCGGGTGATGCAAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACGCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCa                                                                                                                                         <  1:86207/247‑1 (MQ=255)
                                                                                                 aaCACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATaa                                                                                                                                            >  1:461973/1‑248 (MQ=255)
                                                                                                 aaCACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAAcc                                                                                                                                          <  1:680173/250‑1 (MQ=255)
                                                                                                    cccGAGAATATCCGGGGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAgcg                                                                                                                                      <  1:408366/250‑1 (MQ=255)
                                                                                                            aacccggggaggCCAGGGGGTGGCGCCCAACCCCGTGTACCAGCGCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTc                                                                                                                              <  2:216126/240‑1 (MQ=255)
                                                                                                                 ggTGATACAAGAGGGTTGCGCACAAACCCCTGTATCCGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCACGCGCGGCAGTCGGTACTCCATCTATACATCATAATACTCGCGTCCGGCCTGCCAGTCGGGCAGCAGCGCCCGCCCCGGCACGGGGATCACCCCCATATGGGGTGATAACAGCGCACAGCCCGCCCTGGCACAGGGTATGAACATAAACACCAAAATTTTCATccc                                                                                                                         >  2:615036/1‑250 (MQ=255)
                                                                                                                     aTGCCAGAGGGTTGCGCACCATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGAAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCATGCCAATCGGTCAGCAGCGCGCGCCACGGCACGGGGACCACTCCCAAATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGATGATGAAAATAACCACCGCAATTTTCATCCTCGcc                                                                                                                     >  2:608884/1‑251 (MQ=255)
                                                                                                                        ccAGAGGGTTTCGCACAATCCACTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGCGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTAGGTCAGCAGCGCGCGCCACGGCACGGGGATCACACCCATATGGAGTGATAAAAGCGCACAGCCCGCCAAGGCAAGGGTGATGAAAATAACCCGCGAATTTTTGCTCCTCGcctcc                                                                                                                  >  2:154442/1‑251 (MQ=255)
                                                                                                                        ccAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATcatagcatcc                                                                                                                  >  1:354679/1‑242 (MQ=255)
                                                                                                                            agGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCCGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATCACAGCGCCCAACCCGCCAGGGCAAGGGTGATGAAAATAACAAGCGCAATTTACATCCTCCCATCCTcac                                                                                                              >  2:447173/1‑251 (MQ=255)
                                                                                                                             ggttggcgcACACACCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTcaca                                                                                                             <  1:34341/246‑1 (MQ=255)
                                                                                                                              ggTTGCCCACAACCCCCTGTACCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACaa                                                                                                            <  1:635451/251‑1 (MQ=255)
                                                                                                                                ttGCGCACAATCCACTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCACGGCAGTCGGTACTCCATCAAAACATAATAATGCTCGCGTCCGACCAGCCAGTCGGACAGCAGCGCTCGCCACGGGACGGTGATCACTCACATATGGAGGGATAACCACGCACAGCCTGCACGGGCGAGGGGGATGACAATAGACAGCGCAATTATAATCCCCGCCTCCACCCAAGc                                                                                                          >  2:473080/1‑251 (MQ=255)
                                                                                                                                   cgcACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACACACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGCGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGGGATGAAAATAACCAGCGCAATTTTACTCCTCGCCTCCCCACAAGcccg                                                                                                       >  2:668922/1‑249 (MQ=255)
                                                                                                                                      acaATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTAATCCTCGCCTCCTCACAAGCCAGAc                                                                                                     >  1:392629/1‑250 (MQ=255)
                                                                                                                                       caATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACa                                                                                                    >  1:298213/1‑250 (MQ=255)
                                                                                                                                        aaTCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACaaa                                                                                                  <  1:171269/251‑1 (MQ=255)
                                                                                                                                         aTCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACaa                                                                                                   >  1:335666/1‑249 (MQ=255)
                                                                                                                                           ccccTGTATCAGAACGCCCGCCACGGCGAGGGCTGCACCGAAAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCa                                                                                               <  1:50528/251‑1 (MQ=255)
                                                                                                                                            ccccatgtacAGCCCGCCCGCCACGGCGAGGCCTGCACCGACAAACAGTGCCATCAGAAAGTGCAGCAGGCGGTACTCCATCAATCCACAATAATGCCCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCCCGGCACGGGGATCACTCCCATATGGAGTGATAACAGCTCACAGCCCGCCATGGCAAGGGTGATGAAAAGAACCAGCGCAAATTTCAAACTCGCCTCCTCACAAGCCAGTCAAAgga                                                                                               <  2:82280/242‑3 (MQ=255)
                                                                                                                                              cggTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                                                            <  1:543649/249‑1 (MQ=255)
                                                                                                                                                tgtatacacccacccgccACGGCGAGAGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCACTACATAATACTGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCAAGGCCCGGGGATCACACCCATCTGGAGTGATACCAGCGCCCAGCCCGCCAGGGAACGGGTGAAGACACTACCCACCGAACATTTCATCCTCGCCTCATCCCACGCCAGACCCAGCAAggg                                                                                           >  2:631192/13‑251 (MQ=255)
                                                                                                                                                 gatccatcACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAgg                                                                                            <  1:674954/242‑1 (MQ=255)
                                                                                                                                                      aGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGcc                                                                                      >  1:541150/1‑250 (MQ=255)
                                                                                                                                                        cACGTCCGCCGCGGCGAGGGCTGCACCGACAAAAAGGGCCCGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCAGCCGGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCATTCCCAGATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAAt                                                                                   <  2:36877/251‑1 (MQ=255)
                                                                                                                                                          cgcccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCa                                                                                 <  1:658003/251‑1 (MQ=255)
                                                                                                                                                             ccgccACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAgcg                                                                              <  1:539663/251‑1 (MQ=255)
                                                                                                                                                              cgccACGGCAAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATCATGCTCGAGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCAGATGGAGTGATAACAGCGCACAGCCCGCCAGGGCCAGGGTGATGCAAATACACAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCATGGGCTGCCACTCAgcgc                                                                             >  2:189637/1‑251 (MQ=255)
                                                                                                                                                                 cacggatgggtcttacgcGACCAACAGTGCCAGCAGTCTGCGCGGCATTCGGTACTCCACCAATGCAGACTATGGCTTCCGCCCGTGCTGCCAGTCGGCCAGCAGCGCGTGCCACGCCACGGGGATCTCCCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTCTCATCCTCGCCTCCTCGCAAGCCAGACAAAGCAAGGGCTGCCAAGCAgcgcgag                                                                          <  2:531226/234‑4 (MQ=255)
                                                                                                                                                                      cGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTg                                                                     <  1:421164/251‑1 (MQ=255)
                                                                                                                                                                                aCCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCAGCCAATAAGCGCCAGCACTGCGCCTGCggg                                                           >  2:478636/1‑251 (MQ=255)
                                                                                                                                                                                       aaCAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGAtct                                                    <  1:105157/251‑1 (MQ=255)
                                                                                                                                                                                                agaagaaaGCGCGGAAGGCGGTACTCCATCAATACATAATAAGGCGCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGTCACTGGGATCACTCCCATATGGAGTGAGAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTAACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCCCTGAGCCTGCGGGCAGATCTACGGGGAAg                                           <  2:602750/245‑1 (MQ=255)
                                                                                                                                                                                                    gcaAGCGCGGCAGTCGGTACTCCATCAATACCTAATACTGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGAGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGAAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGcc                                        >  2:430545/1‑250 (MQ=255)
                                                                                                                                                                                                     caAGCGCGGATGCCGGTACTCAATCAATACATAATAATGCCCGAGTCCGGCCTGCCAGTCGGTCAGCAGCGCGAGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGACAGGGCAAGGGTGATGAAAATAACCAGCACAAGTTTCATCATCGCCTCCGCGCAAGCCAGACAAAGCACGGGCTGCCAATCAGAGCCAGCACTGCGACTGCGGACAGCTCTACGGGGAATGcccg                                      >  2:99958/1‑249 (MQ=255)
                                                                                                                                                                                                      aaGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTACGGGCAGATCTCCGGGGAAGGCCAgt                                     >  1:145764/1‑250 (MQ=255)
                                                                                                                                                                                                              cAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGc                             <  1:230304/251‑1 (MQ=255)
                                                                                                                                                                                                               aGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGAAAAAGCCAAGGATGCCACTCAGCGCCAGCACTGAGACTGCGGGCAGATCTACGGGGCAGGCCAggccgcggggg                            >  2:565520/1‑241 (MQ=255)
                                                                                                                                                                                                               aGTCGGAACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGAGCACTCCCATATGGAGTGATAACAGAGCACACCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCACTTTTCCTCCCCGCCACATCACAAACCAGACAAAGCAAGGGCGGCCAATCAGCGCCAGCCCTGCGCCTGCGGGCCGATCTCCGGGGAAGGCCAGTGCGCGGGCg                            >  2:21075/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGt                         >  1:267743/1‑251 (MQ=255)
                                                                                                                                                                                                                    gTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCTTCGCTTCCTCACAAGACAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTAc                       <  1:76339/251‑1 (MQ=255)
                                                                                                                                                                                                                             cAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGcc                >  2:133083/1‑249 (MQ=255)
                                                                                                                                                                                                                                   ataataATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCAATCCCATATGGAGTGATAACAGCGCACAGCACGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGACAGCAAGTGCGCGTGCGAGTACATCTGCCAGCAGCAt        >  1:51986/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aatGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCCAGCAGCATCAGc    <  1:501706/250‑1 (MQ=255)
                                                                                                                                                                                                                                        aatGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGGGCGCGTGCGGGTACATCTGCCAGCAGCATCAGCgg  >  1:669453/1‑250 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                     
CCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCGT  >  NC_000913/4512054‑4512536

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: