Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,639,338 T→C 5.2% intergenic (+197/+47) 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,3380TC 5.2% 311.7 / 3.8 94intergenic (+197/+47)pitA/uspBphosphate transporter, low‑affinity; tellurite importer/universal stress (ethanol tolerance) protein B
Reads supporting (aligned to +/- strand):  ref base T (46/43);  new base C (2/3);  total (48/46)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

GCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAA  >  NC_000913/3639067‑3639576
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
gCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTCATCGTTCGCACCCAAACGAGCGGGTCCGCTGGCCCGCTTCCGATTGTGACCTAGTGCGCTTTGCTCATGCAGGATGCGGCGTGAACTCCTTATCCGCCTTACAACAACATGCACATCCACTAGCTTGCAGAGCTTATGTCGCCCTGATATGCGTAGCATCCAACGCAATTTACCCTTTGACATAAGTCTCGACCGGGTcc                                                                                                                                                                                                                                              >  4:459646/1‑273 (MQ=255)
                 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                                                                                                                                                                                                                            <  6:473288/119‑4 (MQ=21)
                                                                                                                                                                             tACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGCCAGTAAACTGACCCGCtcga                                                                                                                                                                                                                            >  5:473288/1‑116 (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                                                                                                                                                                                                       >  5:17305/1‑98 (MQ=255)
                                                                                                                                                                                                                       gCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGa                                                                                                                                                                                                       <  6:17305/98‑1 (MQ=255)
                                                                                                                                                                                                                        ccTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGc                                                                                                                                                                                                          >  4:244873/1‑94 (MQ=255)
                                                                                                                                                                                                                          tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcatc                                                                                                                                                                                  <  8:33857/116‑1 (MQ=255)
                                                                                                                                                                                                                          tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAatcatc                                                                                                                                                                                  >  7:33857/1‑116 (MQ=255)
                                                                                                                                                                                                                          tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                                                                                                    <  8:52344/146‑1 (MQ=255)
                                                                                                                                                                                                                          tGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACca                                                                                                                                                    >  7:52344/1‑146 (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                                                                                                                                              >  5:188171/1‑101 (MQ=255)
                                                                                                                                                                                                                                                                             ggTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAAt                                                                                                                                              <  6:188171/101‑1 (MQ=255)
                                                                                                                                                                                                                                                                              gTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATaa                                                                                             >  6:175848/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                              gTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACg                                      >  4:204098/1‑204 (MQ=255)
                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                              
GCCTGGGTGTTTACCCTTCCGGCTGCGGTACTGCTTTCCGGCGGGCTGTACTGGCTCTCCTTGCAGTTCCTGTAATCGTACGCACCAAAACGAGCGGGTCAGCTGGCCCGCTTCAGATTGTGACATAGTGCGCTTTGTTCATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAAAACATGCAAATTCAATAGATTGCAGAGATTATGTAGGCCTGATAAGCGTAGCGCATCAGGCAATTTTGCGTTTGTAATCAGTCTCGAGCGGGTCAGTTTACTGGCCCGCTTTTTTTATGCGTTATTATTGCTGAAGTGCTTAATGCCAAATCATCAATGCAATCAGGCTGACCACCACCAGACCACACAATGCGCTGGTCAGAATAAACTGCCGACGCACCCGCTCACAGCGGCGAATAAATTCATCATCGTGATGATCGCGATAACGTTGGGCATAGATATACCAAACGAGACGCACCTGTTTGTTGGGTTGGCCATGTGAGGTAAAAAA  >  NC_000913/3639067‑3639576

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