Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 13,320 A→G V386V (GTA→GTG dnaK → chaperone protein DnaK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091313,3200AG98.7% 234.8 / ‑5.4 75V386V (GTA→GTGdnaKchaperone protein DnaK
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (33/41);  minor base T (0/1);  total (33/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

AACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGC  >  NC_000913/13073‑13566
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
aaCTGGAAAGCCTGGTTGAAGAGCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTgaaa                                                                                                                                                                                                                                                     <  2:479118/251‑1 (MQ=255)
    ggAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAATCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGAc                                                                                                                                                                                                                                                  <  1:170743/250‑1 (MQ=255)
           cTGGTTGCAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGCCGTGAAAGACGTActgct                                                                                                                                                                                                                                          >  2:612562/1‑251 (MQ=255)
                 gAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCGGGTCAGGGTGGTGTGCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGa                                                                                                                                                                                                                                    <  2:65167/251‑1 (MQ=255)
                  aaGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGAc                                                                                                                                                                                                                                   <  1:206573/251‑1 (MQ=255)
                         ggTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTAccc                                                                                                                                                                                                                            <  1:216500/251‑1 (MQ=255)
                         ggTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTAc                                                                                                                                                                                                                              <  1:612562/249‑1 (MQ=255)
                             aaCCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCt                                                                                                                                                                                                                        >  1:209473/1‑251 (MQ=255)
                             aaCCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCGGGCCTGGTCGTATCTGATATCGACGACGGTAGCCGCGTTGGTGGTCAGACTCGTTTGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCt                                                                                                                                                                                                                        <  2:431767/251‑1 (MQ=255)
                                 gTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGAAGCTGGCCTGTCCGTATCTGATATCGACGCCGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGGAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATAGGTGCTGCTGTTCAGGGGGGTGTTCTGACTGGTGACGTGCAAGACGTACTGCTGCTGGACGTTCCCCCGCTGtct                                                                                                                                                                                                                    >  2:248469/1‑251 (MQ=255)
                                     cagtGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGCCTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGCCTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTgg                                                                                                                                                                                                                <  2:57527/248‑1 (MQ=255)
                                         gAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTAt                                                                                                                                                                                                            <  1:412961/251‑1 (MQ=255)
                                         gAGCCGATGAAAGATGCACTGAAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTAGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCCGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCAGCTGTCTCTGGGTAt                                                                                                                                                                                                            >  2:286363/1‑251 (MQ=255)
                                             agCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGa                                                                                                                                                                                                         >  2:351608/2‑250 (MQ=255)
                                              gCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGGCTCTGGGTATCGaaa                                                                                                                                                                                                       <  2:250862/251‑1 (MQ=255)
                                                 gAAAGTTGCACTGCAGGACGCTGGCTTGTCCGTATCTGATATCGACGTGGTTATCCTCGTTGGTGGTCAGACGCGTATGCCAATGGTTCAGAAGAAAGTTGCTGTGTTCTTTGGGAAAGAGCCGCGTAAAGACGGTAACCCGGACGAAGCGGTAGCAAGCGGTGCGGCTGGTTAGGGTGGTGTTGTGACTGGTTACTGGAAAGACGTACTGCTGCTGGACGGTACCCCGCTGTCTCTGGGTATCGAAACCa                                                                                                                                                                                                    <  2:299716/251‑1 (MQ=255)
                                                 gAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGGCAGACTCGTATGCCAATGGTTCAGAAGAAAGGTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCa                                                                                                                                                                                                    <  2:35673/251‑1 (MQ=255)
                                                     gTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATggg                                                                                                                                                                                                >  1:68095/1‑251 (MQ=255)
                                                      ttGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGGGACGTGAAAGACGTCATGCTGCTGGAAGTTACCCCGCTGCCTCTGGGTTTCGAAACCAtcgcg                                                                                                                                                                                               >  2:193032/1‑247 (MQ=255)
                                                       tGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGAAGCAATCGGCGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGAATCGAAACCATGGGCg                                                                                                                                                                                              >  1:67117/1‑251 (MQ=255)
                                                        ccgatgcattaCGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTAGCGAAACCATGGGCgg                                                                                                                                                                                             <  2:485689/241‑1 (MQ=255)
                                                         cacggCAGGACGCTTGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACGCGTATGCCAATGGTTCAGAAGAAAGTTGCTTAGTTCTTTGGTTATGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGTTGACGTTAAAGACGTACTGCTGCTGGTCGTGACCCCGCTGTCTCTTGGTAGCGAAACCATGGGCGgt                                                                                                                                                                                            <  1:429450/247‑1 (MQ=255)
                                                          aCTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTGACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGgtg                                                                                                                                                                                           <  2:257232/251‑1 (MQ=255)
                                                                ggACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGtgatga                                                                                                                                                                                     <  1:396952/251‑1 (MQ=255)
                                                                ggACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGtgatg                                                                                                                                                                                      >  1:175415/1‑250 (MQ=255)
                                                                 gACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGCCTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGCAAGACGTACTGCTGCTGGACGTTACCCAGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgac                                                                                                                                                                                    >  2:267763/1‑251 (MQ=255)
                                                                       gtccTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGacgacg                                                                                                                                                                                <  2:484449/247‑1 (MQ=255)
                                                                       ggCCTGTCCGTATCTGATATCGACGACGTTACCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAATAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGGGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTTCTGCTGCTGGACGTTACCCCGCTGTCGCTGGGTATCGAAACCATGTGCGGTGTGATGACGACGCt                                                                                                                                                                              <  2:554252/251‑1 (MQ=255)
                                                                        gCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTg                                                                                                                                                                             <  1:379071/251‑1 (MQ=255)
                                                                         ccTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGa                                                                                                                                                                            <  1:62082/251‑1 (MQ=255)
                                                                           ggTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTTGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATc                                                                                                                                                                          <  1:248469/250‑1 (MQ=255)
                                                                            gTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCAGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGAt                                                                                                                                                                           <  1:267763/249‑1 (MQ=255)
                                                                               cGTATCTGATATCGACGCCGTTATCCCCGTTGGTGGACAGACTCGTCTGCCAATGGTTCAGCAGCAAGATGCTGAGTTCTTTGGCCAAGTGCCGCGTACAGACGATAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAACCCATGGGCGGTGTGATGCCCACGCTGATcgcgc                                                                                                                                                                      >  2:284353/1‑250 (MQ=255)
                                                                                gTATATGATATCGACGAGGTTAGCCACGTTGTGGGTCAGACTCGTATGACAATGGTTCAGAAGAAAGTTGCTGAGTTCTATGGAAAAGAGCCGCGTAAAGACGGGAACGCGGACGACGCTGGAGCAACCGGTGCTGCTGGTCAGGTTCGGGTTGTGGCTGGGGACGTGGAATACGTACGGCTGCTGGACGCTACCTCGGTGCCTCTGGGTATCGAAGCCAGGGGCCGGGTGAGGACGACGGTGATCGCGaa                                                                                                                                                                     <  2:275980/251‑1 (MQ=255)
                                                                                      gATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaaaa                                                                                                                                                                  <  2:517080/248‑1 (MQ=255)
                                                                                       atatCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAAcacc                                                                                                                                                              >  1:67067/1‑251 (MQ=255)
                                                                                            gacgacGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGGTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAGCCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                                                                                                                                         <  1:419170/251‑1 (MQ=255)
                                                                                            gacgacGTTATCCTCGTTGGTGGTCAGACGCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                                                                                                                                         <  1:5916/251‑1 (MQ=255)
                                                                                             acgacgTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                                                                                                                                        >  1:123096/1‑251 (MQ=255)
                                                                                                              ggtggtCAGCCTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGcc                                                                                                                                        >  1:495190/1‑250 (MQ=255)
                                                                                                                  gtCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTGCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtg                                                                                                                                   <  1:543454/251‑1 (MQ=255)
                                                                                                                          cGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAc                                                                                                                           <  1:351608/251‑1 (MQ=255)
                                                                                                                            tATGCCAAGGGGTCAGAAGAAAGTTGGTGAGTTCGTGGGTAAATAGCTGCTTAAAGACGTTAACCCGGACGAAGCTGTATCAATCGGGTCTGCTGTTCAGGGTGGTGTTCTGACTGGTGTCGTGAAAGACGTACTGCTGCTGGACGTTTCCCCGCTGTCTCTGTGGATCGAAACCAGGGGCGGTTTGATGACGACGCTGATCGCGAAATACGCCACTATCCCGACCAAGCACAGCCAGGTGGTCTCTAc                                                                                                                           <  2:375631/249‑1 (MQ=255)
                                                                                                                             aTGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGc                                                                                                                        >  1:210367/1‑251 (MQ=255)
                                                                                                                              tGCCAATGGTTCAGAAGAAAGTTGATGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGAGTTCTCTACCGCt                                                                                                                       >  2:136658/1‑251 (MQ=255)
                                                                                                                                         cAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCa                                                                                                            <  2:79777/251‑1 (MQ=255)
                                                                                                                                                aaGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCg                                                                                                     <  1:581724/251‑1 (MQ=255)
                                                                                                                                                  gTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGt                                                                                                   <  2:594652/251‑1 (MQ=255)
                                                                                                                                                       tGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAcc                                                                                               >  1:66598/1‑250 (MQ=255)
                                                                                                                                                           ttCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAAccatcc                                                                                           >  2:474120/1‑250 (MQ=255)
                                                                                                                                                            tCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGATAAccatccat                                                                                         >  1:463894/1‑251 (MQ=255)
                                                                                                                                                              tttGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCAtgt                                                                                       >  1:117266/1‑251 (MQ=255)
                                                                                                                                                                  gTAAAGAGCCGCGTAAAGTCGTTAACCCGGACGAAGCTGTAGCAATCGGTGGTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGGCTCTGGGTATCTAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGtgctg                                                                                   <  1:621188/251‑1 (MQ=255)
                                                                                                                                                                    aaaGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAAGCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACCCCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTTAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCa                                                                                 <  2:431298/251‑1 (MQ=255)
                                                                                                                                                                      agagCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATTTGCTGCAgg                                                                               <  1:412604/251‑1 (MQ=255)
                                                                                                                                                                       gagCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAggg                                                                              <  1:182391/251‑1 (MQ=255)
                                                                                                                                                                          ccGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTg                                                                            >  2:21415/1‑250 (MQ=255)
                                                                                                                                                                            gcgTAAAGACGTTAACCCGGTCGAAGCTGTAGCAATCGGTGCTGCGGTTCAGGGTGGTGTTCTGACTGTTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAAc                                                                         <  1:86622/251‑1 (MQ=255)
                                                                                                                                                                            gcgTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAAc                                                                         <  2:171869/251‑1 (MQ=255)
                                                                                                                                                                             cgTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGAGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCACTGTGCTGCAGGGTGAACg                                                                        >  2:173734/1‑251 (MQ=255)
                                                                                                                                                                                           cccGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTg                                                          >  1:119372/1‑251 (MQ=255)
                                                                                                                                                                                             cGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGAGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGAGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGAt                                                        >  2:603376/1‑251 (MQ=255)
                                                                                                                                                                                              ggACGAGGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCCCTATCCCGACCACGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGTTGAACGTAAACGAGCGGCTGAt                                                        >  2:151366/1‑250 (MQ=255)
                                                                                                                                                                                                 cGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGACAAACACCACTATCCCGAACACGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTACCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGCTAAc                                                     >  2:425813/1‑250 (MQ=255)
                                                                                                                                                                                                     gCTGTAGCAATCGGTGCTGCTGTTCAGGGGGGAGTTCTGACTGGTGACGTGAAGGACGTACTGATGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACACAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAtc                                                >  2:622855/1‑251 (MQ=255)
                                                                                                                                                                                                            cAATCGGTGCTGCTGTTACGGACGGAGTTCTGACTGGGGAAGCGAAAGACGTACTGCTGCTGGCCGTTACCCCGCTGTCTCTGGCTATCGAAACCCTGGGCGGTGTGATGACGACGCTGGTCGCGAAAAACCCCCCTATCCCGACCAAGCACAGCCCGATGGTCTCTACCGCTGACGCAAACCCACCTGCGTTCACCATCCCTAGGCTGCAGTGAGAACGTAAACGTGCGGCTGATAACAAAtctattgat                                         >  2:185239/1‑245 (MQ=255)
                                                                                                                                                                                                            cAATCGGTGCTGCTGTGCAGTGTGGTGTTCTGACTGCTGACGTGAATGACATACTGCTGCTGGACGTTACCCCGCTGTCTCTTCGTATCGAAACCATGAGCAGTGTGATGACTACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACATCCATATGTGCTCTACCGCTGAAGGCAACCAGGCTGCTGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAAGCTCTGGGt                                         <  1:598513/251‑1 (MQ=255)
                                                                                                                                                                                                               acGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAATACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAg                                      <  1:286363/250‑1 (MQ=255)
                                                                                                                                                                                                                    gctgctGTTCGTGGTTGTGGTCTGGCTGGTGACGTGAAAGACGTACTGGTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCGACCGCTGAAGACAACCAGTCTGCGGTAACCATACATTTGCTGCAGGGTGAACGTAAACGTGCGGCGGATAACAAAACTCTGTGTCAGTTct                                  <  2:558666/250‑2 (MQ=255)
                                                                                                                                                                                                                       gctgTTCAGGGTGGTGTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGCAAAACACCACTCTCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGCAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCCGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCt                              >  1:79662/1‑251 (MQ=255)
                                                                                                                                                                                                                                ggtggtGTTCTGACTGGTGACGTGAAAGACGTGCTGCTCCTGTACGTTACCGCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGCCGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCCACCGCTGAAGACAACCAGTATGCGGTAACCCTCCATGTGCTGAAGGGTGACCGTAAACGTTCGGCGGACAACCAACCTATGGGTCAGttacacctgtatagtat                     >  2:67996/1‑236 (MQ=255)
                                                                                                                                                                                                                                  tggtgTTCTGACTGGTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGGGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCCGTTCAACCTGGATGGTATc                    >  1:56442/1‑250 (MQ=255)
                                                                                                                                                                                                                                        tCTGACTGGTGACGTGAAAGACGTACTGCTTCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCg              <  1:514502/250‑1 (MQ=255)
                                                                                                                                                                                                                                               ggTGACGTGACAGACGTACTGCTGCTGGACGTTCCCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAACACACCCCTATCCCGACCAAGCACAGCAAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCCTGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAACTCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgcg      >  2:584609/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ggTGACGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACcgc       <  2:162301/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   aCGTGAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAACACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTACCCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACCAATCTCTGGGTCAGTTCAACCTGGATGGTATCAACCCGGCACCGCGCGGc  >  2:27652/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
AACTGGAAAGCCTGGTTGAAGATCTGGTAAACCGTTCCATTGAGCCGCTGAAAGTTGCACTGCAGGACGCTGGCCTGTCCGTATCTGATATCGACGACGTTATCCTCGTTGGTGGTCAGACTCGTATGCCAATGGTTCAGAAGAAAGTTGCTGAGTTCTTTGGTAAAGAGCCGCGTAAAGACGTTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACAACCAGTCTGCGGTAACCATCCATGTGCTGCAGGGTGAACGTAAACGTGCGGCTGATAACAAATCTCTGGGTCAGTTCAACCTAGATGGTATCAACCCGGCACCGCGCGGC  >  NC_000913/13073‑13566

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