Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 594,950 T→G 5.1% H165P (CAT→CCT)  cusR ← response regulator in two‑component regulatory system with CusS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913594,9500TG 5.1% 617.9 / 16.2 180H165P (CAT→CCT) cusRresponse regulator in two‑component regulatory system with CusS
Reads supporting (aligned to +/- strand):  ref base T (79/89);  new base G (4/6);  total (83/97)
Fisher's exact test for biased strand distribution p-value = 7.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 182 total aligned reads displayed.

TGCAAAGGCGAAAAACGCCGCGATGGTGGCCAGGCTGATAAAAAAGGTCAGGCGGGTTGCCAGCGAAAACGGGCGCTGAAATGGCTTACTGACCATCCGGCACCTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCAT  >  NC_000913/594675‑595227
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                      
                 ccgcGATGGTGGCCAGGCTGATAAAAAAGGTCAGGCGGGTTGCCAGCGAAAACGGGCGCTGAAATGGCTTACTGACCATCCGGCACCTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCCTTGGTATCGCTGTCAATATTCATGTCCCATACCTGCGATGCCATAAACAGACGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAAcc                                                                                                                                                                                                                                                                        >  4:115896/1‑273 (MQ=255)
                      aTGGGGGCCCGGCTGATAAAAAAGGTCAGGCTGGTTACCAGCGAAAACGGGCGCTGAAATGGCTTACTGACCATCCGGCACCTCAAGCAGGTAACCCCCGCCGCGCCCGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagc                                                                                                                                                                                                                                                                   <  4:465296/274‑1 (MQ=255)
                          tggCAAGGCTGATAAAAAAGGTCAGGCCGGTTGCCAGCGAAACCGGGCGCTGAAATGTCTTACTGACCATCCGGCACCTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGt                                                                                                                                                                                                                                                              <  4:346500/275‑1 (MQ=255)
                            gCCAGGCTGATAAAAAAGGTCAGGCGGGTTGCCAGCGAAAACGGGCGCTGAAATGGCTTACTGACCATCCGGCACCTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTAGCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCGGATGGCGACGGAAAAACTCcacc                                                                                                                                                                                                                                                                   >  1:233388/1‑266 (MQ=255)
                                                   gCGGGTTGCCAGCGAAAACGGGCGCTGAAATGGCTTACTGACCATCCGGCACCTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGACGGAAGCACTCCAGCAGCGTAAACTCCTTACTTGTCCAAGTGAtc                                                                                                                                                                                                                                     >  2:639367/1‑274 (MQ=255)
                                                                                        cTGACCATCCGGCACCTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGCCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGGCGAGATCGAc                                                                                                                                                                                            >  1:109215/1‑279 (MQ=255)
                                                                                                       ctcaacactGTAACCCAGGCCCCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                                                                                                                                                                 <  2:248498/267‑1 (MQ=255)
                                                                                                       cTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGGNGCGACAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGACGGAAGAACTCCCGCAAAGTAAACTCTTTACTGGTNAAAGCTATGCGCGTGCTGCCGCTGGTGAATNNTCTGCTGACGAGAATCACCACCCAAACGGCana                                                                                                                                                                             >  3:105947/1‑277 (MQ=255)
                                                                                                                         gccgcGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACGTCCCCCTGATGGCGGAGGAAGAACCCCGGCAGAGGAAACTCTTTACCGGTCAAAGTGATGCGGGTGCCGCTGCGGGTGACTTTGCGTCTGACATGGTCGACCATCCACCCAGCCACCGCAAACTGAcat                                                                                                                                                                >  4:501611/1‑272 (MQ=255)
                                                                                                                               cACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGAt                                                                                                                                                                                                                                                                                        <  5:784822/149‑1 (MQ=255)
                                                                                                                                   gTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCg                                                                                                                                                                                                                                                                                    >  5:163031/1‑149 (MQ=255)
                                                                                                                                          ttAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCcg                                                                                                                                           <  1:105947/278‑1 (MQ=255)
                                                                                                                                              cTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAAc                                                                                                                                                                                                                                                                         >  8:241756/1‑148 (MQ=255)
                                                                                                                                                  ggCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCCAATCGGCAACCTGAAACTGACTTTCGCTAATCACCGCccccccccccc                                                                                                                                       >  4:14318/1‑269 (MQ=255)
                                                                                                                                                   gCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcag                                                                                                                                                                                                                                                                    <  7:578870/148‑1 (MQ=255)
                                                                                                                                                    cTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCcagca                                                                                                                                                                                                                                                                  <  6:544157/149‑1 (MQ=255)
                                                                                                                                                      cAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCaga                                                                                                                                                                                                                                                                <  6:244989/149‑1 (MQ=255)
                                                                                                                                                       aaaGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCagag                                                                                                                                                                                                                                                               >  6:935289/1‑149 (MQ=255)
                                                                                                                                                       aaaGTCGTTGTCAATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTa                                                                                                                                 <  2:115896/275‑1 (MQ=255)
                                                                                                                                                               tGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActc                                                                                                                                                                                                                                                        <  7:233751/148‑1 (MQ=255)
                                                                                                                                                                  cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCttt                                                                                                                                                                                                                                                     >  8:440791/1‑148 (MQ=255)
                                                                                                                                                                  cGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTa                                                                                                                                                                                                                                                    <  5:956354/149‑1 (MQ=255)
                                                                                                                                                                   gATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                                                                                                                                                                   >  6:1003129/1‑149 (MQ=255)
                                                                                                                                                                    aTTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                                                                                                                                                                   >  8:696944/1‑148 (MQ=255)
                                                                                                                                                                    aTTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACt                                                                                                                                                                                                                                                  <  6:138089/149‑1 (MQ=255)
                                                                                                                                                                     ttttGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACt                                                                                                                                                                                                                                                  >  7:828542/1‑148 (MQ=255)
                                                                                                                                                                     ttttGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTg                                                                                                                                                                                                                                                 >  5:404805/1‑149 (MQ=255)
                                                                                                                                                                          ccGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCcag                                                                                                              <  2:322240/275‑1 (MQ=255)
                                                                                                                                                                               cAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTg                                                                                                                                                                                                                                        <  7:440791/148‑1 (MQ=255)
                                                                                                                                                                                aGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGa                                                                                                                                                                                                                                       <  7:877621/148‑1 (MQ=255)
                                                                                                                                                                                  ccGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGGAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATg                                                                                                                                                                                                                                     <  7:1048033/148‑1 (MQ=255)
                                                                                                                                                                                         accaccACATAAATAGCATTGGTATCGCTGCCAAAATTCATGTCCCATACATGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGGAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGcc                                                                                                                                                                                                                             >  6:103804/1‑149 (MQ=255)
                                                                                                                                                                                         aCCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGcc                                                                                                                                                                                                                             >  6:349743/1‑149 (MQ=255)
                                                                                                                                                                                          ccGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGcc                                                                                                                                                                                                                             >  8:959153/1‑148 (MQ=255)
                                                                                                                                                                                              cacaTCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGc                                                                                                                                                                                                                        <  5:737899/149‑1 (MQ=255)
                                                                                                                                                                                                     aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGa                                                                                                                                                                                                                  >  7:722098/1‑148 (MQ=255)
                                                                                                                                                                                                     aTAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCCCTGCGGGTGa                                                                                                                                                                                                                  <  7:696944/148‑1 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAActct                                                                                                                                                                                                                                                       >  4:687299/1‑105 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGc                                                                                                                                                                                                            >  7:1053729/1‑148 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCg                                                                                                                                                                                                           >  1:121113/1‑149 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCg                                                                                                                                                                                                           <  3:121113/149‑1 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGccc                                                                                                                                             >  4:259717/1‑211 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGccc                                                                                                                                             >  2:432727/1‑211 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTa                                                                             >  2:278807/1‑275 (MQ=255)
                                                                                                                                                                                                           ttGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCCCCCCCCCCCCAGCAGTTCAGCAAAAGCGAATGGCTTCACAAGGTa                                                                             >  4:20060/1‑275 (MQ=255)
                                                                                                                                                                                                             ggTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGc                                                                                                                                                                                                         <  5:1003129/149‑1 (MQ=255)
                                                                                                                                                                                                               tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTAc                                                                                                                                                                                                                                                   <  2:687299/105‑1 (MQ=255)
                                                                                                                                                                                                               tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCcgcg                                                                                                                                         <  2:259717/211‑1 (MQ=255)
                                                                                                                                                                                                               tATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCcgcg                                                                                                                                         <  4:432727/211‑1 (MQ=255)
                                                                                                                                                                                                                aTCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTg                                                                                                                                                                                                       <  7:959153/148‑1 (MQ=255)
                                                                                                                                                                                                                  cGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCcgcg                                                                                                                                         <  4:304865/208‑1 (MQ=255)
                                                                                                                                                                                                                    cTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaagga                                                                                                                                                                                                                                                                               <  3:205810/73‑1 (MQ=255)
                                                                                                                                                                                                                    cTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCgaagga                                                                                                                                                                                                                                                                               >  1:205810/1‑73 (MQ=255)
                                                                                                                                                                                                                     tGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACgaga                                                                                                                                                                                                 <  5:37342/149‑1 (MQ=255)
                                                                                                                                                                                                                      gTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGAt                                                                                                                                                                                                <  5:349743/149‑1 (MQ=255)
                                                                                                                                                                                                                      gTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTCTGCGGCTGACGAGAt                                                                                                                                                                                                >  5:636685/1‑149 (MQ=255)
                                                                                                                                                                                                                       tCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGt                                                                                                                                                                                                                                <  5:547849/116‑1 (MQ=255)
                                                                                                                                                                                                                        cAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAGCTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATc                                                                                                                                                                                               >  8:70891/1‑148 (MQ=255)
                                                                                                                                                                                                                        cAAAATTCATGTCCCAGACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGccccc                                                                <  2:219303/275‑1 (MQ=255)
                                                                                                                                                                                                                          aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc                                                                                                                                                                                            >  5:827354/1‑149 (MQ=255)
                                                                                                                                                                                                                          aaaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc                                                                                                                                                                                            <  5:205964/149‑1 (MQ=255)
                                                                                                                                                                                                                           aaTTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGAc                                                                                                                                                                                            >  8:595270/1‑148 (MQ=255)
                                                                                                                                                                                                                              tCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAAAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTcc                                                          <  2:622534/275‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAt                                                                                                                                                                                    <  7:372146/148‑1 (MQ=255)
                                                                                                                                                                                                                                   tCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATc                                                                                                                                                                                   >  5:330277/1‑149 (MQ=255)
                                                                                                                                                                                                                                     ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAACCTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                                                                                                                                                                 <  6:390042/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                                                                                                                                                                 >  6:930929/1‑149 (MQ=255)
                                                                                                                                                                                                                                     ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                                                                                                                                                                 <  6:562178/149‑1 (MQ=255)
                                                                                                                                                                                                                                     ccATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCgg                                                                                                                                                                                 >  6:66427/1‑149 (MQ=255)
                                                                                                                                                                                                                                      cATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCaa                                                                                                                                                                                      <  7:418696/143‑1 (MQ=255)
                                                                                                                                                                                                                                        tACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCAcc                                                                                  >  1:87401/1‑241 (MQ=255)
                                                                                                                                                                                                                                         aCCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAc                                               <  2:60298/275‑1 (MQ=255)
                                                                                                                                                                                                                                          ccTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAAc                                                                                                                                                                             >  7:237844/1‑148 (MQ=255)
                                                                                                                                                                                                                                          ccTGAGGGGCGCAAAGCGGGCGGGGCAGCACATCGCCCGGGTGGCGAAGGAAGAACTCCAGACGAGTAAACCCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAGCGGCAACCTGACACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAcg                                              <  4:610460/275‑1 (MQ=255)
                                                                                                                                                                                                                                           cTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATc                                                                                                                                                                                        <  5:1016873/136‑1 (MQ=255)
                                                                                                                                                                                                                                            tggtaGGCGATAAGCGGGCGGGGAAACACTTCCACCTGATGGCGAAGGAAGAACCCCAGCATAGTAAAACTTTTAACGGTCAAAGGGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAcgcg                                            <  2:569616/271‑1 (MQ=255)
                                                                                                                                                                                                                                            tGCGAGGCTGCAAGAGAGCGGGGCACGACTTGGCCCTGCGGGGGAAGGAAGAAATCCAGCAGAGTAAAGTCTTTACTGGGCAAAGGGATGCGCGTGCCGCTGCGGGGGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTGCGATAAGCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAcgcg                                            <  4:309900/275‑1 (MQ=255)
                                                                                                                                                                                                                                            tGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAt                                                                                                                                                                                    >  8:859979/1‑139 (MQ=255)
                                                                                                                                                                                                                                            tGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAt                                                                                                                                                                                    <  7:859979/139‑1 (MQ=255)
                                                                                                                                                                                                                                            tGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGAcgcg                                            <  2:186372/275‑1 (MQ=255)
                                                                                                                                                                                                                                              cGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGa                                                                                             >  2:381985/1‑224 (MQ=255)
                                                                                                                                                                                                                                               gAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGaaa                                                                                                                                                                       <  6:404805/149‑1 (MQ=255)
                                                                                                                                                                                                                                                aGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATAAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTACTGTGCGCACCCGCGCCAGCAGTTAAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGACCCCAACTCCAGCCCCTTGACGCGATg                                         >  2:414750/1‑274 (MQ=255)
                                                                                                                                                                                                                                                  gCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATgg                                                                                         <  4:381985/224‑1 (MQ=255)
                                                                                                                                                                                                                                                           gagcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGa                                                                                                                                                            <  7:133320/148‑1 (MQ=255)
                                                                                                                                                                                                                                                            agcgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAACTCGGCAACCTGAAACTGACTATCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGAAGTTCAGCAAAAGAGAATGGCTTCACCAGGGAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAAGGGTGCCa                            >  2:464673/1‑275 (MQ=255)
                                                                                                                                                                                                                                                              cgGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAAt                                                                                                                                                        <  6:163031/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cagcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGcgcg                      <  2:54838/274‑1 (MQ=255)
                                                                                                                                                                                                                                                                   cacccattggCCCGGATGGCGAAGGAACAACGCCGGTAGAGTAAACTCTTTACGGGGAAAGGTGATCGGGGTGCAGCGGCGGGGGACTTTGCGGCTGACTAGATCGACCGTCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGGGCGCACCAGGGCCAGAAGTGCAGCAGAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGcgcg                      <  4:60079/265‑1 (MQ=255)
                                                                                                                                                                                                                                                                    agcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccg                                                                                                                                                   >  8:690280/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     ttccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgcc                                                                                                                                                 >  6:805804/4‑149 (MQ=255)
                                                                                                                                                                                                                                                                     gcaCTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCg                                                                                                                                                             >  8:87293/1‑137 (MQ=255)
                                                                                                                                                                                                                                                                     gcaCTTCGCCCTGATGGCGAAGGAAGAAACCCAGCAGGGTAAACGCTTTACCGGTCAAAGTGAGCCGCGTGCGGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAAGCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTa                   <  4:163400/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                      tccTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATa                                                                                                                                                          <  8:208753/137‑1 (MQ=255)
                                                                                                                                                                                                                                                                        cTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgc                                                                                                                                                  >  6:362953/1‑145 (MQ=255)
                                                                                                                                                                                                                                                                        cTTCGCCATCAGGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCAccgc                                                                                                                                                  <  5:362953/145‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTg                                                                                                                                                                                                                               <  7:993016/66‑1 (MQ=255)
                                                                                                                                                                                                                                                                          tCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAg                                                        <  2:329147/233‑1 (MQ=255)
                                                                                                                                                                                                                                                                          acGCCCGTATGGCGAAGGACAGAACCCGCCAGATTAAACTCTTGACTGTCCCATGTAGTGTGGCTGCCGCTGCTGGTGACTTTGCGGCTGATGAGATCGCTCGGCAGAGCGACAGCCGTAAACTGGCCATCGATAACCACCGCCGCCCCGCGCCGCAGTAATATGCGCACCCGCGCGAGCACTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCAACTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTACGCa               <  4:464673/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                                  aTGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAAACGGCAt  >  1:136755/1‑279 (MQ=255)
                                                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                                                      
