Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,174,077 G→A 7.9% intergenic (‑113/+15) mfd ← / ← ycfT transcription‑repair coupling factor/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,0770GA 7.9% 511.5 / 27.1 164intergenic (‑113/+15)mfd/ycfTtranscription‑repair coupling factor/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (84/67);  new base A (7/6);  total (91/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-02

Only 100 of 166 total aligned reads displayed.

GCAGAGCATTTTGCATATCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAA  >  NC_000913/1173805‑1174348
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
gCAGAGCATTTTGCATATCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCTGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                                                                                                                                               <  2:158549/275‑1 (MQ=255)
    agCATTTTGCATATCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCCAACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGTTATGCCTCAACATATGGGGGCATATCCCTGATAAGCAATGCTGTTCATTATCCGTAATCCACATGGTTTAGCAAATTTAATACCCCCGGCTggg                                                                                                                                                                                                                                                                           >  4:242852/1‑275 (MQ=255)
                     gTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTAATTATCCGGAATCCACATGGTTTAGAAAATTGAAACGCCCCGGCTGGGGCGAGTGCtttaggggt                                                                                                                                                                                                                                                          >  4:447720/1‑268 (MQ=255)
                                       gaccggCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCggg                                                                                                                                                                                                                                        >  2:199371/1‑275 (MQ=255)
                                                                                       cgcCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCTCCGGCTTGGGCGCTGGCTTAGCGGGGGTAACTGACGGCAGCGGGCAGTGCAGGCGGGGAGAACAACAGATCGCTGCATGCGAGCTGTGAAAg                                                                                                                                                                                        >  2:357099/1‑275 (MQ=255)
                                                                                                      ccAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCTCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAAGACAACAGATCGCTGAATGCACGCTGTGAACGTTTCTTTACCAGCa                                                                                                                                                                          >  4:466391/1‑274 (MQ=255)
                                                                                                      ccAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGGGAAGGCGGGGAGAACAACAGATCGATTACTGCACACTGTGAAAGTTTTCTTACCAGCac                                                                                                                                                                         >  4:342576/1‑274 (MQ=255)
                                                                                                         gagagcgcgGCTCACCCCGTTTGCGGTGCAGAGTATAAGTATATGGTTCAGGCATAACGGTTGTCAGTTTTTCTGTAGGATATGCCTCAACATATGGGGGAATATCTCTGATAAGCAATGCTGGTCATTATCCGTAATCGAAATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGAGGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGcc                                                                                                                                                                      <  4:11241/270‑1 (MQ=255)
                                                                                                            gaagctgctCACCGCGTTGGCCGGGCAGAGTATAACGGTATTGGTCAGGCATAACGGTTGGCAGATTCTCGGTAGGATTGGCCGCAACATAGTGGGGCATATCTTTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCg                                                                                                                                                                    <  4:154048/271‑1 (MQ=255)
                                                                                                              cGCTGCTAACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCAGTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCaa                                                                                                                                                                 <  2:242853/275‑1 (MQ=255)
                                                                                                                 tgctTACCCGCTTTGAAGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATGCTCGTTAGGATATGCCTAAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAg                                                                                                                                                              <  2:287337/275‑1 (MQ=255)
                                                                                                                                   ggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                                                                                                                                            <  8:947551/147‑1 (MQ=255)
                                                                                                                                                 atatTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGa                                                                                                                                                                                                                                                                       >  5:959662/1‑138 (MQ=255)
                                                                                                                                                 atatTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGa                                                                                                                                                                                                                                                                       <  6:959662/138‑1 (MQ=255)
                                                                                                                                                   atTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCg                                                                                                                                                                                                                                                                        >  5:108664/1‑135 (MQ=255)
                                                                                                                                                   atTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCg                                                                                                                                                                                                                                                                        <  6:108664/135‑1 (MQ=255)
                                                                                                                                                    ttGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACcagc                                                                                                                           >  2:481128/1‑275 (MQ=255)
                                                                                                                                                      gTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGc                                                                                                                                                                           <  4:299771/225‑1 (MQ=255)
                                                                                                                                                       ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAACGTTTTCTTACCAGCAAGCACGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACcagca                                                                                                                          >  2:162092/1‑273 (MQ=255)
                                                                                                                                                       ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACcagca                                                                                                                          >  1:127645/1‑273 (MQ=255)
                                                                                                                                                               tAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                                                                                                                                            >  7:573308/1‑119 (MQ=255)
                                                                                                                                                               tAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgg                                                                                                                                                                                                                                                                            <  8:573308/119‑1 (MQ=255)
                                                                                                                                                                taCGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                                                                                                                                               <  8:1000739/114‑1 (MQ=255)
                                                                                                                                                                aaCGGTTGTCTGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAAGCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                                                                                                                                               >  7:1000739/1‑115 (MQ=255)
                                                                                                                                                                 aCGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGAGGAACAAACTGACAAAGGGATAAAACAGCAGGAGTTTCAgg                                                                                                               >  1:11336/1‑273 (MQ=255)
                                                                                                                                                                      tGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGACAAAACTGACAAAGGGAGAAACCAGCAGGAGTGTCAGTTCGACt                                                                                                         >  4:138821/1‑275 (MQ=255)
                                                                                                                                                                      tGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGAAACAGTACAGATGAACAACATGACAAAGGGAAAAACCAGCAGGAGTGTCAGTTCCACt                                                                                                         >  4:326317/1‑275 (MQ=255)
                                                                                                                                                                        tCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTTTCCGTATTCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTTAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTTCAGATGACCAAACTGACAAAGGGATAAACCAGCAGGAGTTTCAGTTCGACtt                                                                                                        >  2:361135/1‑274 (MQ=255)
                                                                                                                                                                           gATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTg                                                                                                                                                                                                                                   >  8:392759/1‑148 (MQ=255)
                                                                                                                                                                            aTTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTTTCAGTTCGACTTgcgg                                                                                                    >  2:434534/1‑274 (MQ=255)
                                                                                                                                                                             ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGc                                                                                                                                                                                                                                                                               <  5:553936/102‑1 (MQ=255)
                                                                                                                                                                               ttcGTTAGGTTTTGCTTAAAAATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGAAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGGGTGTCAGTTCGACTTGCGGCGAc                                                                                                <  4:481127/274‑1 (MQ=255)
                                                                                                                                                                                          atGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaaca                                                                                                                                                                                                                    >  8:959635/1‑148 (MQ=255)
                                                                                                                                                                                             ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                                                                                                                                            <  6:562198/137‑1 (MQ=255)
                                                                                                                                                                                             ccTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                                                                                                                                            >  5:562198/1‑137 (MQ=255)
                                                                                                                                                                                                      atGGGGTCATATCATTGATAAAACATGGTGTTCAATATCCATAAACGACATGGTTTAGCAAATTGAATCGCCCCGGNGGGGGCGATGGCGTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGTGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTgcg                                                                          <  3:127645/274‑1 (MQ=255)
                                                                                                                                                                                                      atGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATTGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCCGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCCTTTgc                                                                           >  2:407824/1‑273 (MQ=255)
                                                                                                                                                                                                       tGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGaa                                                                                                                                                                                                       <  7:959635/148‑1 (MQ=255)
                                                                                                                                                                                                          gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGa                                                                                                                                                                                                                                                                       >  6:58373/1‑81 (MQ=255)
                                                                                                                                                                                                          gggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGa                                                                                                                                                                                                                                                                       <  5:58373/81‑1 (MQ=255)
                                                                                                                                                                                                           ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGAc                                                                                                                                                                                                                                                 >  8:224660/1‑102 (MQ=255)
                                                                                                                                                                                                           ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGAc                                                                                                                                                                                                                                                 <  7:224660/102‑1 (MQ=255)
                                                                                                                                                                                                           ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                           >  7:906290/1‑124 (MQ=255)
                                                                                                                                                                                                           ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCa                                                                                                                                                                                                   >  6:525466/1‑148 (MQ=255)
                                                                                                                                                                                                             cATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAgg                                                                                                                                                                                                                               >  2:312660/1‑118 (MQ=255)
                                                                                                                                                                                                             cATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGAACAGCGTGCTGCGTTCATTTGTGCAAGCaga                                                                  >  2:30539/1‑275 (MQ=255)
                                                                                                                                                                                                                 tctcTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                                                                                                                                           <  4:312659/118‑1 (MQ=255)
                                                                                                                                                                                                                  ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                                                                                                                                                     <  5:721414/75‑1 (MQ=255)
                                                                                                                                                                                                                  ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCtt                                                                                                                                                                                                                                                                 >  7:437572/1‑79 (MQ=255)
                                                                                                                                                                                                                  ctctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCtt                                                                                                                                                                                                                                                                 <  8:437572/79‑1 (MQ=255)
                                                                                                                                                                                                                    ctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                                                                                                                                                  <  7:427956/124‑1 (MQ=255)
                                                                                                                                                                                                                    ctGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                                                                                                                                                  >  8:427956/1‑124 (MQ=255)
                                                                                                                                                                                                                      gATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTg                                                         >  4:99036/1‑275 (MQ=255)
                                                                                                                                                                                                                          aGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGtttt                                                                                                                                                                                    <  7:392759/148‑1 (MQ=255)
                                                                                                                                                                                                                              aTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                                                                                                                                                       >  8:814255/1‑109 (MQ=255)
                                                                                                                                                                                                                              aTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                                                                                                                                                       <  7:814255/109‑1 (MQ=255)
                                                                                                                                                                                                                               tGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGACGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGGGAACGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGGACAAACTGACAAAGGGATAAACCAGCAGGAGTTTCCAGTCACATTGCGGCAGCCATGGTGCACAGTTAATTGGCGAAAGCAAGGTTAAGCTGATTAtttac                                                >  4:297359/1‑273 (MQ=255)
                                                                                                                                                                                                                               tGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCa                                                >  2:362024/1‑275 (MQ=255)
                                                                                                                                                                                                                                 cTGTTCATTGTCCGTAAGCGACATGGTTTAGCAAAATGCATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCTGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAATCCCGCAGCAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAAc                                              <  1:62421/275‑1 (MQ=255)
                                                                                                                                                                                                                                 cTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGa                                                                                                                                                                                                               >  2:565317/1‑114 (MQ=255)
                                                                                                                                                                                                                                    ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAg                                           >  2:249888/1‑275 (MQ=255)
                                                                                                                                                                                                                                     tCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                                                                                                                                                  <  5:735495/107‑1 (MQ=255)
                                                                                                                                                                                                                                     tCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                                                                                                                                                  >  6:735495/1‑107 (MQ=255)
                                                                                                                                                                                                                                     tCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGc                                                                                                                                                                                                           <  4:565316/114‑1 (MQ=255)
                                                                                                                                                                                                                                        ttAGCCGTAATCGACATTGTTTCACAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCa                                                                                              <  1:68721/220‑1 (MQ=255)
                                                                                                                                                                                                                                         tATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGc                                      >  2:426191/1‑275 (MQ=255)
                                                                                                                                                                                                                                              gTAATCGACAAGGGTTATCGAATGGAATGGGCCTGGATGGGGCGATGGATTAGCGGGGGTAACGGATGGCAGGGGGCAGGGAAGGCGGGGAGAACAACAGATCGCTGAATCAACGCTGTGAAAGTATTCTTAAAAGCAAGCCCGGAACAGGAAAGATGAACAAACTGACAAAGGGATAAGCCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATgg                                 <  4:184102/275‑1 (MQ=255)
                                                                                                                                                                                                                                               tAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                                                                                                                                                     <  5:391418/46‑1 (MQ=255)
                                                                                                                                                                                                                                               tAAACGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                                                                                                                                                     >  6:391418/1‑46 (MQ=255)
                                                                                                                                                                                                                                                aaTCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGggt                               >  2:138019/1‑275 (MQ=255)
                                                                                                                                                                                                                                                   cGAAATGGTTTAGCAAGTGGAGTGGCCACGGTTGGGGCGAGAGCTTACGGGGATTAACTTACGGCAGAGGGAATTGAAGGCGGGGCGAACAACAGATCGCTGAATGCACGCTGTGAAAGTGTTCTGACAAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGggtggt                            <  4:178788/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      cATGGTTTAGCAAAGTGAATCGCCCCGGCTGGGGGGATGGCTTAGGGGGAGTAACTGACGGAAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTAAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTa                         <  4:64801/275‑1 (MQ=255)
                                                                                                                                                                                                                                                      atgggTTTAACAATTGGAAGCGNNCGGGGNGGGGTGGTGGATTAGCGGGAGTAACGGACGGCAGGGGGGAGTGAAGGCGGGGAGAAAAACAGATCGCTGAAAGCACGATGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTa                         <  3:11336/272‑1 (MQ=255)
                                                                                                                                                                                                                                                         ggTTTGGCAAATTGAACCGCCCCGGCTGGGGGGATGGCTTAGCGGGGGTAACGGACGGCAGCGGGCAGGGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGGGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGa                       <  4:82003/274‑1 (MQ=255)
                                                                                                                                                                                                                                                              aGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCaa                 <  4:362023/275‑1 (MQ=255)
                                                                                                                                                                                                                                                               gCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                                                                                                                                                  >  7:34216/1‑81 (MQ=255)
                                                                                                                                                                                                                                                               gCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaacaac                                                                                                                                                                                                                  <  8:34216/81‑1 (MQ=255)
                                                                                                                                                                                                                                                                  aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGa                                                                                                                                            >  8:714944/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                                                                                                                                                       >  8:942/1‑72 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGgaga                                                                                                                                                                                                                       <  7:942/72‑1 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACa                                                                                                                                                               >  7:389108/1‑128 (MQ=255)
                                                                                                                                                                                                                                                                   aTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACa                                                                                                                                                               <  8:389108/128‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGgtggtg                           >  2:82004/1‑258 (MQ=255)
                                                                                                                                                                                                                                                                        aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTa                                                                                                                                                                                <  7:151061/106‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTa                                                                                                                                                                                >  8:151061/1‑106 (MQ=255)
                                                                                                                                                                                                                                                                        aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAcc                                                                                                                                                                              <  5:93285/108‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAcc                                                                                                                                                                              >  6:93285/1‑108 (MQ=255)
                                                                                                                                                                                                                                                                         tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAgg                                                                                                                                                                                                                               >  7:649514/1‑58 (MQ=255)
                                                                                                                                                                                                                                                                         tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAgg                                                                                                                                                                                                                               <  8:649514/58‑1 (MQ=255)
                                                                                                                                                                                                                                                                         tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg                                                                                                                                     >  7:955980/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                         tCGCCCCAGCCGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGCGAACAACAGATCGc                                                                                                                                                                                                           <  6:612273/78‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGc                                                                                                                                                                                                    <  7:1067878/82‑1 (MQ=255)
                                                                                                                                                                                                                                                                            ccccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGc                                                                                                                                                                                                    >  8:1067878/1‑82 (MQ=255)
                                                                                                                                                                                                                                                                             cccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGTATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATaa                                                                                                                                 >  8:906512/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                             cccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCt                                                                                                                                                                                  >  8:888339/1‑99 (MQ=255)
                                                                                                                                                                                                                                                                             cccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGaa                                                       >  2:521406/1‑222 (MQ=255)
                                                                                                                                                                                                                                                                              ccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTa                                                                                                                                                            >  8:51919/1‑120 (MQ=255)
                                                                                                                                                                                                                                                                              ccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTa                                                                                                                                                            <  7:51919/120‑1 (MQ=255)
                                                                                                                                                                                                                                                                               gggCTGGGGAGATGGCTGAGCGGGAGTAACGGACTGGAGCGGGCAGTGAAGGCGGGGAGAACAACAGCGCGCTGAATGCAAGCTGTGCAAGTATTCTGACCAGCAATCCCGGAAAAGTACAGAGGAGCAAAATGACAAAGGGATAAACCAGCAGGAGGGTCAGTTCGACTTGCGGCGACCAGCGGGCTGCGTTCATTGGCGGAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCaa  <  4:128644/272‑1 (MQ=255)
                                                                                                                                                                                                                                                                                ggCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGTATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAAcc                                                                                                                              <  7:906512/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                                               
GCAGAGCATTTTGCATATCTGGTGCAATGAGTACCACCGGACCGGCGTGACGTTCGGCAATTTCCGCTACCAGCGTTGCACAGGCTGCGCCGGTTAACTCGCCCAGCAGACGCTGCTCACCCGCTTTGACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAGGAGTGTCAGTTCGACTTGCGGCGACCAGCGTGCTGCGTTCATTTGCGCAAGCAGAGTTAAGCTGAATATTTCAACCAGAATGCGATGGGTGGTGTAGATAGCAATGGTGTTGGAGCCAA  >  NC_000913/1173805‑1174348

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