Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 674,992 C→A intergenic (‑209/‑26) leuS ← / → ybeL leucine‑‑tRNA ligase/DUF1451 domain‑containing protein YbeL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913674,9920CA98.8% 307.4 / ‑4.0 85intergenic (‑209/‑26)leuS/ybeLleucine‑‑tRNA ligase/DUF1451 domain‑containing protein YbeL
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (46/38);  minor base G (0/1);  total (46/39)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.02e-01

CTATCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACGATGAACAAGGTTGCTCAATATTACCGTGAACTGGTTGCGTCACTGAGCGAACGCCTGCGCAATGGCGAACGTGATATCGACGCACTGGTGGAACAGGCGCGCGAGCGCGTAATAAAAACAGGGGAGTTAACGCGAACCGAGGTCGATGAGCTGACGCGAGCTGTCAGACGTGACCTGGAAGAGTTCGCCATGAGCTATGAAGAGAGCCTGAAAGAAGA  >  NC_000913/674753‑675235
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
cTATCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGATTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGAcc                                                                                                                                                                                                                                          >  2:626254/1‑251 (MQ=255)
cTATCTCTTCCGGGCGGGATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGAAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGCCGTTGAGATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAAAGGCTATTATTATAGAAAGATAACGAcc                                                                                                                                                                                                                                          >  2:423228/1‑251 (MQ=255)
     tctTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCAgggg                                                                                                                                                                                                                                     >  2:466812/1‑250 (MQ=255)
          cGGGCCGTATTGCCCTTGCATGGCAGCCAGTGGTCCTTTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTGTTTCAGATCCGCATAGCATAGCCCAATCGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCTAGTAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGa                                                                                                                                                                                                                                <  1:688307/251‑1 (MQ=255)
              cGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACg                                                                                                                                                                                                                            <  1:248192/251‑1 (MQ=255)
                             aTGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgcgcgg                                                                                                                                                                                                             >  1:602649‑M1/1‑236 (MQ=255)
                                   gCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgg                                                                                                                                                                                                        >  2:441103‑M1/1‑230 (MQ=255)
                                   gCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccga                                                                                                                                                                                                        >  2:417631‑M1/1‑230 (MQ=255)
                                    ccAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACAGTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagtgtgccgaaa                                                                                                                                                                                                      <  1:382594‑M1/251‑23 (MQ=255)
                                     cAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACAGTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATGATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaat                                                                                                                                                                                                     <  1:347334‑M1/251‑24 (MQ=255)
                                                 tttCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtatt                                                                                                                                                                                         >  1:502337‑M1/1‑216 (MQ=255)
                                                 tttCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtat                                                                                                                                                                                          <  1:103118‑M1/250‑35 (MQ=255)
                                                  ttCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtatt                                                                                                                                                                                         <  1:269183‑M1/250‑36 (MQ=255)
                                                     aaTACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccgg                                                                                                                                                                                     >  1:16682‑M1/1‑212 (MQ=255)
                                                     aaTACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccgg                                                                                                                                                                                     >  1:687287‑M1/1‑212 (MQ=255)
                                                      aTACGGCTACAAATGTAGCGTTGAGGTGTTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggc                                                                                                                                                                                    <  2:671163‑M1/251‑41 (MQ=255)
                                                      aTACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggc                                                                                                                                                                                    <  2:136921‑M1/251‑41 (MQ=255)
                                                      aTACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggc                                                                                                                                                                                    <  1:106297‑M1/251‑41 (MQ=255)
                                                      aTACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCAATCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatggattccggc                                                                                                                                                                                    >  2:106297‑M1/1‑211 (MQ=255)
                                                          ggCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTAGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAAAGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagggtgccgaaataagggatggattacggcacg                                                                                                                                                                                 >  2:663478‑M1/1‑207 (MQ=255)
                                                              aCAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGCACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggagagcgtgccgaaatatgggatgtattccggcacgata                                                                                                                                                                              >  2:249693‑M1/1‑203 (MQ=255)
                                                              aCAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgcta                                                                                                                                                                              >  2:565049‑M1/1‑203 (MQ=255)
                                                              aCAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataag                                                                                                                                                                            <  1:1194‑M1/251‑49 (MQ=255)
                                                              aCAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataag                                                                                                                                                                            <  2:668187‑M1/251‑49 (MQ=255)
                                                              aCAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataag                                                                                                                                                                            <  1:574422‑M1/251‑49 (MQ=255)
                                                              aCAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgata                                                                                                                                                                              >  1:181899‑M1/1‑203 (MQ=255)
                                                                aaaTGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgacaagaa                                                                                                                                                                          <  2:235151‑M1/251‑51 (MQ=255)
                                                                      aGCGTTGAGGTGTGTTTTCAGATCCGCATAGCATAGCCCAAACTGCCGCGGCAAAACAGCCTGTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTGAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggatt                                                                                                                                                                    <  2:197475‑M1/251‑57 (MQ=255)
                                                                            gAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGAAGTTGTAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAACGCAAGGAACAAAGAACGTATATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttac                                                                                                                                                              <  1:59522‑M1/251‑63 (MQ=255)
                                                                               gtggtTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAACCGgcatccggtgcgcgtgccgaaatctgggatgtattccggcacgatcagaagggattatttacgtc                                                                                                                                                           >  1:439549‑M1/1‑186 (MQ=255)
                                                                                     tttCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgac                                                                                                                                                     >  1:612988‑M1/1‑180 (MQ=255)
                                                                                      ttCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgac                                                                                                                                                     <  2:322572‑M1/250‑72 (MQ=255)
                                                                                       acaGATCCGCATCGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTTAAATGATGCCCGGATTACTCATACATTAATTTACAGAGGTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcttgccgaaatatgggatgtattccggcatgataagaagggattattgacgtcgctgacg                                                                                                                                                    <  2:72218‑M1/249‑73 (MQ=255)
                                                                                             ccGCATAGCATAGCCCAAACGTCCCCGTCAAAACAGCCTTTCGCGCACCCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggcttatttacgtcgctgacgcagactc                                                                                                                                             >  2:253185‑M1/1‑172 (MQ=255)
                                                                                                      atagcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattntttacgtcgctgacgcagactcatcaacac                                                                                                                                     >  2:390408‑M1/1‑163 (MQ=255)
                                                                                                      atagcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacac                                                                                                                                     >  1:548447‑M1/1‑163 (MQ=255)
                                                                                                        agcCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagc                                                                                                                                  >  2:72562‑M1/1‑161 (MQ=255)
                                                                                                                   ccGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGATAGGATGCCCGGATGACTCATACATTAATTTACAGAGTTTGTGGGCGTAGTAGCAAAGCAAGGGACAAAGAACGTCTTTTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatggggtgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaa                                                                                                                       <  2:46220‑M1/251‑102 (MQ=255)
                                                                                                                       gTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaaca                                                                                                                    >  2:64739‑M1/1‑146 (MQ=255)
                                                                                                                         aaaaacAGCCTTTCGCGCACGCGACGTTGAANTGATGCCCTGATTACTCATACATTAATTTACAGAGTTGGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataataaaggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatg                                                                                                                 <  2:373405‑M1/250‑108 (MQ=255)
                                                                                                                          aaaaCNGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatc                                                                                                                                          <  2:375334‑M1/225‑83 (MQ=255)
                                                                                                                          aaaaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatc                                                                                                                                          >  1:375334‑M1/1‑143 (MQ=255)
                                                                                                                            aaCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccg                                                                                                              >  2:746042‑M1/1‑141 (MQ=255)
                                                                                                                             aCAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgt                                                                                                             >  2:102633‑M1/1‑140 (MQ=255)
                                                                                                                              cAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtc                                                                                                            <  2:570323‑M1/251‑113 (MQ=255)
                                                                                                                                   tttCGCGCACTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcac                                                                                                       >  2:239393‑M1/1‑134 (MQ=255)
                                                                                                                                            cTCGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcgg                                                                                              >  1:145981‑M1/1‑125 (MQ=255)
                                                                                                                                              cGACGTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggc                                                                                             <  1:240462‑M1/250‑128 (MQ=255)
                                                                                                                                              cGACGTTGAAATGATGCACGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggaacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggcc                                                                                            >  1:751886‑M1/1‑123 (MQ=255)
                                                                                                                                                 cgtcgaaatGATGCCCGGATTACTAATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgccgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagt                                                                                         >  1:263395‑M1/5‑120 (MQ=255)
                                                                                                                                                  gTTGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagtt                                                                                        >  2:699620‑M1/1‑119 (MQ=255)
                                                                                                                                                    tGAAATGATGCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgaagcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgc                                                                                      >  1:590754‑M1/1‑117 (MQ=255)
                                                                                                                                                      aaaTGATGCCCGGATTACTCATACATTAATTTACAAAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTAAACGACCCGGGAGATGAAAAGgcatccggtgagcgtgccgaaatatgggatgtattccggctcgataagaagggattatttacgtcgatgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccg                                                                                    <  2:226454‑M1/251‑137 (MQ=255)
                                                                                                                                                        atgatgCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggacgtattccggcacgataagcagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggt                                                                                  >  1:198717‑M1/1‑113 (MQ=255)
                                                                                                                                                         tgatgCCCGGATTACTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggt                                                                                  >  2:401425‑M1/1‑112 (MQ=255)
                                                                                                                                                                  gATTATGCATACAGTAATTATCAGATTTTGTGGGTGTATTAGCAAAGCAAGGAACAAAGATCGTCTATTATTATAGTAAGTTAACGACCCGGGATATGAAACGgcatcgggtgagcgtgccgaaatatgggatgtataccggcacgataagaagggattatttacgtcgctgacgcagattcatcaacacagcagcaaaaccaaacaatgccgtcagcagccacagcggccagttgccggtacgtgcgtat                                                                        <  2:439549‑M1/251‑149 (MQ=255)
                                                                                                                                                                   attactCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaaagccgtcagcacccacagcggccagttgccggtacgtgcttatg                                                                       >  1:463618‑M1/1‑102 (MQ=255)
                                                                                                                                                                    ttactcataCATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatgg                                                                      >  1:22003‑M1/5‑101 (MQ=255)
                                                                                                                                                                     tactcatACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacg                                                                                                                                                    <  1:395472‑M1/169‑73 (MQ=255)
                                                                                                                                                                     tactcatACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCAGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgatgccg                                                                                                                                                    >  2:395472‑M1/4‑100 (MQ=255)
                                                                                                                                                                                 aaTTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggt                                                          <  1:636734‑M1/250‑163 (MQ=255)
                                                                                                                                                                                    ttACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaccacagcagcaaaaccaaacaatgccgtcagcccccacagcggcccgttgccggtacgtgcgtatggtgtgagtccggtggtc                                                      >  2:148163‑M1/1‑85 (MQ=255)
                                                                                                                                                                                      actgagtttGTTGGCGTATTATCAAAGCAAGGAACAAAGAACGGCTATTATTATAGTAAGTTAACGACCCGGGAGAGTAAACGgcatccggtgagcgtgccgaaagatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgtttccggtggtcgg                                                    <  2:518609‑M1/248‑169 (MQ=255)
                                                                                                                                                                                       cAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatggtatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggc                                                   <  2:181899‑M1/251‑170 (MQ=255)
                                                                                                                                                                                             gtgtggGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcac                                              <  1:152590‑M1/249‑175 (MQ=255)
                                                                                                                                                                                                tgGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatctgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgtt                                           >  1:528020‑M1/1‑73 (MQ=255)
                                                                                                                                                                                                       aTTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggtt                                    <  2:291115‑M1/250‑185 (MQ=255)
                                                                                                                                                                                                       aTTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgaatccggtgggcggcgtcacgttagaggtt                                    >  2:155499‑M1/1‑66 (MQ=255)
                                                                                                                                                                                                       aTTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaacaaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggtt                                    >  1:715189‑M1/1‑66 (MQ=255)
                                                                                                                                                                                                          agcaaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaac                                 >  1:743715‑M1/1‑63 (MQ=255)
                                                                                                                                                                                                           gcaaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacac                               <  1:253185‑M1/251‑190 (MQ=255)
                                                                                                                                                                                                             aaagcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacaccc                             >  2:33688‑M1/1‑60 (MQ=255)
                                                                                                                                                                                                                gcaaGGAACAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgc                          >  1:323629‑M1/1‑57 (MQ=255)
                                                                                                                                                                                                                       aCAAAGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaa                   <  1:441103‑M1/251‑202 (MQ=255)
                                                                                                                                                                                                                         aaaGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaac                                                                                                                          <  1:463249‑M1/146‑99 (MQ=255)
                                                                                                                                                                                                                         aaaGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaac                                                                                                                          >  2:463249‑M1/1‑48 (MQ=255)
                                                                                                                                                                                                                           aGAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgc               <  1:215794‑M1/251‑206 (MQ=255)
                                                                                                                                                                                                                            gAACGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcg              >  1:113284‑M1/1‑45 (MQ=255)
                                                                                                                                                                                                                              aCGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcggg            <  1:72562‑M1/251‑209 (MQ=255)
                                                                                                                                                                                                                              aCGTCTATTATTATAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccgggggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcggg            >  1:570797‑M1/1‑43 (MQ=255)
                                                                                                                                                                                                                               cgtttattagtataggaAGTTAACGGCCCGGGAGATTAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcagcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcgggg           <  2:602649‑M1/247‑210 (MQ=255)
                                                                                                                                                                                                                                    attattatAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcgggatcatc      <  2:22003‑M1/251‑215 (MQ=255)
                                                                                                                                                                                                                                    attattatAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcgggatcat       <  1:746042‑M1/250‑214 (MQ=255)
                                                                                                                                                                                                                                     ttattatAGTAAGTTAACAACCCGGGAGATGAAACGgcagccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtggtcggcgtcacgttagtggttaacacctcgcgggtgaactgcgggatcatc      <  2:590754‑M1/250‑215 (MQ=255)
                                                                                                                                                                                                                                      tattatAGTAAGTTAACGACCCGGGAGATGAAACGgcatccggtgagcgtgccgaaatatgggatgtattccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcagcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtcccgggggcggcgtcacgatagtagttaacacctcgcgggtgaagtgcgggaacatcgc    >  2:463341‑M1/1‑35 (MQ=255)
                                                                                                                                                                                                                                         tataGTAAGTTAACGACCCGGGAGATGAAACGgcatcctttgagcgtgccgaaatatgggatgtatgccggcacgataagaagggattatttacgtcgctgacgcagactcatcaacacagcatcaaaaccaaacaatgccgtcagcacccacagcggccagttgccggtacgtgcgtatggtgtgagtccggtgttcggcgtcacgttagtggttaacacctcgcgggtgaactacgggatcatcgcct  <  1:423228‑M1/250‑219 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
CTATCTCTTCCGGGCGGTATTGCTCTTGCATGGCAGCCAGTGGTCCTGTTTTCAATACGGCTACAAATGTAGCGTTGAGGTGGTTTTTCAGATCCGCATAGCATAGCCCAAACGTCCGCGTCAAAACAGCCTTTCGCGCACTCGACGTTGAAATGATGCCCGGATTATTCATACATTAATTTACAGAGTTTGTGGGCGTATTAGCAAAGCAAGGAACAAAGAACGTCTATTATTATAGTCAGTTAACGACCCGGGAGATGAAACGATGAACAAGGTTGCTCAATATTACCGTGAACTGGTTGCGTCACTGAGCGAACGCCTGCGCAATGGCGAACGTGATATCGACGCACTGGTGGAACAGGCGCGCGAGCGCGTAATAAAAACAGGGGAGTTAACGCGAACCGAGGTCGATGAGCTGACGCGAGCTGTCAGACGTGACCTGGAAGAGTTCGCCATGAGCTATGAAGAGAGCCTGAAAGAAGA  >  NC_000913/674753‑675235

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