TGCAAAGGCGAAAAACGCCGCGATGGTGGCCAGGCTGATAAAAAAGGTCAGGCGGGTTGCCAGCGAAAACGGGCGCTGAAATGGCTTACTGACCATCCGGCACCTCAAGCATGTAACCCACGCCGCGCACGGTCTGAATTAGCTTCGGCTCAAAGTCGTTGTCGATTTTGCCGCGCAGCCGCTTCACCGCCACATCAATAGCATTGGTATCGCTGTCAAAATTCATGTCCCATACCTGCGAGGCGATAAGCGAGCGGGGCAGCACTTCGCCCTGATGGCGAAGGAAGAACTCCAGCAGAGTAAACTCTTTACTGGTCAAAGTGATGCGCGTGCCGCTGCGGGTGACTTTGCGGCTGACGAGATCGACCATCAAATCGGCAACCTGAAACTGACTTTCGATAATCACCGCCGCCCCGCGCCGCAGTAATGTGCGCACCCGCGCCAGCAGTTCAGCAAAAGCGAATGGCTTCACCAGGTAGTCATCTGCCCCCAACTCCAGCCCCTTGACGCGATGTTCAATGGTGCCAAGCGCGGTAAGCAACAGAATCGGCAT  >  NC_000913/594675‑595227

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