Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,684,443 C→T K392K (AAG→AAA pdeK ← c‑di‑GMP phosphodiesterase PdeK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,684,4430CT98.8% 262.1 / NA 82K392K (AAG→AAApdeKc‑di‑GMP phosphodiesterase PdeK
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (45/35);  minor base G (0/1);  total (45/37)
Fisher's exact test for biased strand distribution p-value = 4.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

TAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGA  >  NC_000913/3684199‑3684684
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
tAACCATCAGCCCACAGCACTCAATGCGAACGATTCAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCt                                                                                                                                                                                                                                             >  1:467556/1‑251 (MQ=255)
    cATCAGCCCACAGCACTCAATACGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAAAAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGGCAACCGCCTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGGTTTTTGCCTTTAt                                                                                                                                                                                                                                         >  2:466213/1‑251 (MQ=255)
          cccACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCg                                                                                                                                                                                                                                     >  1:178263/1‑249 (MQ=255)
          cccACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAg                                                                                                                                                                                                                                   <  1:119385/251‑1 (MQ=255)
                gCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTa                                                                                                                                                                                                                              >  1:565731/1‑250 (MQ=255)
                     cAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGcc                                                                                                                                                                                                                        <  1:455943/251‑1 (MQ=255)
                       aTGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGACATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGCACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGc                                                                                                                                                                                                                      >  1:496621/1‑251 (MQ=255)
                          cgatcgatTAAGCCCTCCGGCAGGTCCCTACTGCCATCCGGTTGCTGGATACGCAGTAACACTGCCGCACTGACCAGTTGACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCATTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTGCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTACATGCCGCAGa                                                                                                                                                                                                                   <  2:376974/251‑1 (MQ=255)
                           gatcgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGa                                                                                                                                                                                                                   >  1:307151/1‑250 (MQ=255)
                              cgatTAAGCCATCCGTCAGGTCCCAACTGCCAGCCGGTTGCTGGATGCGCAGTATCTCTTCCGCACTGACCAGTTTTACGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATGTTCCATTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAGATGCCGCAGAAATa                                                                                                                                                                                                               <  1:442358/251‑1 (MQ=255)
                              cgatTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATa                                                                                                                                                                                                               <  1:124736/251‑1 (MQ=255)
                               gatTAAGCCATCCGGCAGGTCACCACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGCCCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGCTATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTACATCTGATTTTTACCTTTATGGCGAGCGGTACATGCCGCCGaaaacg                                                                                                                                                                                                              >  2:478429/1‑248 (MQ=255)
                                atTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAgc                                                                                                                                                                                                             <  1:303252/251‑1 (MQ=255)
                                   aaGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGcgcg                                                                                                                                                                                                          <  1:65540/251‑1 (MQ=255)
                                     gCCATCCGTCAGGTCCCAACTGCCAGCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCGGTGGCTGTAACCAAATAGCAAACTGATGATTGTCCAGTGCATGAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCGATCGGCTGCGTATCAAAGAACTTAATCTGATTTTAGCCGGTATGGCGAGCGGTAAATGCCGCAGAAATAGAGCGAc                                                                                                                                                                                                        <  2:259429/251‑1 (MQ=255)
                                          ccGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGCCTGTa                                                                                                                                                                                                    >  1:327651/1‑250 (MQ=255)
                                           cGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTACCACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTGACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTaaa                                                                                                                                                                                                  <  1:182071/251‑1 (MQ=255)
                                             gCAGGCCCCCACTGCCATCCGGTTGCTGGATTCGCAGTACCACTTCCGCACTTGCCAGTTTATCGCTGGTCGTTTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTGCCAGTTCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTTGGCGGCTTCCATCTGCTGCGTATCAAAGAACTGAATCTGATTTTTGCCGTGAGGGCGAGCGGTAAATGCCGCAGAAATAGCGCGATTGTAAAga                                                                                                                                                                                                <  2:489014/251‑2 (MQ=255)
                                                  tCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGCATCTGATTTTTGCCTTTATGTCGAGCGGTAAATGACGCAGAAATAGCGCGACTGTAAAGATGtt                                                                                                                                                                                            >  2:471072/1‑250 (MQ=255)
                                                   cccAACGGCCATCCGGTTGCTGGATACGAAGTAACGCTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATGTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTcg                                                                                                                                                                                          <  1:208631/251‑1 (MQ=255)
                                                     aatcTGCCAGCGGGTTGCTGGATACGCAGTAACACTTCCGCACTTACCAGTGTACCGCTGGTCATCTCGACCTTTGGCTGTAACCAAATAGCAAGTTTATGTTGTTCCAGTGCATTAAGGATATCGTGCTCTGCCGTCAACCGTTTCTGGGCGGCTTCCATCTGCTGTGGATCAAAGAACTGAATCTGATTTTTGCCTTTATTGCGAGCTTGAAAGGCCGCAGAAATAGCGCGACTGTAAATCTGTTcggc                                                                                                                                                                                        <  2:28280/248‑1 (MQ=255)
                                                          gCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTga                                                                                                                                                                                   <  2:178263/251‑1 (MQ=255)
                                                            cATCCGGTTGCTGGATACGCAGTAACACTTACGCACTGACCCGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAACAAACTGATGCTTTTCCAGTGCATTAAGGATATCGCTCTCTTACGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAACGAACTGAATATGATTTTTCCCTTTATGGCGAGCGGTAAATGCCTCAGAAATAGCGCGCCTGTAAAGCTGTTcggagggggg                                                                                                                                                                                  >  2:246242/1‑243 (MQ=255)
                                                                cgtttgatGGAGACGCAGTAACTCTTCCGCACAGACCAGTTTACCTCTGGTTATCGCTTCCTGTGGCTGTAACCAAATAGCAACCTGATGATTTGCCAGTGCATTAAGGTTATCGCTACCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATAAAATAACTGAATCTGATTTTGGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCTACTATAAAGCTGTTCGGCGGTGAGAt                                                                                                                                                                                <  2:64845/241‑1 (MQ=255)
                                                                cGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAAACAAAAAGCAAACTGATGATTTTCCAGTGCCTTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGCCTGTAAAGCTGTTCGGCGGTTAGATCGc                                                                                                                                                                             >  2:544866/1‑251 (MQ=255)
                                                                        ggaaaCGCCGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACGTATGATGTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGATCTGAATCTGATTTTTGCCTTTATGGCGAGCGGGAAATGCCGCAGAAATAGCGCGACTGTAAATCTGTTCGGCGGTGAGATCGCCGTAGaa                                                                                                                                                                      <  1:263392/246‑1 (MQ=255)
                                                                             cGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATcgc                                                                                                                                                                >  1:100329/1‑251 (MQ=255)
                                                                                     cacTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGAGCGACTGTAAAGCTGTTCGGCGGAGAGATCGCCGTAGAACATcgccacgcca                                                                                                                                                         >  2:224910/1‑250 (MQ=255)
                                                                                         tCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATg                                                                                                                                                      <  1:127839/249‑1 (MQ=255)
                                                                                           cGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAc                                                                                                                                                  >  1:493552/1‑251 (MQ=255)
                                                                                            gCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGt                                                                                                                                                                                             >  1:124364/1‑207 (MQ=255)
                                                                                            gCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGt                                                                                                                                                                                             <  2:124364/207‑1 (MQ=255)
                                                                                              aCTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAACGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTTAGATCGCCGTAGAACATCGCCACGCCAATGCTaaagg                                                                                                                                               >  2:169558/1‑247 (MQ=255)
                                                                                                tGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGc                                                                                                                                             >  1:283904/1‑251 (MQ=255)
                                                                                                     aGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGcggacg                                                                                                                                        >  1:119220/1‑251 (MQ=255)
                                                                                                          aCCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGGAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCCCGCCAATGCTACAGTGCGGACGGAGtt                                                                                                                                   >  1:188172/1‑251 (MQ=255)
                                                                                                             gCTGGTCATCTCGACCGGTGGCTGAAAACAAAGAGCAAACTGATGATTTTCCAGAGCATTAAGGATATCGCTCTCTTCCGTCACCCGCTTCTGGGCGGATTCCATATGCTGCGGCTAAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGCAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTTAGATCGCCGTAGAACATCGCCACGCCACAGCTACacgacggactgagttg                                                                                                                                  >  2:157049/1‑234 (MQ=255)
                                                                                                                ggTCATCTCGACCTGGGGCTGTAACCAAATAGCAAACTGATTGTTTTCCAGTGCATTAAGGATATCGCTTTCTTCCGCCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATa                                                                                                                              <  1:64415/250‑1 (MQ=255)
                                                                                                                  tCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGt                                                                                                                           >  1:148725/1‑251 (MQ=255)
                                                                                                                   cATCTCGACCTGTGGCTGTAACCAAATAGCAAACTTATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTTCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCTCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGt                                                                                                                           <  1:213647/250‑1 (MQ=255)
                                                                                                                   cATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGtt                                                                                                                          >  1:435782/1‑251 (MQ=255)
                                                                                                                       tcGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGGGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGGCGGAGTTGCATACGTTCact                                                                                                                      >  2:229134/1‑249 (MQ=255)
                                                                                                                         gACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCg                                                                                                                    >  1:528955/1‑251 (MQ=255)
                                                                                                                           ccTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAAt                                                                                                                                                       >  1:116559/1‑214 (MQ=255)
                                                                                                                            cTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGc                                                                                                                  >  1:168836/1‑250 (MQ=255)
                                                                                                                                  cTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGACTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGGTCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGgcgcc                                                                                                           >  2:302219/1‑250 (MQ=255)
                                                                                                                                    ttAACCAAATAGCAACCTTATGATTTTCCAGTGCATTAAGGATACCGCTCTTTTCCGTCATCCGCTGCTGGGCGTCTTCCATCCGCATCGGATCAAAGAACGGAATCTTATTTTTGCCATTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCTGCGGTGAGATCGCCGTAGAACATCGCAACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGcgctc                                                                                                          <  2:467556/249‑1 (MQ=255)
                                                                                                                                    gTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTAAGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGAGCGCCGTAGAAAATCGCCACGCCAAAGGTCCAGTGCGGAAGGAGTTGAATAccatcactcggaagggacccg                                                                                                         >  2:276331/1‑231 (MQ=255)
                                                                                                                                          aaaTAGCAAACTGATGATTTTCCCGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTAAAGTGGGGACGGAGTTGAATACGTTCAATCAGCAGGGGCTCGCCCATg                                                                                                   >  2:71949/1‑251 (MQ=255)
                                                                                                                                              aGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGAt                                                                                                                                                                                <  1:439820/170‑1 (MQ=255)
                                                                                                                                              aGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGGGAGAt                                                                                                                                                                                >  2:439820/1‑170 (MQ=255)
                                                                                                                                                 atacggatgatTTACCAGTACATTTAGGATATCGCTCTCGTCCCTCCACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGATCTGAATCTGATTTTTGCCTTTATGGCTAGCGGTAAATGCCTCAGAAATAGCGCGACGGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTGCAATCGGCAGGCGCTCGCTCATGATAGTGa                                                                                            <  2:153546/246‑1 (MQ=255)
                                                                                                                                                       atgatTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGAACtt                                                                                      >  1:335077/1‑251 (MQ=255)
                                                                                                                                                       aggatTTTCCAGTTCTTTAAGGATATCGCTCTCTTCCTTCAACCGCTGCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACtt                                                                                      <  2:437323/249‑1 (MQ=255)
                                                                                                                                                                cAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCACCACGCCAATGCTACAGTGAGGACGGAGTTGAATACGTTCAATCGGCAGGAGCTCGCTCATGATAAGGAGCACTTGCTTACCt                                                                              >  2:210109/1‑250 (MQ=255)
                                                                                                                                                                         aaGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGCTCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCACCGCCAATGCTACAGTGCGGACGGGGTTGAATACGTTCAATCGGCAGGCGCACGCCCATGATAGTGAGCACTTGCTGACCCAAAGTGAttt                                                                    >  2:49534/1‑250 (MQ=255)
                                                                                                                                                                          agggtatCGCTCTCTTCCGGCACCCGCATCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTTAATCTGATTTTTGCCGTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGTTGTGCGGCGGTGAGATCGCCGTATAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATTCGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGAtttc                                                                   <  2:396487/247‑3 (MQ=255)
                                                                                                                                                                             atatCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACtt                                                                                      >  1:512346/1‑229 (MQ=255)
                                                                                                                                                                             atatCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTg                                                                >  1:516957/1‑251 (MQ=255)
                                                                                                                                                                                 tgccggctTCCGTAAACCGCCTCGGGGGGGCTTACAGCTGCTGCTGACCAAAGAATTGAATCTGATTGTTTCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGTGACTGTAAAGCTGTTCGGCGGTGAGATCGCGGTAGAGCATCGCCGCGCCAATTCTACAGTGCGGACGGAGTGGAATACGTTCAATCGGCAGGGGCGCGCTGATGATTGTGATCACTGGCTGACCTAAGGTTATTGGGTGCCa                                                             <  2:245667/244‑1 (MQ=255)
                                                                                                                                                                                   ctctctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACa                                                                                                                                                 >  2:194458/1‑164 (MQ=255)
                                                                                                                                                                                   ctctctTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACa                                                                                                                                                 <  1:194458/164‑1 (MQ=255)
                                                                                                                                                                                     ccctTCCGACAACCGCTTCGGGGCGGCTTCCATCTGCTGCGGATCAAAGAACGGAATCTGATTTTTGCCTTTATGGCGGGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTGCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCGAAGGTGATTGCGTGCCACGGt                                                         <  2:266209/248‑1 (MQ=255)
                                                                                                                                                                                           cGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTa                                                  <  1:468066/251‑1 (MQ=255)
                                                                                                                                                                                                         gggaggcttccagcgtctgaggATCAAAGAACTGAATCTGATTTTTGCCTTGATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGGAAAGCTGTGCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACTCCGTTGGCAAt                                     <  2:493552/230‑1 (MQ=255)
                                                                                                                                                                                                          ggcggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATg                                    >  1:279753/1‑250 (MQ=255)
                                                                                                                                                                                                          ggcggcTTCCATATGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGa                                   <  1:220116/251‑1 (MQ=255)
                                                                                                                                                                                                            cggcTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGAc                                  >  1:274708/1‑250 (MQ=255)
                                                                                                                                                                                                                  ccATCTGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGGCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTACTGATAGTGAGCACTTGCTGACCTAAGGGGATTGCGGGCCACGGGTCCTGGACACCGTTGGCAATGACAGCaa                             >  2:454317/1‑249 (MQ=255)
                                                                                                                                                                                                                  ccATCGGCTGCGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCAAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTCCCTGTACACCGTTGGCAATGACAGCAAAg                           <  1:466213/251‑1 (MQ=255)
                                                                                                                                                                                                                     tCTGCTGCAGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGAGTGCCACGGTTCCTGTACACCGTTGGCAATGCCAGCAAAGTc                         >  1:776/1‑250 (MQ=255)
                                                                                                                                                                                                                     gctgctgCGGATTAAAGAACGGAATCTGAGTTCTGCATTTATGGCTAGCGGTAATTGCCTCATAAATAGCGCGCCTTTAAAGCTTCTCTGCGGTGAGACCGTCGTAGAACATCGCCACGCCAATGCTCCAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCGCGCTCATGATAGTGAGCACTTGCTGACCGAAGGTGATTTCGTGACACGGTTCCTGTACTCCGTTGGCAATGACATCAAAGTc                         <  2:505432/249‑1 (MQ=255)
                                                                                                                                                                                                                       tgctgcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATc                                                                                                                                                                               <  1:106118/98‑1 (MQ=255)
                                                                                                                                                                                                                       tgctgcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATc                                                                                                                                                                               >  2:106118/1‑98 (MQ=255)
                                                                                                                                                                                                                       tgctgcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCAt                       >  1:378958/1‑250 (MQ=255)
                                                                                                                                                                                                                          tgcGGATCAAAGAACTTAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTTAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                    <  1:518246/250‑1 (MQ=255)
                                                                                                                                                                                                                           gcGGATCAAAGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAAc                    <  1:49534/249‑1 (MQ=255)
                                                                                                                                                                                                                                   aaGAACTGAATCTGATTTTTGCCTTTATGGCGAGCGGGAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAAAATCGCCACGCCAATGCTAAAGTGCGGACGGAGTTGAATAAGTTCAATCGGCAGGCGCTCGCTCATGAGAATGAGCACTTACTGCACTAAGGGGATTGCCAGCCACGGTTCCTGTACCACGTTGGCAATGACAggcaatacctaacccataatccg          >  2:529232/1‑229 (MQ=255)
                                                                                                                                                                                                                                     gAACTGAAGCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACAGCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTgcg        <  1:357555/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aaCTGAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGt                                                         <  1:185363/201‑1 (MQ=255)
                                                                                                                                                                                                                                          gAATCTGATTTTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAgg   >  1:469356/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aTCTGATTTTTGCCTTTATGACGAGCGGTAAATGCCGCAGGAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAAGTCGCCGTAGAACATCGCCACTCCAATGATACAGGGAGGAAGGAGTTGAATACGTTCAATCGGCAGGCGCTCGATCATAATAGTGAGCACTTGCTGACCTAAAGAGATTTCGTACCACGGTTCCTGTACACCGATGGCAATGACAGCAACGTCTGAACAGCTAATCTGCGCGAGGa  >  2:248090/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                 
TAACCATCAGCCCACAGCACTCAATGCGATCGATTAAGCCATCCGGCAGGTCCCAACTGCCATCCGGTTGCTGGATACGCAGTAACACTTCCGCACTGACCAGTTTACCGCTGGTCATCTCGACCTGTGGCTGTAACCAAATAGCAAACTGATGATTTTCCAGTGCATTAAGGATATCGCTCTCTTCCGTCAACCGCTTCTGGGCGGCTTCCATCTGCTGCGGATCAAAGAACTGAATCTGATTCTTGCCTTTATGGCGAGCGGTAAATGCCGCAGAAATAGCGCGACTGTAAAGCTGTTCGGCGGTGAGATCGCCGTAGAACATCGCCACGCCAATGCTACAGTGCGGACGGAGTTGAATACGTTCAATCGGCAGGCGCTCGCTCATGATAGTGAGCACTTGCTGACCTAAGGTGATTGCGTGCCACGGTTCCTGTACACCGTTGGCAATGACAGCAAAGTCATAACCGCTAATCTGCGCGAGGA  >  NC_000913/3684199‑3684684

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