Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,363 G→A S210S (TCG→TCA)  thrC → threonine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,3630GA93.4% 209.4 / ‑0.1 78S210S (TCG→TCA) thrCthreonine synthase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base A (31/40);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

ATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCG  >  NC_000913/4118‑4609
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
aTTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTaa                                                                                                                                                                                                                                                   <  1:259999/251‑1 (MQ=255)
     gACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTTAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGa                                                                                                                                                                                                                                              <  1:467670/251‑1 (MQ=255)
        cgcgACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTa                                                                                                                                                                                                                                           <  1:190235/251‑1 (MQ=255)
                  ggTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCg                                                                                                                                                                                                                                 >  1:316306/1‑251 (MQ=255)
                   gtgcgaCCGGAGAGACAGTGGCTAATGCTTTCTACGGTTTACCGAATGGGAAAGAGGTTATCCTCTATACACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGTCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCgg                                                                                                                                                                                                                                <  2:120943/246‑2 (MQ=255)
                   gTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGACATATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCTGTTGATGATGA‑AGGGCTGACAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGt                                                                                                                                                                                                                                >  2:364255/1‑251 (MQ=255)
                           ggAGCCGCAGGGGCTCATGCTTTCTACCGTTTACCGAAGGTGAAAGTGGTTATCCTGTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGGACATTGGGCGGCAATATCGAAACTTTGGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGTTGAGGCAGGCGTTTGATGATGA‑AGAGCTTAAAGTGTCGCTGGGGTTAATCACAGCTAACTCGCTTAACATGAGTCGTTTGCTGgc                                                                                                                                                                                                                        <  2:655131/251‑1 (MQ=255)
                               cGGCAGAGGCTCATGCTATCTACGGTTTACCGAATGTGCAAGTGGTTATACTCTATCCACGAGGCAAAATCAGTCCCCTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAAATGTTGCCATAGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCAGAAAGGGGCGCTGTGGTTAAACTCAACTAACTCGATTAACATCAGCTGTTTGCTGgcgac                                                                                                                                                                                                                     >  2:560991/1‑248 (MQ=255)
                                ggCAGTGGCTCATTATTTCTCCGGATTACCGAATGTGAAAGGGGTTATCCTCTATCCAAGCGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCAAAACTGTTGCAATAGACGGCGATTTAGATGACTGTCAGGCGCTGGAGAAGCAGGAGTTTAATGATGA‑CGAGCTGAAAGGGGCGCTGGGGGTAAAATCCGCTAACTCGTTTAACAACAGCAGTTTGCGGGCGCAgg                                                                                                                                                                                                                   >  2:468910/1‑250 (MQ=255)
                                   aGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGAttt                                                                                                                                                                                                                <  1:325572/251‑1 (MQ=255)
                                                   tACGGTTTACCGCATGTGAACGAGGTTATCCTCTATCCACGCGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGCTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGCGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGCTTTGCTCCTACTTTGAAGc                                                                                                                                                                                                >  2:503454/1‑251 (MQ=255)
                                                    aCGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAg                                                                                                                                                                                                 <  1:262686/249‑1 (MQ=255)
                                                    aCGGTTTACAAAAAGTGAAAGTGTTTAGCCTGTATCCACGAGGCAAATTCAATCCACTGCAAGCAATACTGTTCTGTACATTGGGCGGCAATATCGAACGTGATGCCACCGACGTCGATTTCGATGCCTGTCAGGCGCTGGTGGTGCAGGCGTTTGGTGATGA‑AGAACGAAAAGTGGCGCTAGGGTTCAACTCGGATAACTCGATTAAAATCAGCCTATCGCTGGCGCAGATGTGCTAAGACTTTGAAGcc                                                                                                                                                                                               >  2:268457/1‑250 (MQ=255)
                                                       gTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                            >  1:4511/1‑251 (MQ=255)
                                                        tttACCGAATGTGAAAGTGGTTATCCGCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                            <  2:456319/250‑1 (MQ=255)
                                                              gAATGTGAAAGTGGTAATACTCTATCCACGAGGCAAAATCAGTCCAATGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGCAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGAAGGCGTTTGCTGATGA‑AGACCTGAAAGAGGCGCTAGGGTTAAACACGGCTAAATCGATTAAAAACAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTgcgc                                                                                                                                                                                        >  2:593325/1‑248 (MQ=255)
                                                                 tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCCGCTAACTCGATTACCATCAGCCGTTTGCTGGCGCGGATTTGCTACTACTTTGAAGcag                                                                                                                                                                                              >  2:418725/1‑237 (MQ=255)
                                                                 tgtgAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGcag                                                                                                                                                                                              <  1:418725/239‑3 (MQ=255)
                                                                                       ccACGAGGCAAAATCATGCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCt                                                                                                                                                            <  1:114863/251‑1 (MQ=255)
                                                                                               cAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAACGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtct                                                                                                                                                    >  2:279846/1‑251 (MQ=255)
                                                                                               cAAAATCAGGCCCCTGAAAGAAAAACTGGTCTGTACATTGGGAGGCAATATCGAAACTGTTGCCACCGACGGCGATGTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGTCACTAGGGTTAAAATCGGCTAACTCGATGAACATCAGCCGTTTGCTGGCGCAGATTTGCTCCTACTTGGAAGCTGTTGCGCGGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGtct                                                                                                                                                    <  2:51029/251‑1 (MQ=255)
                                                                                                   aTCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                                                <  1:537250/251‑1 (MQ=255)
                                                                                                      aGTCCACTGCAAGAAAAACTGTTCTGGACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGACAGTGGCGCTGGGGTTAAACTCAGCTCACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCATTTGCGCAGCGGCCGCAGAAAGCGCGCACCCACCTGGTTGtcgcgtgtcc                                                                                                                                             >  2:581461/1‑243 (MQ=255)
                                                                                                               cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaa                                                                                                                                     >  1:539105/1‑250 (MQ=255)
                                                                                                               cAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCAGCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGaaa                                                                                                                                    <  1:286116/251‑1 (MQ=255)
                                                                                                                aaGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAAATCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAAc                                                                                                                                   >  1:379782/1‑251 (MQ=255)
                                                                                                                 aGAACAACAGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCTATCAACGGCGCTTTCGATGTCTGACAAGCCCTGGTGAAGCAGACGCTTGATGATGC‑GGAGCTGAAAGGGGCGGTGGGGGTAAACTCAGCTAACTCGATTACCATCAGCGGTGTGCTGGCGCAGGTTGGATACTACTTTGAAGCAGTTGCGTCGCTCCCGCAGGAAGCGCACAACCAGCTGGATGTACCGGTGcccacaggaaact                                                                                                                                  >  2:13155/1‑240 (MQ=255)
                                                                                                                 aGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACt                                                                                                                                  >  2:536611/1‑251 (MQ=255)
                                                                                                                        cTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGCCGGCGATTTCGATGCCTGTAAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCgg                                                                                                                           >  2:34376/1‑250 (MQ=255)
                                                                                                                                tACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGTGTTGAcgg                                                                                                                   >  2:372779/1‑251 (MQ=255)
                                                                                                                                 aCATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                                                     >  1:137070/1‑216 (MQ=255)
                                                                                                                                 aCATTGGGCGGCAATATCGAAACTGTGGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGtc                                                                                                                                                     <  2:137070/216‑1 (MQ=255)
                                                                                                                                  cATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGAcggcg                                                                                                                 >  1:397026/1‑251 (MQ=255)
                                                                                                                                     tGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGt                                                                                                                                                >  1:487639/1‑217 (MQ=255)
                                                                                                                                      gggCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGt                                                                                                              <  1:313348/250‑1 (MQ=255)
                                                                                                                                      gggCGGCAATAACGAACCTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGTAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTc                                                                                                             <  2:122206/251‑1 (MQ=255)
                                                                                                                                          ggcAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGACTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCACACGGAAACTTCGGCGATTTGCCGGCgggggtgtt                                                                                                         >  2:319302/1‑245 (MQ=255)
                                                                                                                                            cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCCGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGCAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGt                                                                                                                                                                                             >  1:645176/1‑165 (MQ=255)
                                                                                                                                            cAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATAAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGtacgctgg                                                                                                       >  1:380436/1‑244 (MQ=255)
                                                                                                                                               tatCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTggggg                                                                                                    >  1:393309/1‑248 (MQ=255)
                                                                                                                                                         gTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCg                                                                                           <  1:104744/250‑1 (MQ=255)
                                                                                                                                                            gCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCt                                                                                       <  1:536611/251‑1 (MQ=255)
                                                                                                                                                             ccatggccggagagttcggggCCTGTCAGGCGCTGTTGAAGCGTGCGTTTGAGGA‑GAGAGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACGCGAGTAACTTCAGCCGTTGGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCTCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGGCACTCGGTCTg                                                                                      <  2:588811/231‑1 (MQ=255)
                                                                                                                                                             ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGGAGGAGTTTGATGATGA‑AGAGCTGAAAGTGGCGATGGGGTTAAACTCAGCTAACTCGATTAAAATCAGCCGTTTGCTGGCGCAGAAGTGCTACTACTTTAAAGAAGTTGCGAAGATGCAGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAACATTCGGCGATTTGACGGCGGGActgtgggcgggagggacaggggtg                                                                                      >  2:659677/1‑230 (MQ=255)
                                                                                                                                                             ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCCGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCCGGtgggctgggggaggaccgtggatgg                                                                                      >  2:514265/1‑227 (MQ=255)
                                                                                                                                                             ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                                                        <  2:452215/201‑1 (MQ=255)
                                                                                                                                                             ccATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCg                                                                                                                                        >  1:452215/1‑201 (MQ=255)
                                                                                                                                                               aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGcc                                                                                    <  1:245566/251‑1 (MQ=255)
                                                                                                                                                               aTCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCGGCTGGCggcgacaacccgcgcggc                                                                                    >  2:312113/1‑234 (MQ=255)
                                                                                                                                                                  gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                            <  2:146957/144‑1 (MQ=255)
                                                                                                                                                                  gACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGtt                                                                                                                                                                                            >  1:146957/1‑144 (MQ=255)
                                                                                                                                                                   aCGGCGATTTCGATTCCTGTCAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTg                                                                                <  1:330595/251‑1 (MQ=255)
                                                                                                                                                                    gggCGATTTGAATGCCTGTCAGGCGTGGGTTAAGCAGGCGTTGGAGGAGGA‑AGAGCGTAAAGTGGCGCTGGGTTTAAACTCAGCTAACGCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGa                                                                               <  2:314301/250‑1 (MQ=255)
                                                                                                                                                                    cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                                               >  1:401063/1‑187 (MQ=255)
                                                                                                                                                                    cGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTGGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTg                                                                                                                                               <  2:401063/187‑1 (MQ=255)
                                                                                                                                                                     ggCGATTTCGATGCCTGTCAGGCGCTGGCGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGTGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAAACAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGaa                                                                              <  1:259989/251‑1 (MQ=255)
                                                                                                                                                                       cGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGATGCTGGCggaggtccaccgtctgccggggaga                                                                             >  1:232509/1‑226 (MQ=255)
                                                                                                                                                                         aTTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGt                                                                          <  1:617040/251‑1 (MQ=255)
                                                                                                                                                                            tCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGttt                                                                        <  1:431307/250‑1 (MQ=255)
                                                                                                                                                                                  ccTGGCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACGCAGCTAAGCCGATGAACATCAGCCGGTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCTATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATtgct                                                                 <  2:397026/251‑1 (MQ=255)
                                                                                                                                                                                       cAGGCGCTGGTGAAGCAGGCGGTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGAc                                                            <  2:316306/251‑1 (MQ=255)
                                                                                                                                                                                                 ggagaccggCGTTTGATGATGA‑AGAGCTGAAAGGGGCGCTGGGGTAAAACTAAGCTAACTCGATTAACAACAGCGGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAGCGTGAAc                                                  <  2:4511/244‑1 (MQ=255)
                                                                                                                                                                                                       aGGCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGTTTAAACTCAGCTAACTCGCTTAACATCGGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCACCCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTa                                                                      <  2:542452/225‑1 (MQ=255)
                                                                                                                                                                                                        ggCGTTTGATGATGA‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCg                                           <  1:268960/251‑1 (MQ=255)
                                                                                                                                                                                                              tgatgatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTTTTTCTGCGACCAACGTGAACGATACCGTGcccc                                     >  1:193019/1‑249 (MQ=255)
                                                                                                                                                                                                                 ggatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTct                                                                                                            <  2:19485/176‑1 (MQ=255)
                                                                                                                                                                                                                  gatga‑AGAGCTGAAAGTGGCGCTGGGGTTAAACCCAGCTAACTCGATTAACAGCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCGGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGttt                                 <  2:193019/251‑1 (MQ=255)
                                                                                                                                                                                                                        agaaCTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGc                            <  1:312113/251‑1 (MQ=255)
                                                                                                                                                                                                                         gagcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                                          <  1:156575/121‑1 (MQ=255)
                                                                                                                                                                                                                         gagcTGAAAGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTgg                                                                                                                                                          >  2:156575/4‑124 (MQ=255)
                                                                                                                                                                                                                         gagcTGAAAGTGGCGCGGGGGTTAACCTCAGCGAACTCGATGAACAGCAGCCGTTTGCTGGCGCAGATTTGGTACCACTTTGAAGCAGTTGCGCAGCTGCCGTAGGAAGCGCGCAAACAGATGGTTGTCTCGGTGCCAAGCGGAAACGTCGGCGATTTTACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCa                           <  2:384063/248‑1 (MQ=255)
                                                                                                                                                                                                                               aaaGTGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACgg                     <  1:36973/251‑1 (MQ=255)
                                                                                                                                                                                                                                   tGGCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTAACTCGGTCTGCCGGTGGAAAGTTTTTTTGTTGCGACCCACGTGAAAGATAACGTGCCCCGTTTTCTGcccgccgggccg                 >  1:428663/1‑240 (MQ=255)
                                                                                                                                                                                                                                    ggCGCTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCTGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGGCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCATGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAg                 <  1:222900/250‑1 (MQ=255)
                                                                                                                                                                                                                                      cgcTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTgg              <  1:610749/251‑1 (MQ=255)
                                                                                                                                                                                                                                      cgcTGGGGTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGGCTGCTGGCGAAGGCACTCGGGCTGCCGGTGGAACGGTTTATTGCTGGGACCAAAGTGAAAGACACCGGGGCACGCTTCCTGCCCGACGGTCAGTgg              >  1:104366/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tggggTTAAACTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCCGGCGAAGTCACCCGGGCCGCCGGTGGAACGTTTTTTTGCTTCGACCAACGGGAAAGAAACCGTGCCAAGTTTTCTGCACGACCGTCCGTgggcc           >  1:587991/3‑248 (MQ=255)
                                                                                                                                                                                                                                                  aCTCAGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCAGTTGCGCAGCTGCCGCAGGAAGCGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCg  <  1:500122/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      
ATTCTGACCGCGACCTCCGGTGATACCGGAGCGGCAGTGGCTCATGCTTTCTACGGTTTACCGAATGTGAAAGTGGTTATCCTCTATCCACGAGGCAAAATCAGTCCACTGCAAGAAAAACTGTTCTGTACATTGGGCGGCAATATCGAAACTGTTGCCATCGACGGCGATTTCGATGCCTGTCAGGCGCTGGTGAAGCAGGCGTTTGATGATGA‑AGAACTGAAAGTGGCGCTAGGGTTAAACTCGGCTAACTCGATTAACATCAGCCGTTTGCTGGCGCAGATTTGCTACTACTTTGAAGCTGTTGCGCAGCTGCCGCAGGAGACGCGCAACCAGCTGGTTGTCTCGGTGCCAAGCGGAAACTTCGGCGATTTGACGGCGGGTCTGCTGGCGAAGTCACTCGGTCTGCCGGTGAAACGTTTTATTGCTGCGACCAACGTGAACGATACCGTGCCACGTTTCCTGCACGACGGTCAGTGGTCACCCAAAGCG  >  NC_000913/4118‑4609

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: ‑