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.5% 229.6 / ‑2.4 80H41R (CAC→CGC) fecDferric citrate ABC transporter membrane subunit FecD
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (46/31);  total (48/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.93e-01

CAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCG  >  NC_000913/4512046‑4512535
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
aacccccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCAAAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTc                                                                                                                                                                                                                                                 <  1:463990/248‑1 (MQ=255)
     ccGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGcc                                                                                                                                                                                                                                            >  1:286600/1‑251 (MQ=255)
      cGCCAGGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCAAGAGGGTTGCGCACAAGCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCt                                                                                                                                                                                                                                           <  1:262983/251‑1 (MQ=255)
       cccATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTg                                                                                                                                                                                                                                          <  1:181282/250‑1 (MQ=255)
             ccgcccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTc                                                                                                                                                                                                                                    >  1:113173/1‑251 (MQ=255)
              cgcccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCACGTCCGGCCTGCCAGTCg                                                                                                                                                                                                                                   >  1:331918/1‑251 (MQ=255)
                cccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGt                                                                                                                                                                                                                                 >  1:8788/1‑251 (MQ=255)
                 ccgcAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTc                                                                                                                                                                                                                                >  1:389935/1‑251 (MQ=255)
                  cgcAAAGGCCAGCAGCGGCAGCCCCACCTCGGGCGGTGCCGGCATAAGAAGTAGACCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCCCAACCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTca                                                                                                                                                                                                                               <  2:389935/251‑1 (MQ=255)
                      aaGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGCCGAGGGCTGCACCGACAAACAGTGCCAGCCGCAAGCGCGGCAGTCGGGACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTcagc                                                                                                                                                                                                                             >  1:427636/1‑249 (MQ=255)
                          ccAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCCGCAgcgcg                                                                                                                                                                                                                       >  1:461802/1‑251 (MQ=255)
                          ccAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCGGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAgcgcg                                                                                                                                                                                                                       <  1:115565/251‑1 (MQ=255)
                             gaagcgcaagcgccATCACGGGCAGTGACTGCATAAGAAGTAAGGCCCCCACAGAGGCCAGGCTGGCAGCTGGGTTAACACCGTGAAGATCGGGTGAGGCCAAAGGGTTGCGCACAATCCCCTGTATCACCGCTCCCGCCACGGCGAGGGGTGCACCGACAAACAGTGCCAGCAGCATGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAgcgcgcgc                                                                                                                                                                                                                    <  1:248248/239‑1 (MQ=255)
                                    cagcaCCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGAGCACAATCCCCTGTATCAGCACGCACGCCACGGCGAGGGCTGCACCGACAAACAGTGACAGCAGCAAGCGCGGCAGTCGGTACTCCCTCACTCAATAATAATGATCGCGGCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCcacggcc                                                                                                                                                                                                             >  1:159281/1‑250 (MQ=255)
                                         ccATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGTGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACgggg                                                                                                                                                                                                        >  1:404264/1‑251 (MQ=255)
                                          cATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACgggg                                                                                                                                                                                                        >  1:134053/1‑250 (MQ=255)
                                           atcccGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTGGCGCACAATCCCCTGTACCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGAt                                                                                                                                                                                                      <  1:186881/247‑1 (MQ=255)
                                               cGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCACTACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCCAGGCACGGGGCTCACt                                                                                                                                                                                                  >  2:363065/1‑251 (MQ=255)
                                                       gACGTCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCGGGGGGCCTGGTTAACCCCGAGAATACCCGGGGAGGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCCCGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATg                                                                                                                                                                                          <  2:120364/251‑1 (MQ=255)
                                                             aTAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAATCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTg                                                                                                                                                                                     >  1:276832/1‑250 (MQ=255)
                                                              aaagaagaagAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGAt                                                                                                                                                                                   <  1:297857/250‑1 (MQ=255)
                                                               aagaagGAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGGGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACGCCCATATGGAGTGATa                                                                                                                                                                                  <  1:17444/251‑1 (MQ=255)
                                                                agaagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAgcgc                                                                                                                                                                                                                        >  2:396028/1‑212 (MQ=255)
                                                                agaagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAgcgc                                                                                                                                                                                                                        <  1:396028/212‑1 (MQ=255)
                                                                agaagTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATa                                                                                                                                                                                  >  2:366407/1‑250 (MQ=255)
                                                                    gaggagCCCCCACGGAGGCCAGGGGGGCGGCATGGTAACCACCGAGAATACCGGGGGATGCAGGAGGGTTGCACACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAgc                                                                                                                                                                             <  2:124582/248‑1 (MQ=255)
                                                                        agCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGGCCGGCCTGCCAGTCGGTCAGCAGAGCGCGCCACGGCCCGGGGATCACTCCCCTATGGAGTGATAAAAGCGcac                                                                                                                                                                         >  2:276560/1‑251 (MQ=255)
                                                                          cccccACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACCTACTAATACTCGCGTCCGGCCGGCCAGTCGGTCAGCCGCGCTCGCACCGGCACGGGGCTCACTCCCATATGGAGTGATAACAGCGcacc                                                                                                                                                                        >  2:122845/1‑249 (MQ=255)
                                                                          cccccACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGcac                                                                                                                                                                         >  1:348938/1‑249 (MQ=255)
                                                                            ccaccgggggCCAGGCGGGCGGCAGGGTTAACCCCGAGAAAATCCGGTGAGGCCAGAGGGTTGCCCAAAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACAGAATAAGGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCTCACAGcc                                                                                                                                                                     <  2:195896/243‑1 (MQ=255)
                                                                                     gCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGc                                                                                                                                                            <  1:375142/251‑1 (MQ=255)
                                                                                     gCCAGGATGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGGATCAGCACGCACGCCACGGCGAGGGCTGCACCGACAAACAATGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACCTCATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCAAGGGGTGCACTCACATATGGAGAGATAACAGCGCACAGCCCGCAAGGGc                                                                                                                                                            >  2:206060/1‑251 (MQ=255)
                                                                                           cgggcggcATGGTTAACACCGAGAATATCCGGGGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGt                                                                                                                                                      <  1:420004/249‑1 (MQ=255)
                                                                                               gggcATGGTTAACACCGAGATTACCCGGTGATGCCAGAGGGTTGCGCCCAACCCCCTGAATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatg                                                                                                                                                  <  1:474417/250‑1 (MQ=255)
                                                                                               cggcATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGtgatg                                                                                                                                                  >  1:132107/1‑251 (MQ=255)
                                                                                                 gcATTGTTAACCCCGAGAATCCCCGGTGATGCCAGAGGGTTGCGCAAAACCCCCTGTAGCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAAGAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGaa                                                                                                                                                <  2:204492/251‑1 (MQ=255)
                                                                                                      gTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCAGCCACGGCGCGGGCTGCACCGACAAACAGTGCCAGCAGCACGCGCGGCCGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCAAGGGGATCACTCCCATATGGAGTGATAACAGCGCACACCCCTCCAGGGACAGGGtgctgcaaaaa                                                                                                                                            >  2:481277/1‑241 (MQ=255)
                                                                                                         aaacccGAAAATACCCGGTGATGCCAGAGGTTTGCGCACAACCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCa                                                                                                                                        <  1:151354/247‑1 (MQ=255)
                                                                                                                 gaATACCCGGTGAGGCCAGGGGGTTGCGCACAATCCCCTGTACCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAAtt                                                                                                                                <  1:353508/251‑1 (MQ=255)
                                                                                                                    tatCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCACGCGCGGCAGTCGGTACTCCATCACTACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTc                                                                                                                             >  1:47734/1‑251 (MQ=255)
                                                                                                                         ggTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCt                                                                                                                        >  1:402938/1‑251 (MQ=255)
                                                                                                                          gTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTc                                                                                                                       >  1:467131/1‑251 (MQ=255)
                                                                                                                                       gTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCACTTTTCATCCTCGCCTCCTcaca                                                                                                            >  1:244510/1‑249 (MQ=255)
                                                                                                                                                    cccTGTATCAGCACGCCCGGCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAAAACATAATAATGCTCGCGTCCGGCCTACCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGACAACAGCGCACAGCCCGCCAAGGCAAGGGTGATGAACATCACCAGCGCAATGTTCATCCTCGCCTCCCcacccaccagacacagaa                                                                                              >  2:281218/1‑235 (MQ=255)
                                                                                                                                                    cccTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGGACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATCTGGAGTGATAACAGCGCACCGCCCGCCAGGGCAAGGGGGATGAAAATAACCAGAGCCATTTTCATCCCCGCCTCCCCACAAGCCAGACAAAgaa                                                                                              >  2:385440/1‑248 (MQ=255)
                                                                                                                                                         tATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACCAAGAAAGGGcc                                                                                        >  2:161978/1‑250 (MQ=255)
                                                                                                                                                          aTCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTg                                                                                       <  1:87608/251‑1 (MQ=255)
                                                                                                                                                                       ccACGGCGAGGGCTTCACCGACAAACAGTGCCAGCGGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCCCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAAGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCGCAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCa                                                                          <  1:259082/251‑1 (MQ=255)
                                                                                                                                                                              gAGGGCTGCACCGACAAAAAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgc                                                                   >  1:119522/1‑251 (MQ=255)
                                                                                                                                                                               aGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTgcg                                                                  >  1:154039/1‑251 (MQ=255)
                                                                                                                                                                                  gCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGGACTCCATCAATCAATAATAATGCTCGCGTCCGGCCTGCCAGTCGGGCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGGGTGATAACATCGCACAGCCCGCCAGGGCAAGGGTGATGCAAATCACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCAGCCAATCCGCGCCAGCACTGCGCCt                                                               >  2:17703/1‑251 (MQ=255)
                                                                                                                                                                                  ccgacaccgcaagaCAGCGCCAGCGGCAAGCGCCTCTGTCGTTCCTCCACCAATCCATAATAAATCTGGCGGCCGGCCTGCCAGTCGGGCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATAGGGAGTGATAACAGCGCACAACCCGCCAGGGCAAGGGTGAGGAAAATAACCAGCGCAATGTACATCCTCGCCTCCTCACACGCCAGACAAAGCAAGGGCTGCCAATCATCGCCATCACTGCTCCt                                                               <  2:101537/238‑1 (MQ=255)
                                                                                                                                                                                      cACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCgg                                                           >  2:99965/1‑251 (MQ=255)
                                                                                                                                                                                       cccGACAAACAGGGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCggg                                                          <  2:170938/250‑1 (MQ=255)
                                                                                                                                                                                       aCCGACAAACAGGGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATAAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCCCGGGGCTCACTCCCATCTGGAGTGATAACAGCGCCCCGCCAGCCAGGGCAAGGGTGATGAAAATACCCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAACGCAAGGGCTGCCACTCAGCGCCAACACCGCCCCTTCggg                                                          >  2:475042/1‑251 (MQ=255)
                                                                                                                                                                                          gACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAg                                                       <  1:52532/251‑1 (MQ=255)
                                                                                                                                                                                                     ccAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCgggg                                             >  1:472307/1‑250 (MQ=255)
                                                                                                                                                                                                     acagcagCAAGCGCGGCAGTCGGTACTCCATCAATACATAATACTGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCCGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATATACgggg                                             >  2:291197/2‑250 (MQ=255)
                                                                                                                                                                                                      cagcagCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa                                           <  2:393459/251‑1 (MQ=255)
                                                                                                                                                                                                      cagcagCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGaa                                           <  1:377418/251‑1 (MQ=255)
                                                                                                                                                                                                       agcagcAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGa                                            >  1:442310/1‑249 (MQ=255)
                                                                                                                                                                                                         cccctacccgacgctcttccgatcaaTCAATACATAATAATGCGCGCGTCCGGCCGGCCAGGCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGc                                        <  2:43656/226‑1 (MQ=255)
                                                                                                                                                                                                            caAGCGCGGCAGTCGGTACTCGATCAATACATAATAATGCTCGCGTCCGGCATGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGTGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCAATTCGCCTGCGGGCAGATCTCCGGGGAAGGcc                                       >  2:316880/1‑249 (MQ=255)
                                                                                                                                                                                                              aGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCCCTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGCTGACAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAgtg                                   >  1:204858/1‑249 (MQ=255)
                                                                                                                                                                                                                cgcgGCAGGCGGTACTCCATCAATACATATTAATTCTCGCTTCCGGCCTGCCAGTCGGTCAGCGGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAgtgcg                                 <  1:441158/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggCAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGt                              >  1:300064/1‑251 (MQ=255)
                                                                                                                                                                                                                     cAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGc                            >  2:444647/1‑251 (MQ=255)
                                                                                                                                                                                                                     cAGTCGGTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGc                            <  1:148590/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTAc                      <  1:380962/251‑1 (MQ=255)
                                                                                                                                                                                                                           gTACTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTAc                      >  1:75444/1‑251 (MQ=255)
                                                                                                                                                                                                                             aCTCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACAt                    >  1:460399/1‑251 (MQ=255)
                                                                                                                                                                                                                               tCCATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCATCAGCGCGCGCCACGGCCCGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCATCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCCGCGCCAGCACTGCGCCCGCGGGCAGATCTCCGGGGAAGGCCAgtgcccggtcgaggacatcc                  >  2:135703/1‑232 (MQ=255)
                                                                                                                                                                                                                                ccATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGCCAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGAGAGCGGGCGAGTACATCTg                 >  2:104111/1‑251 (MQ=255)
                                                                                                                                                                                                                                ccATCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTg                 >  1:374916/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cATCAATACATAATAATGCTCGAGTCCGGCCTGACAGTCGGTCAGCAGCGCGCGCCCCGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCACCCGGGCACGGGTGATGAGACTAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAACCAAACCAAAGGATGCCAAGCAGAGCCAGCACTGCGCCTGCGGGACGCTCTCCGGGGAAGGCCCGTTCGCGGGCGCGTACATCTGc                >  2:87057/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aTCAATACATAATAATGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGc                                        >  1:43656/1‑226 (MQ=255)
                                                                                                                                                                                                                                   tCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTc                                                                                                                       <  1:98219/146‑1 (MQ=255)
                                                                                                                                                                                                                                   tCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTc                                                                                                                       >  2:98219/1‑146 (MQ=255)
                                                                                                                                                                                                                                     aacaCATAATAATGCCAGCGCCCGGCCTTCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGTAGTGATACCAGCGCACAGCCCGCCAGGGCTAGGGTGATGAAAATAACCAGCGCAATTTTCCTCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCCCTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACATCTGCcagc            <  2:353621/248‑1 (MQ=255)
                                                                                                                                                                                                                                               aatGCTCGCGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCAAGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTCCGGGCAGATCTCCGGGGAAGGCCAGTGCGCGTGCGAGTACCTCTGCCAGCAGCATCAGCg  >  2:428283/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
CAACCCCGCCATGCCGCCCGCAAAGGCCAGCAGCGGCAGCACCATCACGGGCAGTGACGGCATAAGAAGTAGAGCCCCCACAGAGGCCAGGCTGGCGGCATGGTTAACACCGAGAATATCCGGTGATGCCAGAGGGTTGCGCACAATCCCCTGTATCAGCACGCCCGCCACGGCGAGGGCTGCACCGACAAACAGTGCCAGCAGCAAGCGCGGCAGTCGGTACTCCATCAATACATAATAATGCTCGTGTCCGGCCTGCCAGTCGGTCAGCAGCGCGCGCCACGGCACGGGGATCACTCCCATATGGAGTGATAACAGCGCACAGCCCGCCAGGGCAAGGGTGATGAAAATAACCAGCGCAATTTTCATCCTCGCCTCCTCACAAGCCAGACAAAGCAAGGGCTGCCAATCAGCGCCAGCACTGCGCCTGCGGGCAGATCTCCGGGGAAGGCCAGCGCGCGTGCGAGTACATCTGCCAGCAGCATCAGCG  >  NC_000913/4512046‑4512535

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