Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,587,531 G→A L195F (CTC→TTC)  ugpQ ← glycerophosphodiester phosphodiesterase UgpQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,587,5310GA98.7% 261.5 / NA 78L195F (CTC→TTC) ugpQglycerophosphodiester phosphodiesterase UgpQ
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (39/38);  minor base C (1/0);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCG  >  NC_000913/3587293‑3587775
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
gcaTATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGAtt                                                                                                                                                                                                                                          <  2:444782/251‑1 (MQ=255)
      gctgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATg                                                                                                                                                                                                                                    <  2:143607/251‑1 (MQ=255)
       ctgcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATg                                                                                                                                                                                                                                    >  2:398098/1‑250 (MQ=255)
         gcCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAAt                                                                                                                                                                                                                                 <  1:46790/251‑1 (MQ=255)
                 ggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGc                                                                                                                                                                                                                         >  2:581503/1‑251 (MQ=255)
                 ggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGc                                                                                                                                                                                                                         <  1:269951/251‑1 (MQ=255)
                 ggagaAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACGATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCTGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGc                                                                                                                                                                                                                         <  2:203074/251‑1 (MQ=255)
                               aTGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAAATCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGAAGCCCAGCcgcgcg                                                                                                                                                                                                           >  1:281160/1‑251 (MQ=255)
                                        cAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTc                                                                                                                                                                                                   <  1:45349/250‑1 (MQ=255)
                                           cccAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCATCCCAGCCGCGCGGTCAGTTCgcgc                                                                                                                                                                                               >  1:25251/1‑251 (MQ=255)
                                                    ccgcAGCATGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTACGGCGTCTTTCAACTGCATCACTAGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCAGCGCGGTCAGTTCGCGCCAGtcgtcg                                                                                                                                                                                      >  1:128335/1‑251 (MQ=255)
                                                        aGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGcc                                                                                                                                                                                  >  2:581336/1‑251 (MQ=255)
                                                               cATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATc                                                                                                                                                                           >  1:212620/1‑251 (MQ=255)
                                                                   ccGTTGAAAACTCTTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACAAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGc                                                                                                                                                                       >  1:415436/1‑251 (MQ=255)
                                                                   ccGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAGTCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAg                                                                                                                                                                        <  1:582550/250‑1 (MQ=255)
                                                                    cGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGc                                                                                                                                                                       >  2:144848/1‑250 (MQ=255)
                                                                         aaaaCTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACaaa                                                                                                                                                                 <  2:142528/251‑1 (MQ=255)
                                                                          aaaCTATAGGGCCGTAAAGTGCGGCCCAATCCCGACAATAGCATCGGTACAGATGCAATCAACAACCCAGCGCAGCAACTATGATGCGAGATGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGCGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCCACAAAc                                                                                                                                                                >  2:69530/1‑251 (MQ=255)
                                                                          aaaCTAGTGGGCCGTAAAGGTCGGACCAATCACGTCAATCGCATCGGTACAGAGGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAAGCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGAGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGGGCGCGCCATTCGTCGCGCCACTCATCCAGCAACAAAc                                                                                                                                                                <  2:422542/251‑1 (MQ=255)
                                                                            aCTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACcg                                                                                                                                                              >  1:470023/1‑251 (MQ=255)
                                                                            aCTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAAATGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCCTCCAGCACCAAACcg                                                                                                                                                              >  2:347952/1‑251 (MQ=255)
                                                                                     ccGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGt                                                                                                                                                     <  2:25251/251‑1 (MQ=255)
                                                                                        tAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCAACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCCCTCATCCAGCAACAAACCGCGCGGCAGTTcc                                                                                                                                                  >  2:261042/1‑251 (MQ=255)
                                                                                         aaaGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCg                                                                                                                                                 <  2:270315/251‑1 (MQ=255)
                                                                                                ggACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACAACCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGAGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTgccgc                                                                                                                                          >  2:240886/1‑251 (MQ=255)
                                                                                                  aCCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCCTCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAACCCGCGCGGCAGCTCCGGTgccgccc                                                                                                                                        >  2:158546/1‑250 (MQ=255)
                                                                                                    cAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCtgt                                                                                                                                      <  1:581336/251‑1 (MQ=255)
                                                                                                         aCGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGc                                                                                                                                 <  2:457066/251‑1 (MQ=255)
                                                                                                           gTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGc                                                                                                                                 >  2:167881/1‑249 (MQ=255)
                                                                                                               aTCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTc                                                                                                                           >  1:225340/1‑251 (MQ=255)
                                                                                                                      cGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGc                                                                                                                     <  1:265047/250‑1 (MQ=255)
                                                                                                                       ggTACAGATGCAATCCACACCCCAGCGCAGCCACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCAt                                                                                                                   >  1:480333/1‑251 (MQ=255)
                                                                                                                       ggTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCAt                                                                                                                   <  1:144848/251‑1 (MQ=255)
                                                                                                                        gTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAAGCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATc                                                                                                                  <  2:460263/251‑1 (MQ=255)
                                                                                                                               tGCAATCCACCCCCCAGCGCAGCAACTCTGATGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAAtctc                                                                                                            >  1:178155/1‑250 (MQ=255)
                                                                                                                                        cacCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGAc                                                                                                  >  2:121200/1‑251 (MQ=255)
                                                                                                                                          ccccAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGAcag                                                                                                <  1:229410/251‑1 (MQ=255)
                                                                                                                                            ccAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGAcagc                                                                                               >  1:473567/1‑250 (MQ=255)
                                                                                                                                              aGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACagcagc                                                                                            <  1:458346/251‑1 (MQ=255)
                                                                                                                                               gcgcAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGGTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAgcg                                                                                           <  2:520497/251‑1 (MQ=255)
                                                                                                                                                cgcAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAgcgg                                                                                          >  1:55851/1‑251 (MQ=255)
                                                                                                                                                   agcaTCTCTGCTGCGCGCTGGGGTTTAGTGACGGTATAAACCAGAATCCGCAGTTCTGCGTCGTTCAACAGCATCACTCGCGCGGTGTGGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGGTCGCGCCAGTCGTCGCGCCGCTCATCCAGCAACAAACCGCGCGGCAGTGCCGGTGCCGCCTGTGGTGCAGCTTCTAAAGCATCAAGCTCAAACGATGTCAGCAgcggcgg                                                                                       <  1:231116/251‑1 (MQ=255)
                                                                                                                                                   agcaACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAgcggcgg                                                                                       >  1:554461/1‑251 (MQ=255)
                                                                                                                                                      aaCTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGt                                                                                    >  1:50438/1‑251 (MQ=255)
                                                                                                                                                            gctgcGCGCTGTGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCTCTTTATTGAATAACGTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCAt                                                                                 <  2:564225/248‑1 (MQ=255)
                                                                                                                                                                cgcgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGcc                                                                          <  2:321046/251‑1 (MQ=255)
                                                                                                                                                                cgcgcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCGCAAACTATGACAGCAGCGGCGGCGTCATACCGGcc                                                                          <  2:244473/251‑1 (MQ=255)
                                                                                                                                                                   gcTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCcac                                                                       >  2:181522/1‑251 (MQ=255)
                                                                                                                                                                      ggggTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCcaca                                                                      <  2:288111/249‑1 (MQ=255)
                                                                                                                                                                      ggggTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGc                                                                    >  1:62307/1‑251 (MQ=255)
                                                                                                                                                                         gTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGCTTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTcc                                                                 >  2:426946/1‑250 (MQ=255)
                                                                                                                                                                           ttATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTcgcg                                                               >  1:105888/1‑251 (MQ=255)
                                                                                                                                                                                gACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTgcc                                                           <  2:62307/250‑1 (MQ=255)
                                                                                                                                                                                  cGGTATAGACCAGAAGCCGCAGTCCGGCGGCTTGCAACGGCAGCAACCGCGCTTTATTGAATAACTGATGATGGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGGTCGCGCCAGTCGTCGCGCCGCTCATCCGGCAACAAACCGCGCGGCAGTTCCGGGGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCGCACTGCTCGCGTgccgcc                                                        <  2:73841/251‑1 (MQ=255)
                                                                                                                                                                                   ggTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAAAGATGACAGCAGCGGCGGCGTCATACCGGCCCACACCTAGCGAgccgcc                                                        >  1:511311/1‑250 (MQ=255)
                                                                                                                                                                                     tataAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAgc                                                     >  1:188757/1‑251 (MQ=255)
                                                                                                                                                                                      ataAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCAGCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCGCGCGTGCCGCCAgc                                                     <  1:66327/250‑1 (MQ=255)
                                                                                                                                                                                        aagccAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCATCCGCGCGGTCAGTTCGCGCAAGTCGTCGCGCCACTCATCCCGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCATCTAAAGCCTCAATCTCAAACGATGACAGCAGCGGCGGCGTCATAACGGCCCACAGCTAGCGAGCCGCCAGCGcc                                                  >  2:444611/4‑251 (MQ=255)
                                                                                                                                                                                          aCCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGccac                                                <  2:310739/251‑1 (MQ=255)
                                                                                                                                                                                             aGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGACTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCCGATTCTAAAGCATCAATCTCAACCGATGCCAGCAGCGGCGGCGTCATCCCGGCCCACAGCTCGCGTGCCGCCCGCCccaccct                                             >  1:202260/1‑249 (MQ=255)
                                                                                                                                                                                                aTCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCAtttt                                          <  1:138398/251‑1 (MQ=255)
                                                                                                                                                                                                  ccGCAGTCCGGCGTCGTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTg                                         <  2:178113/250‑1 (MQ=255)
                                                                                                                                                                                                   cGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGcc                                       >  2:216957/1‑251 (MQ=255)
                                                                                                                                                                                                     cAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGCTTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACACACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCCGCGGCGGCGTCCTACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCg                                     >  2:74367/1‑251 (MQ=255)
                                                                                                                                                                                                         ccGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTaa                                 <  1:216957/251‑1 (MQ=255)
                                                                                                                                                                                                         ccGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGACGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCCCAGCTCGCGTGCCGCCAGCGCAACCATTTTGCCCGTTa                                  >  1:69330/1‑250 (MQ=255)
                                                                                                                                                                                                           ggCGTCTTTCACCTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTAGCGCCAGTCGTCGCGCCACTCATCCAGGAACAAACCGCGCGGCAGTTCCGGTGACGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGCCAGAAGCGGCGGCGGCATACCGGCCCACAGCTCGCGTGCAGCCAGCGCCACCCTTTTGCCCGTTCATg                               >  2:73024/1‑251 (MQ=255)
                                                                                                                                                                                                                 tttCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCCTTTTGCCCGTTAATGGCccgg                         >  2:139742/1‑251 (MQ=255)
                                                                                                                                                                                                                   tCAACTGCATCACTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCccggt                        <  1:398098/250‑1 (MQ=255)
                                                                                                                                                                                                                              aCTCGCGCTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCATTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTggg            <  2:388287/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cgcgCTTTATTGAATAACTTATGATTGAGATGAATAGAGCCGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGGCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGt           <  1:294585/249‑1 (MQ=255)
                                                                                                                                                                                                                                  gcgcTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCCGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTg        <  2:470023/251‑1 (MQ=255)
                                                                                                                                                                                                                                   cgcTTTATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGa       <  1:347952/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tttATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGAttt    >  1:218157/1‑251 (MQ=255)
                                                                                                                                                                                                                                       ttATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGAAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGAtt     <  1:35201/249‑1 (MQ=255)
                                                                                                                                                                                                                                       ttATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTAGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCAAGCGCCACCATTTTGCCCGTTAATGGCCCGGGGCCGGTGGTGGGTTTGAGTTc   >  2:475310/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tATTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTc   <  1:451288/250‑1 (MQ=255)
                                                                                                                                                                                                                                         aTTGAATAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCCGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCg  >  1:500515/1‑250 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                    
GCATATGCTGCCGGAGAGGAGAAACGGCGACATGCTGAATCAGCCCAGCACGCCGCAGCCTGACATTCCGTTGAAAACTATTGGGCCGTAAAGTTCGGACCAATCACGTCAATCGCATCGGTACAGATGCAATCCACACCCCAGCGCAGCAACTCTGCTGCGCGCTGGGGTTTATTGACGGTATAAACCAGAATCCGCAGTCCGGCGTCTTTCAACTGCATCACTCGCGCTTTATTGAGTAACTTATGATTGAGATGAATAGAGACGCAGCCCAGCCGCGCGGTCAGTTCGCGCCAGTCGTCGCGCCACTCATCCAGCAACAAACCGCGCGGCAGTTCCGGTGCCGCCTGTTGTGCAGCTTCTAAAGCATCAATCTCAAACGATGACAGCAGCGGCGGCGTCATACCGGCCCACAGCTCGCGTGCCGCCAGCGCCACCATTTTGCCCGTTAATGGCCCGGTGCCGGTGGTGGGTTTGATTTCG  >  NC_000913/3587293‑3587775

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