Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 989,943 T→G 22.4% S138S (TCA→TCC pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9430TG22.4% 152.6 / 4.9 83S138S (TCA→TCCpncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (20/40);  new base G (13/10);  total (33/50)
Fisher's exact test for biased strand distribution p-value = 7.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989707‑990177
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
gtTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                                                             <  1:52679/238‑1 (MQ=255)
gtTCCCAGCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                                                             <  2:1099742/238‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATTGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  2:110450/239‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCTCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  1:950384/239‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  1:611882/239‑1 (MQ=255)
    ccATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  1:529084/239‑1 (MQ=255)
    ccATCCGCTGGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTTCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcga                                                                                                                                                                                                                                        <  2:961014/239‑1 (MQ=255)
     cATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgat                                                                                                                                                                                                                                       <  2:635802/239‑1 (MQ=255)
       tCGGCGTGAGGGAAAGCCGACGCACCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatag                                                                                                                                                                                                                                     <  1:249339/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  2:1023546/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  2:211184/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  2:1220634/239‑1 (MQ=255)
         ggCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGagcgatagcg                                                                                                                                                                                                                                   <  1:1075988/239‑1 (MQ=255)
           cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCGAGAAAGTGACCAATTTTCTTTCCAACGGGTCGGGGGCGTTGGCAAACTCGGCCTGCGGTGagcggtagcgat                                                                                                                                                                                                                                 >  1:414676/1‑239 (MQ=255)
           cGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGCTAACACCCAGAAGACGAGTAAATTGCTTTTCAGCGCGGAGGGGGCGTGCGCACCGCCGCGCCGCCGGGatcggtagcgag                                                                                                                                                                                                                                 >  2:548900/1‑238 (MQ=255)
            gTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATg                                                                                                                                                                                                                                <  1:182974/239‑1 (MQ=255)
                gggAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTAc                                                                                                                                                                                                                             <  1:992390/238‑1 (MQ=255)
                     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACCAATTTGCTTTCCAGCGTGTCGGAGGGCTGGGCCAACTCGGCCTGCGGTGAGCGCTAGCGCTGTACCActt                                                                                                                                                                                                                        >  2:897698/1‑238 (MQ=255)
                     aGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTAGACTAATTTGCTTTCCAGCGTGTCGAGGGGTTGCGCCAAGTCGGCCTGCGGGGAGCGAGAGCGATGTAcccattt                                                                                                                                                                                                                       >  1:148373/1‑234 (MQ=255)
                                    tCGTATTTGCTGGTGCCCACAAACCAGGATTCCTGTGGCAGACGTTTAACTATGGTTTCTTTTACTACGAGAAAAAAGCGGCGACGGGTGCCAAAACCCATCAGATGGAAGCGCGACATATCAAGCCCGGGGGTTAACGCCGAGAGGGCGACTAATTTGCTTCCCAGCGTGTCGAGGGTTGGCGCAACGCCGGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGcc                                                                                                                                                                                                         <  2:270343/237‑1 (MQ=255)
                                            gcggtggCCCACAAACCAGGATTCCGGTTGCAGACGTTTAACGATGTTTTCTTTTACTTCGCGAAAAAAAGGGCGACGGGTCCCAAAATCCACCAGATGGAAGCCCGACATACAAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGaa                                                                                                                                                                                               <  2:148373/233‑1 (MQ=255)
                                             cTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCCGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGGGATGGACCATTTTACCGATAACCGCCAGCAAAAGa                                                                                                                                                                                                >  1:1209956/1‑237 (MQ=255)
                                              tGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGGGGGCTTGCGCCACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTTACTGATAACCCCCCGCACAAGAACt                                                                                                                                                                                             >  2:428772/1‑239 (MQ=255)
                                               ggTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTGCTTTCCAGCGTGTCGGGGGGGTGCGCAAACTCGGCCTGCGGGGGGCGATAGCGCTGTCGCAATTTACTGATAAACGCCCACCAAGGGAACt                                                                                                                                                                                             >  2:44251/1‑239 (MQ=255)
                                                  gCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGGACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTACGCCGACAAATCGACTAATTTTCTTTCCCGCCTGGTGAGGGCTTAAGCAAAGTCCGCCTGCGGTGAGCGAAAAGGGTTGACCCTTTTCCTGAGACACGCCAGCAAAGGGAATTcgc                                                                                                                                                                                         >  2:649585/1‑237 (MQ=255)
                                                  gCCCACAAACCAGGATTCCTGTGGCAGACGCTTAACGATGTTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAAGCCATCAGATGGAAGCGCGACATATCAAGCCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTcc                                                                                                                                                                                          <  1:1212552/238‑1 (MQ=255)
                                                       caAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTTCAACGGGGCGAGGGCCTGCGCAACGCCGGCCCGCGGTGAGCGATAACGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCCGAgg                                                                                                                                                                                    >  1:136552/1‑239 (MQ=255)
                                                         aaCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCAAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTACCGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGAt                                                                                                                                                                                  <  1:428772/239‑1 (MQ=255)
                                                          aCCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAAGTGACTAATTTGCTTTTCAACCGGTCGAGGGGCTGCGCAAACGTGGGCCGCGGGGGGTGGAAGAGGTGTACCATTTAACTGATAACCCCCAGGAAAGAAACATCCCAGAGGaag                                                                                                                                                                                 >  2:934169/1‑237 (MQ=255)
                                                               gATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTTCAACGTGTCGAGGGCTTGCGCAACGTCGGGCCGCGGGGGGCGATAGCGGTTTACCATTTCACCGATAAACGCCCGCAAAGGAAATTTCCCGAGGGTGACtt                                                                                                                                                                             >  2:519363/1‑238 (MQ=255)
                                                                aTTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGGGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                            <  1:484931/238‑1 (MQ=255)
                                                                 tttccgTTGCAGGCCGCTTAAGGATGGTTTCTTCTACGTCGCGAAAAAACCGGCGACGGGTGCAAAAAGCCATAAGATGGAAGGCCGAAATATCAAGACCGGGGTTTAACGCCGAAAAGTCGATTAATTTGTTTCCCAGCGTGTCGAGGGCTTGCGCACCGTCGGCCTGCGGTGGGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTc                                                                                                                                                                            <  2:784158/233‑1 (MQ=255)
                                                                    cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTTCTTTTCAGCGGGTCGAGGGCTTGCGCAACGGCGGCCTGCGGTGAGCGGTAACGATGTAACATTTTACTTATAAACGCCCGCAAAGGAACTTCCCCGAGGGTGACTTTACGc                                                                                                                                                                        >  1:532900/1‑238 (MQ=255)
                                                                    cTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCCAGAAGTCGACTCATTTTCTTTCCCGCGGGTCGGGGGGTTGCGCAACGCCCGCCCGCGGGGGGGGGTAGCGGTGTACCATTTTCATGATAAACCCCCGCAAAAGGACATCCCAGGGGGAGACtttccgc                                                                                                                                                                        >  2:182688/1‑233 (MQ=255)
                                                                                           ttttCTTGTACTGCGGAAGAAAAAGGGCGAGGGGTGCCAAAATCCACAGAAGGGAAGCGCGACATATCAAGACGGGCGTTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGAAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAatat                                                                                                                                                  <  1:967773/237‑1 (MQ=255)
                                                                                                ttgtgcggcgaaagaAAAAGGGGGAGGGGTCCCAAAACCAAGAAAGGGGAGGCGCGAAATATCAAGACGGGCGGTAAACACCGAGAGGCCGACAATTTGCTTTTCGAGCGTGTCGAGGGCTTGCAAAACGGCGGCCGGGGGTGAGCGATAGCGATGAACTATTCCACAGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCt                                                                                                                                            <  2:1185405/227‑1 (MQ=255)
                                                                                                 tGTCCTTGGGGAAAAAAGCGGGGACGGGGCCAAAAACCCATCAGATGGAAGCCGGACATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTGGCTTTCCAGCGTGTCGAGGGCTTGCGAACCGTCGGCCTGCGGGGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCtt                                                                                                                                           <  2:961553/239‑1 (MQ=255)
                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTTAAGCCGAGAAGTCGACTAATTTGCTTTTCAGCGGGTCGAGGGGTTGCGCAACGGCGGCCCGCGGGGAGCGATAGCGGTGTACCATTTCACCGATAACCCCCAGCAAAGGAACTTCCCAGAGGGTGACTTTACGCCCCGGGCCGCCTAAAcagaaatacagatttccctta                                                                                                                                    >  2:1117046/1‑219 (MQ=255)
                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTCGGGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCAATGATAACCGCCCGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGc                                                                                                                                                                 >  1:1235500/1‑210 (MQ=255)
                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCCGCGGTGGGCGATAGCGATGTAACATTTCACTGCTAAACGCCAGCAAAAGAACTTCCCCCAGGATGACTTCACGCCACCGGCCGCTTAAACGCATATCCAACTTTCCCttt                                                                                                                                    >  1:1001106/1‑238 (MQ=255)
                                                                                                        cgcgAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGGACAACAAATTGCTTTACAGAGGGTCGGAGGGCTGCGCGAAATCGGGCTTCAGGGGGCGATAGACATTTCACATTTCCTCGCAAAACGCCCGCCAAGGGACTTACCAGAAGAAGGGTTCCAGCCACCGGGCGCCTAAACCGATATTCCACTTGCCACTa                                                                                                                                    >  2:210005/1‑239 (MQ=255)
                                                                                                                aaccGGCGAGGGGGGCAAAAACCCATCAGATGGAAGCGCGACATACCAAGACCGGCGGTAAACGCCGAGAAGTCGACTAATTTGCTTCCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGGGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGa                                                                                                                            <  1:296234/236‑1 (MQ=255)
                                                                                                                   cgcgggcggtccccAAAACCCATCAGAGGGAACCGCAACATAGCAAGCCCGGCTGTTAACGCAGAAAGGTCGATTATTTTGTTTCCCAGCGGGTCGAGGCCTTGCGCAACGTCGGCCTGGGGTGAGCGATAGCGATGAACCATTTCACTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGacac                                                                                                                         <  2:672391/227‑1 (MQ=255)
                                                                                                                              gcaaaaacacATCAGATGGAAGCGCGACATATCAAGCCCGGCGTTTAACGCCGAGAAGTCGACTAATTTGCTTCCCGGCGTGTCGAGGGCTTGCGCAACGTCGGCCGGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                             <  2:532900/230‑1 (MQ=255)
                                                                                                                              gacAAAACCCATAAAGTGGAAGGCCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTCTTTTCCAGCGTGTCGGGGGTTGGCGCACCGCCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                             <  2:849192/237‑1 (MQ=255)
                                                                                                                               ccAAAATCCATCAGAGAAAAGCGGAAAATATAAGGACCGGCTGTTAACGCGAATAAGTCGACTATTTGGCTTCCCGGCGTGTGTGGGGCTTGCGCACCGTGCGCCGGCGGTGGTCGAGAGCGAGTTCCCATTTCACTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                             <  2:943213/238‑1 (MQ=255)
                                                                                                                               ccAAAACCCATCAGATGGAAGCGCGACATATCAAGGCCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGtt                                                                                                             <  1:519363/238‑1 (MQ=255)
                                                                                                                                         tCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCCTTTCACTGATAACCGCCCGCCAAGGAACTTCCCAAGAGATGACTTCACGCCAACGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCATTGGACACAGGGGCTTGTTCCGGGGTaaa                                                                                                   >  2:572585/1‑238 (MQ=255)
                                                                                                                                              aTGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTg                                                                                                                                                                                                                                             <  2:950375/98‑1 (MQ=255)
                                                                                                                                                ggAAGCGCGAAAGATCAAGACCGGGGGTTAACGCCGAGAAGCCGACTAATTTGCTTCCCAGCGTGTCGAGGTCTTGCCCAACGTCGGCCGGGGGTGAGCAATAGCGATGTACCTTTCCCCTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGaa                                                                                             <  2:1209956/238‑1 (MQ=255)
                                                                                                                                                        aaCATATCAAGCCCGGGGGTTAACGCCGGGAGGTCGACTATTTTGTTTCCCAGCGTGTCGGGGGTTTGCGCAACGGCGGCCTGCGGGGGGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                    <  2:289057/238‑1 (MQ=255)
                                                                                                                                                         aCATATCAAGCCCGGCGGTTAACGCGGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa                                                                                    <  1:44251/238‑1 (MQ=255)
                                                                                                                                                            tatCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAAATTTCTTTCCACCGTGTCGAGGGCTTGCGCCACGTCGGCCTGCGGTGAGCGATAGCGATGTACCCTTTCACTGATAACCGCCCGCAAACGCACTTCCCCGAGGATGACTTCCCCCCACCGGCCGCTTAAACGAATATCCCGCCTGCCCTTATTCTTGTACACGGTGGCTTGTTTCCGGGTTAAACGGACAACGCGGAACCac                                                                                >  2:689178/1‑238 (MQ=255)
                                                                                                                                                               cAAGCCCGGGGTTTAACGTCGAAACGTCGGCTAGTTGGCTTCCCAGGGTGTGGAGGGTTTGCCCAAGGTCGGCGCGGGGGGACGGATAGCGTGGTCCCATTTCTCGAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTa                                                                             <  2:1226803/239‑1 (MQ=255)
                                                                                                                                                                      gggggTAAACGCCGGGAAGGCGAATAATTTGCTTTCCAGCTGGTCGAGGGCATGCCAACCGTCGGCTGGCGGTGACGAATAGCGATGTACAATTGCCCTGAGAACCGCCAGCAAAGGAACTGCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCGTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGt                                                                      <  2:157560/236‑1 (MQ=255)
                                                                                                                                                                       gcggTTTAAGCCGAGAAGTCGACTAATTTGCTTTCCAGCGGGTCGGGGGCTTGCGCAACGTCGGCCTGCGGGGAGCGATAGCGATGTACCCTTTTACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCCCGGGCCGCTTAAACGAAAATCCAGCTTGCCATTATCGTTGGACACGGGGGATTGTTCCGGGGTAAAGCGGAACTCCGGGAAACAGTTaaaaaaaac                                                                     >  2:686493/1‑232 (MQ=255)
                                                                                                                                                                       gcggTTAACGCCGAGAAGTCGACTAATTTTCTTTCCAGCGTGTCGGGGGCTTGCGCAACGGCGGCCTGCGGGGAGCGAGAGCGGTGTACCATTTTACTGATAAACGCCCGCCAAGGAACTTCCCCGAGGGTGACTTTCCCCCACCGGCCGCTTAAACGAAAATACACCTTGCCCTTATCGTTTGACACGGTGACCTGTTTCGGGTTCAAGCGGAACTCCCGTAAACCGTTACGATAGt                                                                      >  2:631283/1‑238 (MQ=255)
                                                                                                                                                                        cggTTACCGCCGAGAAGTCGACTAATTTGCTTCCCAGGGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGAGATTTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                     <  1:740722/238‑1 (MQ=255)
                                                                                                                                                                        cggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCGGCGTTGCGTGTGCTTGCGAACCGGCTGCGTGCGGGGATCGATAGCGATGTCCCATTTCCCTAAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTc                                                                     <  2:136552/238‑1 (MQ=255)
                                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                        >  2:379034/1‑239 (MQ=255)
                                                                                                                                                                                    agaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTaaagaaa                                                        >  2:116645/1‑239 (MQ=255)
                                                                                                                                                                                      aaGTCGACTAATTGGCTTTCCAGCGTGTCGGGGGCTTGCCAACCGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                      <  1:901034/239‑1 (MQ=255)
                                                                                                                                                                                       agtcatagaaTTTGTTTCCCACTGTGCCGAGGGTTTGCGCCCCGTGGGCCTGCGGGGAGAGTTAGTGATGTACCATTGCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAgg                                                      <  1:1158402/230‑1 (MQ=255)
                                                                                                                                                                                        gTCGGCTAATTTGCTTTCGTTGGTGGGGAGGTCTTGCGCATGGCCGCCCTGCGTGGAGGGAGAGCGATGTCCTTTCTCATTGATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggc                                                     <  2:379445/238‑1 (MQ=255)
                                                                                                                                                                                         tCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcag                                                   >  2:1218794/1‑239 (MQ=255)
                                                                                                                                                                                            aCTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAggcagggc                                                >  1:204556/1‑239 (MQ=255)
                                                                                                                                                                                                   ttttctttccaggcttcggggggTTTGCGCACCGTCGCCCTGCGGTGCGCGATATCAATGGACCATTCCACTGAAAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGaaa                                           <  1:1023389/219‑1 (MQ=255)
                                                                                                                                                                                                    tgcttccgggCGTGTGGGGGGCTTGCGAACGGCGGGCGAGCGGTGGGGGATGGCGGTGTTCCATTTCGCTAATACCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGc                                         <  2:20286/230‑1 (MQ=255)
                                                                                                                                                                                                        ttCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGcc                                        >  1:707889/1‑236 (MQ=255)
                                                                                                                                                                                                           cAGTGTGTCGGGGGTTTCGGCAATCGCGGCCTGCGTGGGCCAATAGCTATGTACCATTTCACTTATATCCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGat                                  <  1:865714/239‑1 (MQ=255)
                                                                                                                                                                                                               gtgtCGAGGGTTTGCGCAACGCCGGCCGGCGGGGAGCGATAGCGATGTACAATTTCACAGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                               <  1:416752/238‑1 (MQ=255)
                                                                                                                                                                                                               gtggCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCCCTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATAtt                               <  1:1218794/238‑1 (MQ=255)
                                                                                                                                                                                                                tgtCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATtca                             >  1:819315/1‑239 (MQ=255)
                                                                                                                                                                                                                            tGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  1:731040/1‑239 (MQ=255)
                                                                                                                                                                                                                             gcgcAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  2:249316/1‑238 (MQ=255)
                                                                                                                                                                                                                             gcgcAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAg                  >  1:1091819/1‑238 (MQ=255)
                                                                                                                                                                                                                                          gCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgct       >  2:334793/1‑238 (MQ=255)
                                                                                                                                                                                                                                            cggCGGTGAGCGATAGCGATGTCCCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:1152536/237‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGGGGTGAGCGATACCGTTGTGCCATTTCACTAATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:680435/239‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGCGGTGGGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCa    <  1:182688/239‑1 (MQ=255)
                                                                                                                                                                                                                                            cTGCGGTGAGGGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGtgctgc     <  1:1117046/238‑1 (MQ=255)
                                                                                                                                                                                                                                               cGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc  >  2:678394/1‑238 (MQ=255)
                                                                                                                                                                                                                                               cGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATc  >  1:467246/1‑238 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                            
GTTCCCATCGGCGTGAGGGAAAGCCGACGCGCCAGATCGTAGTTGCTGGTGCCCACAAACCAGGATTCCTGTTGCAGACGCTTAACGATGGTTTCTTGTACTTCGCGAGAAAAACGGCGACGGGTGCCAAAATCCATCAGATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAACCAGTTAAGATAGTCGGCCTTAAAGAAAGGCAGGGCAGAAAGCCACTGATATTCATCATCCTGCAGGCGCAGGTGCTGCATC  >  NC_000913/989707‑990177

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