Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,945 A→G 17.2% S138P (TCA→CCA)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9450AG17.2% 123.5 / 2.5 71S138P (TCA→CCA) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base A (24/32);  new base G (7/6);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTG  >  NC_000913/989709‑990182
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
tCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagc                                                                                                                                                                                                                                               <  2:66118/239‑1 (MQ=255)
  ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagcga                                                                                                                                                                                                                                             <  1:956577/239‑1 (MQ=255)
  ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagcga                                                                                                                                                                                                                                             <  1:1240384/239‑1 (MQ=255)
   cATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagcgat                                                                                                                                                                                                                                            <  2:4742/239‑1 (MQ=255)
     tCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagcgatag                                                                                                                                                                                                                                          <  1:532502/239‑1 (MQ=255)
     tCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagcgatag                                                                                                                                                                                                                                          <  1:675118/239‑1 (MQ=255)
       ggCGTGAGGGAAAGCCGACGCGGCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                        <  1:800664/239‑1 (MQ=255)
         cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGagcgatggcga                                                                                                                                                                                                                                       <  2:1266383/238‑1 (MQ=255)
                      cGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCAc                                                                                                                                                                                                                          <  1:333845/238‑1 (MQ=255)
                                gATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCAACAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCTATAGCGATGTACCATTTCAc                                                                                                                                                                                                                          <  1:190770/228‑1 (MQ=255)
                                          gCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGGGGGCGGTAGCGATGTACCATTTCACTGATAACC‑GCCAGCaaaagga                                                                                                                                                                                                    >  1:1288675/1‑235 (MQ=255)
                                          gCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGGCGAGAAGTGGACTAATTTGCTTTCCTGCGGGGCGGGGGGGTGGGCC‑ACC‑TCGGCCAGCGGAGGGCGATAGGCATGGAACATTTCTCTGATACAC‑Gccccggcactta                                                                                                                                                                                                    >  2:190770/1‑229 (MQ=255)
                                              gTGCCCACAAACAAGGATTCGTGTGGAAGACGCTTAACGAGGGTTTTTTGTACTTCGCGAAAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATACCAAGACCGGCGGTTACCCCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCGGCGGTGGGCGATAGCGATGTACCATTCCACGGATAACC‑GCCAGCAAAGGAACtt                                                                                                                                                                                                 <  2:194986/238‑1 (MQ=255)
                                                 cccACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGACGACTAATTTTCTTTCCCGCGTGTCGGGGGGTTGCGCA‑ACC‑TCGGCCTGCGGGGAGCGGTAGCGATGTACCCTTTCACTGATAACC‑GCCCGCAAAGGAACTTcccc                                                                                                                                                                                             >  1:964644/1‑238 (MQ=255)
                                                   cacaAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCCGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCaga                                                                                                                                                                                           >  1:213677/1‑239 (MQ=255)
                                                      aaaccgggtctctggTTGCAGATGCTTACGTGTGGTTTGTGGTACACCGCGAGAAAACCGGCGCCGGGTGCCAAAACCCTCCAGATGGAACCGCGACATACCAAGCCCGGCGTTTCGCGCCGGGAGGGCGACTATTTTGCTTTCCCGCGTGTGGAGGGCTTGCAAA‑ACG‑TCGGCCTGCGGTGGGCGACGGAGATGTTCCATTGCCCGGATAGCC‑GCCAGCAAAGGAACTTCCCAGAGGa                                                                                                                                                                                        <  2:708095/225‑1 (MQ=255)
                                                             gctgccggttgagggcgtttacagATGGTTTCTTGTCCTTCTGAAGAAAAACGGGGACGGGTGCCAAAACACATAAGAGGGAGGCGTGACATATCAAGCCGGGGGGTTAACGCGGAGAAGGCAACAAATTGGCTTTCCGGCGTGTAGAGGGCTTGCGCA‑ACG‑TCGGCCTGGGGTGAGCAATAGCGATGTCCCATTTCACTGAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACtt                                                                                                                                                                                  <  2:689661/215‑1 (MQ=255)
                                                                 ccTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAACCGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGCCCGGCGGTTACCGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGaa                                                                                                                                                                                                    <  1:583861/216‑1 (MQ=255)
                                                                  cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGACGACTAATTTGCTTTCCCACGTGGCGAGGGGGTTCGCA‑ACG‑TCGGCCCGCGGGGAGCGATAGGGATGTACCATTTCCCTGATAACCACCCCGCAAGGAAATTCCCCAAAGGATGA‑TTCAcac                                                                                                                                                                             >  1:1358054/1‑236 (MQ=255)
                                                                    gTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGCGAAGTCGACTAATTTGCTTTTCAGCGGGGCGGGGGGTTGCGCC‑ACG‑TCGGGCTGCGGTGAGCGGTAGCGATGTACCATTTTACTGATAACC‑GCCCGCAAAGGAACTTTCCAGAGGATGACTTCACGcccc                                                                                                                                                                          >  1:501511/1‑237 (MQ=255)
                                                                           aCGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGGGGGCTTGCGCA‑ACG‑TCGGCCCGCGGGGAGCGATAGGGGTGTAACCTTTCACTGATAACC‑CCCAGCAAAAGAACTTTCCCGAgggcggctcttatacacatctgt                                                                                                                                                                    >  2:583861/1‑217 (MQ=255)
                                                                           aCGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAACGCCGAGAAGACGACTAATTTTGTTTCAAGCGTCTCGGGGGGTTGCGCC‑ACG‑TCGGCCTGCGGCGGGCGATAGCGCTGGAACATTTCCCAGGAAACA‑GCCAGCAAAAGGAAGTCCCAGGGGATGACTTTACCCCTCCGGccc                                                                                                                                                                    >  2:1241610/1‑237 (MQ=255)
                                                                           aCGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGGTAAAGCCGAGAAGTCGACTAATTTTCTTTTCAGCGTGTCGAGGGATTGCGCA‑ACG‑TACGCCTGCAGTGAGTGGTAGCGGTGGACCATTTTCCTGAGAACA‑CCCAGCCAAAGAAGTTTCCAGGGGAGAGCTTCACAcccacggccg                                                                                                                                                                    >  2:605814/1‑230 (MQ=255)
                                                                                tAACGATGGTTTCTTGTCCTCCGCGAGAAAACCGGGGACGGGTGCAAAAATCCATCAGGGGAAACCGCGACATATCAGGACCGGGGGTTACCGCCGAGAGGTCGACTAATTGGCTTTCCAGCGTGTGGAGGGCTTGCGCA‑ACG‑TCGGCCGGCGGTGAGCGATAGCGATGTACCATTTCACTGAAAGCC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTaaa                                                                                                                                                              <  2:1165655/239‑1 (MQ=255)
                                                                                  acggtggtttcgtagcgtcaccgAAAAAAACGTGGACGGGTCCCAAAACCCAGCAGTGGGAACCGCGACTAATCAGCGCCGGCGTTTACCCCCGAGATGTCGACTAATTTGCTTCCCGGGTGGTCGAGGGCTTGCCCA‑GCG‑TCGCCCGGGGGGGAGCGATAGGGATGTACCATTTCACTGAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAAc                                                                                                                                                             <  2:501511/217‑1 (MQ=255)
                                                                                       ggTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAATCGACTAAATTGCTTTCCAGCGTGTCGGGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGCTGGACCATTTCACTGATAAAC‑GCCAGCAAAGGAAATTCCCAGAGGATGACTTCACGCCACGGGCCGCCTAAACGaaaa                                                                                                                                                        >  1:1206731/1‑236 (MQ=255)
                                                                                         tttCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGAATAATTTTCTTTCCAGCGGGTCGGGGGCTTGCGCA‑ACG‑TCGGCCCGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAAGAACATCCCCGAGGATGACTTCACGCCCACGGCCGCTTaaaagcaaaaa                                                                                                                                                      >  2:1274740/1‑230 (MQ=255)
                                                                                           tCTTGTACTTCGCGAGAAAACCGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCACGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGCGGGGTTGCCCA‑AAGTTCGGCC‑ACGGGTAGCGATACGGGTGTACCATTTCTATGACAAAC‑GCCCCCAAAAGAACTTCCCCGACGAGAGCCTCCCGCCCCGGGCCggctacactaaactacc                                                                                                                                                    >  2:1326521/1‑222 (MQ=255)
                                                                                            cTTGTACTTCGCGAAAAAAACGGCGACGGGTGCAAAAACCCATCAGATGAAAGCGCGACATATCAAGACCGGCGTTTACCGCCGAGAATTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGGAATAGCGATGTCCCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGc                                                                                                                                                  <  1:71086/239‑1 (MQ=255)
                                                                                             ttGTACTTCGCGAGAAAAAGGGCGACGGGGGCCAAAATCCATCAGAGGGAAGCGCGACATATCAAGCCCGGGGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑t                                                                                                                                                <  1:38408/239‑1 (MQ=255)
                                                                                              gggtctgggggagaAAAACGGGGGCGGGTGCAAAAACCCATAAAGGGGAAGCGTGACATAACAAGCCCGGCGGTTACCGCCGAGAAGTCGCCTTATTTGCTTCCCAGGGTGTGGGGGGTTTGCGCA‑CCG‑CCGGCCGGCGGGGGGCGATAGCGATGTCCCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑tt                                                                                                                                               <  2:1162615/229‑1 (MQ=255)
                                                                                                   ttCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCCAGAAATCGACTAATTTGCTTTCCACCGTGTCGGGGGCTTGCGCA‑ACG‑CCGGCCCGCGGTGGGCGGTAGCGATGTTCCATTTCACCGATAAAC‑CCCCGCAAAGGAACTTACCACAGGAATACTTCACGCCACGGGGCCCTTAAACGAATATCCACC‑CTGccct                                                                                                                                          >  2:719420/1‑237 (MQ=255)
                                                                                                   tccgcgAAAAAAACGGGGGGGGGTCCAAAAATCCATAAGAGGAAAGCGCGACATAACAAGCCGGGCGGTTACCGCGGAGAATCCGACTAATTTGCTTTCCAGCGTGTAGAGGGTTGGCGAA‑CCG‑TCGCCCTGGGGGGAGAGATAGGGATGTCCCATTTCCCTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCAt                                                                                                                                          <  2:1318142/237‑1 (MQ=255)
                                                                                                     cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGTTGGAAGCGCGACATATCAAGACCGGCGGCTAAAGCCGAGAAATCTGCTAATTTGCTTTCCAGCGTGTAGAGGGCCTGCGGA‑ACC‑GCGGGCCGCGGTGTGCGAAAGCGATGTGCCCTTTTTCTCACAACA‑GCCCGCAAAGGGACTTGCCCGAGGATGTCATTACGCCACCCGCCCCCTTAAAGAATAACACAC‑TTGCCCTTa                                                                                                                                        >  2:160932/1‑239 (MQ=255)
                                                                                                     cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAAGCCGAGAAGTCGGCAAATTTGCTTTTCAAAGGGTCGCGGGATTGCGCC‑ACG‑TGGGCGTGCGGGGGGCGAAAGGGAGGTACCATTTCAATGATAACC‑GCCAGCAAAGGAACTTCCCCGAGGATGAGTTCACGCCACGGGCGCCGTAAACGaaactcccacatgcgatta                                                                                                                                         >  2:588592/1‑222 (MQ=255)
                                                                                                     cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGGCGAGGGCTTGCGCA‑ACG‑GCGGCCTGCGGTGAGCGGTAGCGAGGTACCATTTCACTGATAACC‑GCCCGCCAAAGAAATTCCCAGAGGATGACTTCACGCCCCGGGCCGCTTAAACGAATATACAGC‑TTTCCCTTa                                                                                                                                        >  2:761944/1‑239 (MQ=255)
                                                                                                     cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGGCGACTAATTTGCTTTCCAGCGTGTCGCGGGCTTTCGCC‑ACG‑GCGGCCTGCGGTGAGCGGTAGCGGTGTACCCTTTTACTGGTAACC‑GCCAGCAAAGGAACTTCCCCGAGGATGACTTTCCCCCACGGGCCGCCTAAAAGAATATCCCGC‑CTgtccttt                                                                                                                                        >  2:1374695/1‑233 (MQ=255)
                                                                                                     cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACCCCCAGAAGTCGACAAATTTGCTTTCAAACGTGCCGTAGGCTCGCGCA‑ACG‑TCCGCCTGCGCGCAGGCATAGCGGAGTTCCATTTCACCGATAACA‑CCCCGCCAAAGAACTTCCCCAAGGAGAACTTCCCCCCCCCGGCCGCGTAAAACAACATCCACCATTGCCAtt                                                                                                                                         >  2:738154/1‑239 (MQ=255)
                                                                                                     cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAACGCCGGGCAGACGACTAATTTTATTTTCCGCGTGTCGGGGGCGTGCGCACACG‑‑CCGCCTGCGGTGGGCGCTAGCGATGTAGCATTTTTCTGATAAGC‑GCCCGCAAAGGAACTTCCCTGAGGCTGACTTCACGCCCCGAGCCGCTGTAAAGAAAATCCAAC‑CTGGCATTa                                                                                                                                        >  2:1051256/1‑239 (MQ=255)
                                                                                                                        ggTGCCAAACACAAACAGAAGGAAGCGCGACATAGAAAGCCCGGCGGTTAACGCTGAGAATCCGACTAATTGCCTTTCCAGCGTGGGGAGGGCTGGGGAA‑ACG‑TGGGCCTGGGGAAAGGGATAGCAATGTCCCATTTCGCTAAAACCC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGAc                                                                                                                       <  2:1218247/237‑1 (MQ=255)
                                                                                                                                                 cgcgACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTCGAGGGGTTGCGCA‑ACG‑GCGGCCTGCGGTGAGGGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGGACTTCCCAGAGGATGACTTCACGCCACGGGGCGCTTAAACCAATATc                                                                                                                                                      >  1:141332/1‑182 (MQ=255)
                                                                                                                                                    gACATATCAAGACCGGCGGTTTACGCCGGGAAGTCGACTAAATTGCTTTCCAGCGTGTCGAGGGCCTTCGCC‑ACG‑CCGGCCCGCGGTGAGCGAAAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGGACTTCCCCGAGGATGACTTCACCCCACGGGCCGCCTCAACGAATATCCCAC‑TTTCCATTATCGTGTTGCACGGGGGCTTTTTCCGGGGTAAAACGGGACTcgcgga                                                                                          >  2:1049122/1‑236 (MQ=255)
                                                                                                                                                    aaagtaccAAGACCGGGGGTTAACGCCGAGACGCCCACTATTTTGCTTTCGAGCGTGCCGAGGTCTTGCGAA‑ACG‑GCGCCCTGCGGTGTGGGATAGGGATGTCCCATTTCACTAATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                         <  1:801429/232‑1 (MQ=255)
                                                                                                                                                         atCAAGACCGGCGGTTAACGCCGAGAAGTCGACTATTTTGCTTCCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGAt                                                                                                                                                                                       <  1:1260426/141‑1 (MQ=255)
                                                                                                                                                                   ggggTTAGCGCGGAGAGGGCGACTAATTTGCTTTCCAGCGTGGGGGGGGCTTGCGCA‑CCG‑TCGGCGTGCGGGGGGGAATAGCGATGTACCATTTCACTGAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                           <  2:979529/236‑1 (MQ=255)
                                                                                                                                                                   ggggTTAACGCCAAGAGGCGAACTAATTTTCTTTCATGTTTTTCGGGGTCTTGCTAA‑ACG‑TCGGCCGGCGGTGAGGGATAGCGATGTACCATTCCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                           <  2:1206731/236‑1 (MQ=255)
                                                                                                                                                                   gcgttaaaCGCCGGAAATTGGACTATTTTGCTTCCAGGGGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGGGGTGAGCGATGGCGATGTACCATTTCACTAATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                           <  1:647794/232‑1 (MQ=255)
                                                                                                                                                                   gcgtTTAACGCCGGTAAGTCGACTAATTTGCTTTCCGGCGTGTCGGGGCCTTGCGCA‑ACG‑TCGGCGTGCGTTGAGCGATAGCGATGTACCATTTCCCGGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                           <  1:767947/238‑1 (MQ=255)
                                                                                                                                                                   gcggTTAACGCCGCGAAGTCGACTAATTTGCTTTCCAGCGGGTCGAGGGCTTGCGCC‑ACG‑TCGGCCTGCGGTGGGCGATAGCGATGTACCATTTCACTGATAAAC‑GCCAGCAAAGGAACTTCCCAGAGGGTGACTTCACGCCACCGGCCGCCTAAACGAAAATCCCGC‑TTGCCATTATCGTTGGACACGGGGACTTGTTTCGGGTTAAAAAGGAACTCTCGCAACAAGTGAAGATaatg                                                                          >  2:530600/1‑236 (MQ=255)
                                                                                                                                                                   cggttaaaggcggaaaggTCGACTAATTTGCTTTCCAGCGGGCGGGGGGCTGGCGAA‑AGG‑TGGGCCTGCGGGGAGCGATAGCGATGTCCCTTTCCACTGAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                           <  2:579767/221‑1 (MQ=255)
                                                                                                                                                                    gggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                          <  1:1374695/237‑1 (MQ=255)
                                                                                                                                                                    cggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTCGAGGGCTTTCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTTACTTATAACC‑GCCAGCAAAAGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCCTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGGGACTTGTTTCCGGTTTAAGCGGGACTCCCCTAACCAGCTAATATTGTCg                                                                         >  2:80632/1‑239 (MQ=255)
                                                                                                                                                                              ggaaaagcggattagtttGCTTTCCGGGGTGTGGGGGCCTGGCGCA‑ACG‑GCGGCGTGCGGGGAGCGATAGCGATGTACCATTTCACTAATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaag                                                                <  2:1288675/223‑1 (MQ=255)
                                                                                                                                                                                 gaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaag                                                            >  2:362211/1‑239 (MQ=255)
                                                                                                                                                                                  aaGTCGACTATTTGCCTTTCCGGCGTTTCGAGGGGTTTGGCA‑ACG‑TGGGTCGGTGGTGAGCGATAGCGATGTCCCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                           <  1:400397/239‑1 (MQ=255)
                                                                                                                                                                                  aaGTCGAATAATTGTGTTTCCGGCGTTTCGGGGGCTTGCGAA‑CCC‑CCGGTCTGCGGTTAGCGATAGCGATGTACCATTTCACTGAAAGCC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                           <  1:918647/239‑1 (MQ=255)
                                                                                                                                                                                   aGTCTACTAATTTGCTTCCCAGCGTGTCGAGGGTTGGCGAA‑ACG‑CCGGCCGGCGTTGAGCGATAGCGATGTCCCATTTCCCGAATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                          <  1:294285/239‑1 (MQ=255)
                                                                                                                                                                                        aCTAATTTGCTTTCCAGCGCGGCGAGGGCGTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGGCCCATTTCACTGATAACC‑GCCAGCTAAGGAACTTCCCAGAGGATGACTTCACGCCACGGTCCGCTTAAACGAATATCCGGC‑CTGCCTTTATCCTCGGACACGGTGAGTTGGTCCGGGTTAAACCGGAACTTGACTCACCAGTTTAGGAAGTGGGCCTTGACAGAAggcagggc                                                     >  2:1111617/1‑239 (MQ=255)
                                                                                                                                                                                                 cTTTCCAGCGGGTCGGGGGCTTGCAAA‑ACG‑TGGGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCa                                            <  1:605814/239‑1 (MQ=255)
                                                                                                                                                                                                  tttCCGGCGGGTCGAGGGCTGGCGCA‑ACG‑TCGCCCGGGGGTGAGCGATAGCGATGTACCATTTCACGGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCa                                            <  1:80632/238‑1 (MQ=255)
                                                                                                                                                                                                      cAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                       <  1:362211/239‑1 (MQ=255)
                                                                                                                                                                                                          gtgtAGAGCGCTTGGGCA‑GCG‑CGGCCCTGGGGTGAGGGATAGAGATGTACCTTTCCCCGGAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                                    <  1:1324783/238‑1 (MQ=255)
                                                                                                                                                                                                           tgtCGGGGTCTTGCGCA‑AGG‑TCGGCCTGCGGTGAGCGAGAGGGATGTGCCATTCCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtc                                   <  2:650234/238‑1 (MQ=255)
                                                                                                                                                                                                           tgtCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtca                                  >  2:594722/1‑239 (MQ=255)
                                                                                                                                                                                                                     ttGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCa                        >  2:291113/1‑239 (MQ=255)
                                                                                                                                                                                                                       gcgcA‑ACG‑TCGGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAgg                      <  1:325606/239‑1 (MQ=255)
                                                                                                                                                                                                                       gcgcA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                       >  1:634120/1‑238 (MQ=255)
                                                                                                                                                                                                                         gcA‑ACG‑CCGGCCTGCGGTGGGCGATAGCGATGTACCATTCCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGgcg                    <  1:1049122/239‑1 (MQ=255)
                                                                                                                                                                                                                              cG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGA‑TTGCAATTATCGTTGGACACGGTGACTTGTTCCGGGCTAAAGCTGAACTCGCGTAACCAGATAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCAATGATATTTATCATCCTGCATgctccgg                >  2:1278089/1‑234 (MQ=255)
                                                                                                                                                                                                                               g‑TCGGCCTGCGGTGAGCTATAGCAATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGt               <  1:849248/239‑1 (MQ=255)
                                                                                                                                                                                                                                  cggccggcGGTGAGCGATACCGAGGTACCATTACACGAAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTAAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtg              <  1:1036121/238‑1 (MQ=255)
                                                                                                                                                                                                                                      cggCGGAGGGCGATGGCGTTATCCCTTTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa         <  2:1170241/237‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGGGGTGAGCGATAGCGATGTACCATTTCACTAATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa         <  1:1241610/239‑1 (MQ=255)
                                                                                                                                                                                                                                      cTGCGGTGAGCGATAGCGATGTCCCATTTCACTGAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgctgc          <  1:761944/238‑1 (MQ=255)
                                                                                                                                                                                                                                       tGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa         >  1:1011007/1‑238 (MQ=255)
                                                                                                                                                                                                                                         cGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                         <  1:1234413/156‑1 (MQ=255)
                                                                                                                                                                                                                                         cGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                         >  2:1234413/1‑156 (MQ=255)
                                                                                                                                                                                                                                         cGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc       >  2:1190366/1‑238 (MQ=255)
                                                                                                                                                                                                                                             gAGCGATAGCGATGTACCATTGCACAGAAAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGACTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTg  <  1:1051256/239‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                               
TCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCA‑ACG‑TCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACC‑GCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGC‑TTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATCGCCTG  >  NC_000913/989709‑990182

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