Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,344 T→A 5.3% intergenic (+203/+41) pitA → / ← uspB phosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,639,3440TA 5.3% 332.6 / 7.1 96intergenic (+203/+41)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base T (44/46);  new base A (2/3);  total (47/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.02e-01

TCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAA  >  NC_000913/3639084‑3639576
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
tCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCCAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGCTGCGGCGTGCACGCCTTATCCGGCCTACAAACACATGCAAATTCAATAGATTGCAGAGATTATGTAGGACTGATAAGCGGCGCGCAACAGGGTAATTTGCGTTTGTCATTAGTCTCGAGCGGGTCGGGTTAAGGGACCGGtttg                                                                                                                                                                                                                            >  4:406557/1‑274 (MQ=255)
            aCTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttat                                                                                                                                                                                                                <  2:417178/275‑1 (MQ=255)
               gCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGAATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTtattat                                                                                                                                                                                                             <  2:309279/275‑1 (MQ=255)
                          ggCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCACAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTCACATAGTCCGCTTTGTTCATGCCGGATGCGGCGTGAACGCATTCTCCGGCCTACAAAAACATGCAAATTCATTCGATGGCAGAGATTATTTAGGCCTGATAAGCGTAGCGCACCAGGCCATTTTGCGTTTGTAATCCGTCTCGATCGGGTTCGTTTTGTGGCCAGCTTTTTTTATGCTTTTTTTTTgtgtaaattt                                                                                                                                                                                                  >  4:498287/1‑266 (MQ=37)
                                        tcCTTGCAGTTCCTGTAATCATACGCACCAAGACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTGTCCGGCATACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTa                                                                                                                                                                                                                                         <  2:203764/222‑1 (MQ=255)
                                         ccTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatc                                                                                                                                                                                     >  2:470726/1‑273 (MQ=255)
                                            tGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCaa                                                                                                                                                                                <  2:51466/275‑1 (MQ=255)
                                              cAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGAGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGAACATGCAAATTCAATAGATTGCAGAGCTTCTGTAGGCCTGATAAGCGTACCGAATCAGGCAATTTTGCGTTTGTAATCAGTTCCGAACGGGTCAGTTTCCTGGCCCGCTTTTTTTCTCCTTTATTATTTCTTACATGCTTAATTCCAAatctctttc                                                                                                                                                                               >  4:2550/1‑268 (MQ=25)
                                                  tCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAGAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTa                                                                                                                                                                                               <  2:2550/254‑1 (MQ=255)
                                                    cTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCa                                                                                                                                                                        >  2:234274/1‑275 (MQ=255)
                                                      gTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAgg                                                                                                                                                                      >  2:215098/1‑275 (MQ=255)
                                                      gTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGGTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTAAGCGTTATTATTGCTGAAGTTCTTAATGCCAAATCAACAATGCAATCAgg                                                                                                                                                                      >  2:236444/1‑275 (MQ=255)
                                                       aaaTCTTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGGGAACGCCTTAGCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGGAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGc                                                                                                                                                                     <  2:331440/274‑1 (MQ=255)
                                                               cGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAATGCCTTATCCGGCCTACAAAAATATTCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGCCTTGAGCGAGTCAGTACACTGGCCCGCTTTTTTTATGAGTTATTATTGCTGAAGTGCTTAATGACAAATCATCAATACAATCAGGCTGaccccc                                                                                                                                                             >  2:308748/1‑275 (MQ=37)
                                                                gCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccacca                                                                                                                                                            <  2:215337/275‑1 (MQ=255)
                                                                                 tCAGCTGGCCTGCTTCAGATGGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCAGCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCCTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgcg                                                                                                                                           <  2:28126/275‑1 (MQ=255)
                                                                                 tCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgcg                                                                                                                                           <  2:406562/275‑1 (MQ=255)
                                                                                             cTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGAAGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCCGGCAATTTTGAGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACGGGCCCGCTTTTTTTCTGCGTTATTATTGCTGAAGTGCTTAATGCCAAATAATAAATGCAATCAGGCTACCCACCACCAGACCACCCAATCCGCTGGTCAGAAt                                                                                                                                >  4:334353/1‑274 (MQ=37)
                                                                                                          cATAGGGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCa                                                                                                                  <  2:334357/275‑1 (MQ=255)
                                                                                                                gcgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGttt                                                                                                                                                                                                                                          >  8:405987/1‑149 (MQ=255)
                                                                                                                 cgcTTTGTTCATGCCTGATGCGGCGTGAACGCCTTATCCGGCCTACACAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCCGAATAAACTGCCGACGCCCCCGCTc                                                                                                           >  2:519299/1‑275 (MQ=255)
                                                                                                                 cgcTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTc                                                                                                           >  2:487985/1‑275 (MQ=255)
                                                                                                                    tttGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTg                                                                                                                                                                                                                                      >  7:137989/1‑149 (MQ=255)
                                                                                                                          cATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAACACATTGCAGAGAGTATGTAGGCCTGATAAGCGTAGCGCATCAGGGAACTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGc                                                                                                                                                                                                                                <  6:186172/149‑1 (MQ=18)
                                                                                                                                ggATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCtttttt                                                                                                                                                                                                                          <  7:405987/149‑1 (MQ=255)
                                                                                                                                  aTGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCCATGCAATCAGGCTGACCACCACCAGACCACACAATTCGCTGGTCAGAATAAACTGCCGACGCACCCGCCCACAGCGGCGAAGAAAtt                                                                                          >  2:35079/1‑275 (MQ=255)
                                                                                                                                  aTGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCACTGCAATCAGTCTGACCACCACCAGACCACACAATGCGCTGTCCAGAATAAACTGCCGCCCCACCCGCTCACAGCGGCACATAAAtt                                                                                          >  2:122851/1‑275 (MQ=255)
                                                                                                                                   tGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTAt                                                                                                                                                                                                                       <  7:230317/149‑1 (MQ=255)
                                                                                                                                     cggcgTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGc                                                                                                                                                                                                                     <  7:16752/149‑1 (MQ=255)
                                                                                                                                            aagctctaaTCCGTACTAAAAAAAAAGGGAAATTCAATAGTTTTGAGGGATTATGTAGGCCTGATAAGGGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGtg                                                                                <  4:487980/267‑1 (MQ=255)
                                                                                                                                            aaCGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttatta                                                                                                                                                                                                              <  8:315290/149‑1 (MQ=255)
                                                                                                                                              cGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTattatt                                                                                                                                                                                                            <  8:428683/149‑1 (MQ=255)
                                                                                                                                                    aTCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCACTCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACCCACCCGACCACAGCGGCGAATAAATTCATACACGTGATGATcca                                                                        >  4:215583/1‑273 (MQ=255)
                                                                                                                                                    aTCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGAGAAAAACTTCATCATCGAgatgataga                                                                        >  1:13819/1‑272 (MQ=255)
                                                                                                                                                       cGGCCGAAAAAAACATGCAAATTCAATGGATGTCAGGGGTTATGTGGACCTGATAAGATGGGCGCATCAGTCAACTTTGGTTTTGTAATCGGTATCGAGGGGGTCAGTTTACTGGCCCGCGTTTTTGATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCCGACCACACAATGCGCTGGTCAGAATAAACTGCCTACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGa                                                                      <  4:215095/274‑1 (MQ=255)
                                                                                                                                                        ggCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGc                                                                                                                                                                                                  >  7:307219/1‑149 (MQ=255)
                                                                                                                                                           atataaaaaCAGGCAAAGTCAAAAGAGTTCAGAGTGTATGTAGGCCTGATGAGCGTAGCGCATCAGGGAATTTTGCGTTTGCAATCAGTCTCGAGAGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCGGAATTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGGTGATCGCGATAACg                                                                 <  4:122850/271‑1 (MQ=255)
                                                                                                                                                            tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTaa                                                                                                                                                                                              <  8:137989/149‑1 (MQ=255)
                                                                                                                                                            tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtcga                                                                                                                                                                                                                            >  5:473288/1‑116 (MQ=21)
                                                                                                                                                            tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtcga                                                                                                                                                                                                                            <  6:473288/119‑4 (MQ=21)
                                                                                                                                                                aaaaCATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGATCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTggg                                                            <  2:498292/275‑1 (MQ=255)
                                                                                                                                                                      tGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                                                                                                    >  4:33902/1‑181 (MQ=255)
                                                                                                                                                                       gCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcat                                                                                                                                                                                   >  6:315184/1‑149 (MQ=255)
                                                                                                                                                                        cAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCCCAGCGGCGAATAAATTCATCATCGTGATGATc                                                                          >  1:18840/1‑253 (MQ=255)
                                                                                                                                                                          aaTTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACcacaca                                                                                                                                                <  2:33902/181‑1 (MQ=255)
                                                                                                                                                                             tCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGc                                                                                                                                                                             <  5:315184/149‑1 (MQ=255)
                                                                                                                                                                               aaTAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcat                                                                                                                                                                                   <  5:234701/141‑1 (MQ=255)
                                                                                                                                                                               aaTAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcat                                                                                                                                                                                   >  6:234701/1‑141 (MQ=255)
                                                                                                                                                                                   gATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAg                                                                                                                                                                       <  6:190452/149‑1 (MQ=255)
                                                                                                                                                                                    aTTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACgaga                                        <  2:50800/275‑1 (MQ=255)
                                                                                                                                                                                      tGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGacc                                                                                                                                                                >  2:538819/1‑153 (MQ=255)
                                                                                                                                                                                       gCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACg                                      <  2:431917/274‑1 (MQ=255)
                                                                                                                                                                                         agagATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTCATGCCAAATCATCAATGCAATCAGGCTGac                                                                                                                                                                 >  8:219604/1‑149 (MQ=255)
                                                                                                                                                                                          gagaTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccacca                                                                                                                                                            <  4:538814/153‑1 (MQ=255)
                                                                                                                                                                                          gagaTTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCAcc                                  <  2:459651/275‑1 (MQ=255)
                                                                                                                                                                                              ttATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGTCCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCCCACAGCGGCGAATAAATTCATCATTGTGATGCTcg                                                                         >  4:298044/1‑232 (MQ=255)
                                                                                                                                                                                                ttGTAGGCCTGGTAAGCGTAGCGAATCAGTCAATTTTGCGTTTGTAATCAGTCTCGTGCGGGTCAGTTTCCTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGAAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCtgtttg                            <  4:234270/274‑1 (MQ=255)
                                                                                                                                                                                                  gTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGAt                                                                     <  2:298048/232‑1 (MQ=255)
                                                                                                                                                                                                      gCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGa                                                                                                                                                                                                       <  6:17305/98‑1 (MQ=255)
                                                                                                                                                                                                      gCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGa                                                                                                                                                                                                       >  5:17305/1‑98 (MQ=255)
                                                                                                                                                                                                       ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGc                                                                                                                                                                                                          >  4:244873/1‑94 (MQ=255)
                                                                                                                                                                                                         tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcatc                                                                                                                                                                                  >  7:33857/1‑116 (MQ=255)
                                                                                                                                                                                                         tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcatc                                                                                                                                                                                  <  8:33857/116‑1 (MQ=255)
                                                                                                                                                                                                         tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                                                                                                    >  7:52344/1‑146 (MQ=255)
                                                                                                                                                                                                         tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                                                                                                    <  8:52344/146‑1 (MQ=255)
                                                                                                                                                                                                         gtaanngCGTAGCGCATCAGGCAAGTTTTCGTTTGTATTCGGTCTCAAGCGGGTNAGTTTACTGGCCCGCTATTTTAATGCGTGATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCTGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCAGCTCACAGCGGCGAATAAATtcat                                                                                       <  4:509048/201‑1 (MQ=255)
                                                                                                                                                                                                           aTAAGCGTGGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGaa                                                                                                                                                                                                      <  2:244877/94‑1 (MQ=255)
                                                                                                                                                                                                           aTAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCtttt                                                                                                                                                                                                                            >  4:486948/1‑72 (MQ=255)
                                                                                                                                                                                                           aTAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAAAAGCACCTGTTTGTTGGGTTGGcc                 >  2:351557/1‑275 (MQ=255)
                                                                                                                                                                                                               gCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTa                                                                                                                                                                                                                        <  2:486953/72‑1 (MQ=255)
                                                                                                                                                                                                                cGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACaa                                                                                                                                               >  1:81789/1‑144 (MQ=255)
                                                                                                                                                                                                                   aGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAAAGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTgg                                                                                                                                       >  6:53506/1‑149 (MQ=255)
                                                                                                                                                                                                                    gcgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATgcg                                                                                                                                           <  3:81789/144‑1 (MQ=255)
                                                                                                                                                                                                                     cgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCg                                                                                                                                                                                                                    >  4:249114/1‑70 (MQ=255)
                                                                                                                                                                                                                     cgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGaccacca                                                                                                                                                            >  1:84341/1‑126 (MQ=255)
                                                                                                                                                                                                                     cgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGACTAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGAACCGGTTTGTTGGGTTGGCCATGTGAGGTa       >  1:151509/1‑275 (MQ=255)
                                                                                                                                                                                                                     cgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATa                                                                    >  2:190333/1‑214 (MQ=255)
                                                                                                                                                                                                                     cgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGt        >  4:183844/1‑274 (MQ=255)
                                                                                                                                                                                                                     cgcATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACACGCTCACAGCGGCGAATAAAt                                                                                           >  2:509053/1‑191 (MQ=255)
                                                                                                                                                                                                                        tgtgggCAATTTTTCGTTTGTAATTGTTGTCGAGAGGGTAAGTTTACTGGCCGGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTaaaa    <  4:351552/271‑1 (MQ=255)
                                                                                                                                                                                                                         tCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGttat                                                                                                                                                                                                                <  2:249118/70‑1 (MQ=255)
                                                                                                                                                                                                                         tCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAg                                                                                                                                                        <  3:84341/126‑1 (MQ=255)
                                                                                                                                                                                                                          cAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTaaaaaa  >  2:238198/1‑275 (MQ=255)
                                                                                                                                                                                                                                     gCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGAc                                                                                                                     <  7:37313/149‑1 (MQ=255)
                                                                                                                                                                                                                                        tttGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCa                                                                                                                  >  5:9254/1‑149 (MQ=255)
                                                                                                                                                                                                                                        tttGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGC‑‑TTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCAcc                                                                                                                >  7:344667/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     aaaaGCGGGCCAGTAAACTGACCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAAt                                                                                                                                                                               <  6:414122/72‑1 (MQ=25)
                                                                                                                                                                                                                                                     aaaaGCGGGCCAGTAAACTGACCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAAt                                                                                                                                                                               >  5:414122/4‑75 (MQ=37)
                                                                                                                                                                                                                                                       gAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAAAGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAgcggc                                                                                                   <  5:53506/149‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAAt                                                                                                                                              <  6:188171/101‑1 (MQ=255)
                                                                                                                                                                                                                                                            ggTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAAt                                                                                                                                              >  5:188171/1‑101 (MQ=255)
                                                                                                                                                                                                                                                             gTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATaa                                                                                             >  6:175848/1‑149 (MQ=255)
                                                                                                                                                                                                                                                             gTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACg                                      >  4:204098/1‑204 (MQ=255)
                                                                                                                                                                                                                                                               cAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAAt                                                                                           >  5:8783/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                 gTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCAcc                                  <  2:307882/204‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCAcc                                  <  2:204101/204‑1 (MQ=255)
                                                                                                                                                                                                                                                                 gTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACAACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCAGCGTGATGATCGCGATAACg                                                                 <  2:75719/173‑1 (MQ=255)
                                                                                                                                                                                                                                                                  tttACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATtca                                                                                        >  8:463057/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                        
TCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAA  >  NC_000913/3639084‑3639576

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