Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,388,502 G→C 7.5% A290G (GCC→GGC)  mscM ← mechanosensitive channel protein, miniconductance

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,388,5020GC 7.5% 408.0 / 20.8 120A290G (GCC→GGC) mscMmechanosensitive channel protein, miniconductance
Reads supporting (aligned to +/- strand):  ref base G (56/54);  new base C (4/5);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 7.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-02

Only 100 of 131 total aligned reads displayed.

AGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGT  >  NC_000913/4388235‑4388775
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                  
aGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTCCCCCCGCAGCGCTTCGCCGACCAGCTTGGCAGATCCCCACCATTGCGACTGTTCACGCACCGTATTCCAAGCCCGCCGGACCGGTAACGTCCGGCCGTCACCCt                                                                                                                                                                                                                                                                              >  4:482386/1‑274 (MQ=255)
  attACGGTTTTGCTCGGCAGTCAGCGGCTGACCGCCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCTTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGa                                                                                                                                                                                                                                                                            <  2:378549/272‑1 (MQ=255)
        gTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTggt                                                                                                                                                                                                                                                                     >  2:153706/1‑273 (MQ=255)
            tGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCCGCCTGACGCTGTTGCGa                                                                                                                                                                                                                                                                 >  2:291213/1‑275 (MQ=255)
             gCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAg                                                                                                                                                                                                                                                                <  2:181968/275‑1 (MQ=255)
                 gggattcaggtgCTGACGGTCGGCCTGAGGAGTTTTCGGGAGGAGCGGCCTTTTATTAAGACAAACCTGATAACGTACCCGGGGAACAGGCAACTGGGCCATGGCGGAATCATGCGGGGGTGGTGTCGGAATGTCCGGCCGCCGTGCCACCTGCGCCCGAAGCGCTTCGCCGAGCAGATTGGACGATCCCGGAAATTGCGACTGTTCACGCAGCGTATTCAACGCCGGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGGGAGGGaa                                                                                                                                                                                                                                                            <  4:242518/264‑1 (MQ=255)
                                         caggtgtaTTTGACGGAGGAGGGGCTGTTGTTTAAGCATATCCTCATGATGTAACGGTGTCACACGCAAGTGGGCGATTGCGTGCTCAAGCTTTTGTGGTTGCGGCATTGCCGGCAGCCGTGCAACCTGCGCCCGCAGCGCGTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGt                                                                                                                                                                                                                                     <  4:290354/267‑1 (MQ=255)
                                         cTGATGAATTTGCCGTAGCATCGGCTGTTTATTAAGCAGCTCCTCATAACGTAACCGTGGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCCGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTGCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGt                                                                                                                                                                                                                                     <  1:39906/274‑1 (MQ=255)
                                            atgaATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAt                                                                                                                                                                                                                                 >  2:191016/1‑275 (MQ=255)
                                                             agagacagtttttaAACAGATCCCTATAACGTACCCGTTGCACACGGAACTGCGCCATTTCGGTATCAGGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGGGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                                                                                                                                                                <  4:142713/264‑1 (MQ=255)
                                                                                aTCCTCATAACGTAATCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTGGATTCAAAGCCGCCGAACGTTCGCGGTTTACTTTGAAtt                                                                                                                                                                                             >  2:92758/1‑275 (MQ=255)
                                                                                aTCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGa                                                                                                                                                                                                                    >  2:142714/1‑252 (MQ=255)
                                                                                    gacagaCGGAACCTGTGCACACGCAACTGCGCCCTTTCGTTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgcg                                                                                                                                                                                          <  1:1596/269‑1 (MQ=255)
                                                                                          cGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGATGCACCGAGATCCATCCGCTGCGCCTGTTCATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGCATTGCGcgacga                                                                                                                                                                                    >  1:98243/1‑274 (MQ=255)
                                                                                          cGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgcc                                                                                                                                                                                                                         >  2:150254/1‑237 (MQ=255)
                                                                                            tAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTACTTTTGAATTGCGCGACGatat                                                                                                                                                                                 >  2:258957/1‑275 (MQ=255)
                                                                                              aCCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGccgccg                                                                                                                                                                                                                     <  4:150253/237‑1 (MQ=255)
                                                                                              aCCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGCTGTGGTTTCGGCATTTCAGGCAGCCGTGCCACCCGTCCCCGCAGCGCTTCGCCGAGCAGATTGCACGATCCCAGCCATTGAGACTGTTCACGCAGCGTACTCAACGCCTGCCGGACCTGTAACGTCCGGCTGGCCACCTGAGTCGGCTGCTAGGCACCGGGACCCCCCCGCCGCCGCCGGTGGTTCAAGACCGCCACTAGTTCCGGGTACATGTTGGATTgcccacccataga                                                                                                                                                                                >  4:27855/1‑263 (MQ=255)
                                                                                                     aactcggaaCTGTGCCATGTCGGTATCAAGCTGTTGTGGTGTCGGCAGTTCCGGCAGCCTTGCCACCTGAGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGTGACTGTTCACGCAGCGTAGTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCa                                                                                                                                                                        <  2:487194/268‑1 (MQ=255)
                                                                                                                      tttCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGCCGGTATCTTTCGGCCAATCGGCGCTGGTTTCt                                                                                                                                                      >  2:148626/1‑275 (MQ=255)
                                                                                                                           gtttCAAGCTGTTGTGGTTTCGGCATGTCCGGCAGCCGGGCCACCTGGGCCCCCCGGGCTTCGCCGGCCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAg                                                                                                                                                 <  4:103914/272‑1 (MQ=255)
                                                                                                                           gTATCAAGCTGGTATGGTTTCGCCACTTCCGGCAGCCGTGCCCACTGCGCCCGGAGCGCTACGCCGAGCAGTGTGGACGATCCCAGCCATTGCGACTGTTCACGCCGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAg                                                                                                                                                 <  4:367104/275‑1 (MQ=255)
                                                                                                                               cAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAGCGTCAGGCTGCCAGccag                                                                                                                                                                                                                                                                             >  5:578788/1‑146 (MQ=255)
                                                                                                                                aaGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGa                                                                                                                                                                                                                                                                            <  6:502440/148‑1 (MQ=255)
                                                                                                                                 aGCTGTTGTGGTTTCGGCATTTCTGGCAGCCGTGCCACCTGCGCCCGCCGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCCGTAACGTCTGGCTGNCAGCCTTCCGCTGTTGCCAGGCCACGAGATCCATCCGCTGCGCCTCTTGTTTCCAAGCCCCCGCTATTTCGCGGTTAATTNGGAATTTCGCGCCGACACCTTTCGGNAAAGACGCGCTCTTTTCTTccccccatt                                                                                                                                            >  3:91835/1‑267 (MQ=255)
                                                                                                                                     gttgtGGTTTCGGCGTTTCCGGCCGCCGTGCCACCTGCGCCCGTCGCGTTTCGCCGGGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCGGCGGATTCAACGCCTGCCGGACCTGTAACGCCTGGCTGGCAGCCTGGCGCTGTTGCGAGGCAACGGGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTa                                                                                                                                       <  4:191015/275‑1 (MQ=255)
                                                                                                                                     gttgtGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAGCGTCAGGCTGCCAGCCAGACgtta                                                                                                                                                                                                                                                                       <  6:578788/148‑4 (MQ=255)
                                                                                                                                          ggTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGt                                                                                                                                                                                                                                                                      >  5:365055/1‑144 (MQ=255)
                                                                                                                                          ggTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGt                                                                                                                                                                                                                                                                      <  6:365055/144‑1 (MQ=255)
                                                                                                                                             ttCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCt                                                                                                                                                             >  4:41005/1‑245 (MQ=255)
                                                                                                                                              tCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTTATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAg                                                                                                                               >  4:484855/1‑274 (MQ=255)
                                                                                                                                                ggCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCTGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCCACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATCGTTCGCGGTTAATTTTGAATTTCGCGCCGCTATCTTTCGGCACATCGGCGCTGTTTTCTGCCCGCCATTCTGTACTTTCCCgc                                                                                                                              >  1:124884/1‑273 (MQ=255)
                                                                                                                                                 gCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGttt                                                                                                                                                         <  2:41006/245‑1 (MQ=255)
                                                                                                                                                          gCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGccag                                                                                                                                                                                                                                                                             <  5:368413/121‑3 (MQ=255)
                                                                                                                                                                 tGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTgc                                                                                                                                                                                                                                           >  6:328914/1‑148 (MQ=255)
                                                                                                                                                                    cACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGcc                                                                                                                                                                                                                                        <  7:144739/148‑1 (MQ=255)
                                                                                                                                                                        tGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                                                                                                                                                                                                                  <  7:8744/134‑1 (MQ=255)
                                                                                                                                                                        tGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                                                                                                                                                                                                                  >  8:8744/1‑134 (MQ=255)
                                                                                                                                                                          cgcCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCAATTGCGACTGTTCACGCAGCGTATTCAAGGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTgg                                                                                                   <  1:75324/274‑1 (MQ=255)
                                                                                                                                                                              cGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGa                                                                                                                                                                                                                                                                            <  7:213776/102‑1 (MQ=255)
                                                                                                                                                                              cGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGa                                                                                                                                                                                                                                                                            >  8:213776/1‑102 (MQ=255)
                                                                                                                                                                                  gcgcTTCCGCGGGCAGATGGGGCAGTCCCAGCGACGGCGACGGTTCACGAAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCGGACGCTGTTGCGAGGCAACGCGATCCATCCGCGGCGCCGGTTGATTCAAAGGCGCCGATAGTTCGCGGTTAATTTTGAATTGCGAGACGATATCTTTCGGCAAATCGGCGCTGTTTTCGGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAAt                                                                                          <  4:302952/275‑1 (MQ=255)
                                                                                                                                                                                           cGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAg                                                                                                                                                                                                                 >  6:729300/1‑148 (MQ=255)
                                                                                                                                                                                             aGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTg                                                                                                                                                                                                                                                                             >  6:1013556/1‑86 (MQ=255)
                                                                                                                                                                                             aGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTg                                                                                                                                                                                                                                                                             <  5:1013556/86‑1 (MQ=255)
                                                                                                                                                                                             aGCAGATTGGACGATCCCAGCCATTGCGACTGGTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGatccat                                                                                                                                                                                                                                                  >  5:320368/1‑113 (MQ=255)
                                                                                                                                                                                                   ttGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGt                                                                                                                                                                                                                                                                      >  6:8815/1‑87 (MQ=255)
                                                                                                                                                                                                   ttGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGt                                                                                                                                                                                                         >  8:773806/1‑148 (MQ=255)
                                                                                                                                                                                                     ggACGATCCTAGCAATGGGGCCGGTTACCGCGGCGGATTAAACGCCTGCCGGACCTGTAACGTCGGGCTTGCAGCTTGACGCGGTTGCGAGGAAACGGGAGCCATCCGCTGCGCCGGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCa                                                                       <  4:363501/275‑1 (MQ=255)
                                                                                                                                                                                                          aTCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCaa                                                                                                                                                                                                                                                            >  6:1075000/1‑90 (MQ=255)
                                                                                                                                                                                                          aTCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCaa                                                                                                                                                                                                                                                            <  5:1075000/90‑1 (MQ=255)
                                                                                                                                                                                                          aTCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAg                                                                                                                                                                                                                           >  5:822362/1‑123 (MQ=255)
                                                                                                                                                                                                            cccAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCaa                                                                                                                                                                                                                             >  6:446196/1‑119 (MQ=255)
                                                                                                                                                                                                              cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGccaga                                                                                                                                                                                                                                                                            >  7:192004/1‑67 (MQ=255)
                                                                                                                                                                                                              cAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCAGGCTGCCAGCCAGa                                                                                                                                                                                                                                                                            <  8:192004/70‑1 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                                                                                                                                                                       <  7:489190/90‑1 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                                                                                                                                                                       >  6:53892/1‑90 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                                                                                                                                                                       >  7:795593/1‑90 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                                                                                                                                                                       <  5:61984/90‑1 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                                                                                                                                                                       <  5:53892/90‑1 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                                                                                                                                                                       <  8:795593/90‑1 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACgaga                                                                                                                                                                                                                                                       >  8:489190/1‑90 (MQ=255)
                                                                                                                                                                                                               aGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAAtt                                                                                                                                                                                             >  6:254114/1‑148 (MQ=255)
                                                                                                                                                                                                                 ccATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGc                                                                                                                                                                                                                                              <  7:89816/97‑1 (MQ=255)
                                                                                                                                                                                                                     tGCGACTGTTCTCGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGAAAAGAGCGCCCTCCGCTGCGCCTGTTGATTCAAAGCCGCCGCTAGTTCGCGGTTAATTTTGCAGTGCGCGACGATAGCTTACGGCAAAGCGGACGGTTTTTCTTGCAGCAATACGGTACTTTCCAGCGCCCGCTCAGACCCAAGCTGACGTTGTCTGTTTAATTGATTGCGTATGGCCGGCAAATAAGCATCCCGtt                                                        >  4:326227/1‑274 (MQ=255)
                                                                                                                                                                                                                     tGCGACTGTTCACGCAGCGTCTTCAACGCCTGCCGGACCTCTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCCACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCCACGATAACTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCCCCGCCTCAAGCTGACGTTGGCTGTTTAATTGTTTACGAAAGGCCTGCAAATACGCATCCAGtt                                                        >  4:466093/1‑274 (MQ=255)
                                                                                                                                                                                                                     tGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAAt                                                                                          >  2:248948/1‑240 (MQ=255)
                                                                                                                                                                                                                     tGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTg                                                       >  2:362347/1‑275 (MQ=255)
                                                                                                                                                                                                                             ttCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTggatg                                                >  2:485141/1‑271 (MQ=255)
                                                                                                                                                                                                                              tCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGccgccg                                                                                                                                                                                                                     >  5:75184/1‑109 (MQ=255)
                                                                                                                                                                                                                              tCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGccgccg                                                                                                                                                                                                                     <  6:75184/109‑1 (MQ=255)
                                                                                                                                                                                                                                      cGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCg                                      >  2:103619/1‑275 (MQ=255)
                                                                                                                                                                                                                                       tagttCAAGGACGGCCGGACCGGTAACGTCGGGCTGGCAGCCTGACCGGGTTGGTAGGAAACGGCAGCCTTCCGCTGGGCCTGTTGATTAAAAGCCGCCGATAGGTCGAGGTTAAGTTTGAATGGCGCGAAGATAGCTTTCGGCAAATCGGCGCTGTTGTCTGCCAGCAATTCGGGACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCg                                      <  4:509478/271‑1 (MQ=255)
                                                                                                                                                                                                                                       gTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCACATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGCAGTTGGCTGTTTAATTGATTACCCAAGGCCTGCACATCCGCATCCAGTTCCTTGCCCTCTTTTTCCGc                                     >  4:218138/1‑275 (MQ=255)
                                                                                                                                                                                                                                            cAACGCCTGCCGGACCTGTAACGTCTGGTGGGCAGCCTGACGCGGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGct                                <  2:327934/275‑1 (MQ=255)
                                                                                                                                                                                                                                            cAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGct                                <  4:396941/275‑1 (MQ=255)
                                                                                                                                                                                                                                            aaccGCCGGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGATGTCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGct                                <  2:56930/272‑1 (MQ=255)
                                                                                                                                                                                                                                              aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTg                                                                                                                                                                                                                                   >  8:618612/1‑79 (MQ=255)
                                                                                                                                                                                                                                              aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTg                                                                                                                                                                                                                                   <  7:618612/79‑1 (MQ=255)
                                                                                                                                                                                                                                              aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                                                                                                                                                                >  5:1054092/1‑98 (MQ=255)
                                                                                                                                                                                                                                              aCGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGt                                                                                                                                                                                                                <  6:1054092/98‑1 (MQ=255)
                                                                                                                                                                                                                                               cGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgc                                                                                                                                                                                           <  6:241549/118‑1 (MQ=255)
                                                                                                                                                                                                                                               cGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTgc                                                                                                                                                                                           >  5:241549/1‑118 (MQ=255)
                                                                                                                                                                                                                                                  cTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGCTTACGCAAGGCCTGCAAATACGCCTCCCATTGCTGGCTCCTTTTTTccgcccgcctctaac                          >  4:515113/1‑268 (MQ=255)
                                                                                                                                                                                                                                                       ggaccgtgaaCGTCTGGGTGGCAGCCGTCCGCTGTTGCGAGGCAACGGGATCCATCCTCTGCGCCTGTTGATTCAAANCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAg                     <  3:124884/267‑1 (MQ=255)
                                                                                                                                                                                                                                                       ggACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAg                     <  1:91835/275‑1 (MQ=255)
                                                                                                                                                                                                                                                            tGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAtt                                                                                                                                                                                                                                >  6:35825/1‑68 (MQ=255)
                                                                                                                                                                                                                                                            tGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGAtt                                                                                                                                                                                                                                <  5:35825/68‑1 (MQ=255)
                                                                                                                                                                                                                                                               aaCGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                                                                                                                                                             >  7:342589/1‑116 (MQ=255)
                                                                                                                                                                                                                                                               aaCGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCttt                                                                                                                                                                             <  8:342589/116‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cGTCTGGCTGGCAGCCTGACGCTGTCGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCgg                                                                                                                                                                  <  7:692989/125‑1 (MQ=255)
                                                                                                                                                                                                                                                                   tcaggctgccAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccg                                                                                                                                                                                                                        <  5:66800/66‑1 (MQ=38)
                                                                                                                                                                                                                                                                   tcaggctgccAGCCAGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAgccg                                                                                                                                                                                                                        >  6:66800/4‑69 (MQ=38)
                                                                                                                                                                                                                                                                    ctggctggCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt                                                                                                                                             >  7:546375/1‑142 (MQ=255)
                                                                                                                                                                                                                                                                    ctggctggCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAAt                                                                                                                                             <  8:546375/142‑1 (MQ=255)
                                                                                                                                                                                                                                                                      ggctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAttt                                                                                                                                                                                                   >  5:872044/1‑87 (MQ=255)
                                                                                                                                                                                                                                                                      ggctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAAttt                                                                                                                                                                                                   <  6:872044/87‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gctggcAGCCTGACGCTGTTGCGAGGCAACGCGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCAGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATc                                                             <  2:209576/219‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTg                                                                                                                                                                                            >  6:20289/1‑93 (MQ=255)
                                                                                                                                                                                                                                                                       gctggcAGCCTGACGCTGTTGCGAGGCAACGAGATCCATACGCGGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTg                                                                                                                                                                                            <  5:20289/93‑1 (MQ=255)
                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                  
AGAATACGGTTTTGCTCGGCAGTCAGCGGCTGACCGTCGGCCTGATGAATTTGCCGTAGCAGCGGCTGTTTATTAAGCAGATCCTCATAACGTAACCGTTGCACACGCAACTGCGCCATTTCGGTATCAAGCTGTTGTGGTTTCGGCATTTCCGGCAGCCGTGCCACCTGCGCCCGCAGCGCTTCGCCGAGCAGATTGGACGATCCCAGCCATTGCGACTGTTCACGCAGCGTATTCAACGCCTGCCGGACCTGTAACGTCTGGCTGGCAGCCTGACGCTGTTGCGAGGCAACGAGATCCATCCGCTGCGCCTGTTGATTCAAAGCCGCCGATAGTTCGCGGTTAATTTTGAATTGCGCGACGATATCTTTCGGCAAATCGGCGCTGTTTTCTGCCAGCAATTCGGTACTTTCCAGCGCCCGCTCCGCCTCAAGCTGACGTTGGCTGTTTAATTGATTACGCAAGGCCTGCAAATACGCATCCAGTTGCTGGCTCTCTTTTTCCGCCAGCTCTGAGCGTAAGCGCGCTAATTCCTGGCGGT  >  NC_000913/4388235‑4388775

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