Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,479,020 A→T 5.3% intergenic (+131/+14) rraB → / ← yjgM protein inhibitor of RNase E/GNAT family putative N‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,479,0200AT 5.3% 461.5 / 12.0 132intergenic (+131/+14)rraB/yjgMprotein inhibitor of RNase E/GNAT family putative N‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base A (67/58);  new base T (5/2);  total (72/60)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.75e-03

Only 100 of 137 total aligned reads displayed.

CCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCC  >  NC_000913/4478747‑4479293
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
        gAACAGCTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTGTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCa                                                                                                                                                                                                                                                                            <  2:677036/274‑1 (MQ=255)
               tgatgaCGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGACCACGCGGGCATTTTTTTa                                                                                                                                                                                                                                                                    >  2:274777/1‑275 (MQ=255)
                        tGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGACGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                                                                                                                                                                                                                                                            >  2:146816/1‑274 (MQ=255)
                        tGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGACGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTcgc                                                                                                                                                                                                                                                            >  2:138386/1‑274 (MQ=255)
                                  aattGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                                                                                                                                                                                                                                                  <  2:107797/272‑1 (MQ=255)
                                            gAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                                                                                                                                                                                                        <  4:146815/274‑1 (MQ=255)
                                                acgacgGATGGGGCACTTACTTTGAAGATCCTAAAGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGAc                                                                                                                                                                                                                                   <  2:64756/275‑1 (MQ=255)
                                                       aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTATTGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                                                                                                                                                                                                      >  4:33155/1‑249 (MQ=255)
                                                           ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCg                                                                                                                                                                                                                                                  <  2:33156/249‑1 (MQ=255)
                                                               cttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAACACGTACGTTTTGAATGTTCGAATGCCCGCGAGTTGAACACACGGGCATTTTTTTAGAGTTCGCCCAACACCCGCACTTCCCAAACCAAATGGccacgtcagcc                                                                                                                                                                                                                    >  4:424531/1‑266 (MQ=255)
                                                                      tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCTCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCa                                                                                                                                                                                                             >  2:540850/1‑275 (MQ=255)
                                                                                       aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCCCTTCACAATCGACATGGCCCGTGCAGCCAAGCACATATTCGATATTCTCaaa                                                                                                                                                                                            >  1:70830/1‑275 (MQ=255)
                                                                                                                    gAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                                                                                                                                                                                                                >  2:300084/1‑162 (MQ=255)
                                                                                                                             gACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGcc                                                                                                                                                                                                                             >  2:219436/1‑204 (MQ=255)
                                                                                                                                ggAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcg                                                                                                                                                                                                                                                                                  <  5:717093/148‑1 (MQ=255)
                                                                                                                                 gAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTg                                                                                                                                                                                                                         <  4:219435/204‑1 (MQ=255)
                                                                                                                                  aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                                                                                                                                                                                                                <  5:241384/148‑1 (MQ=255)
                                                                                                                                  aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGACCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTCAAAAAAGCg                                                                                                                                                >  2:27719/1‑275 (MQ=255)
                                                                                                                                      cGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCCGCCAAGTGCATAGTCGAAATGCTCAAAGCCCAAACGCTCATAACGCGCAATGGCCTCCTTTAAAAAAGCGGtt                                                                                                                                             >  4:25262/1‑273 (MQ=255)
                                                                                                                                          aCTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                                                                                                                                                                                                        >  6:498074/1‑148 (MQ=255)
                                                                                                                                           cTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttttt                                                                                                                                                                                                                                                                      <  8:438114/149‑1 (MQ=255)
                                                                                                                                           cTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAttttt                                                                                                                                                                                                                                                                       >  6:138653/1‑148 (MQ=255)
                                                                                                                                                      tAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACAcgcg                                                                                                                                                                                                                                                                                >  2:740493/1‑128 (MQ=255)
                                                                                                                                                       aGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                                                                                                                                                                                                  >  1:32032/1‑173 (MQ=255)
                                                                                                                                                         cgcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCAACATGGCCCGTGCAGCCAAGCGCATAGTCCATATGCCCAAAGCCCAAATGCTTAATACAGGCAATGGCTTCCTTTAAAAAACGGGGTGTTTCCaagccttgtctttt                                                                                                                         >  4:568504/1‑262 (MQ=255)
                                                                                                                                                          gcATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCa                                                                                                                                                                                                                                                                            <  4:740491/128‑1 (MQ=255)
                                                                                                                                                           cATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAANTTGTAGNCCGNATAAGGCNTTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGc                                                                                                                                                                                                                              <  3:32031/173‑1 (MQ=255)
                                                                                                                                                            aTCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACAt                                                                                                                                                                                                                                                     >  7:319218/1‑149 (MQ=255)
                                                                                                                                                             tCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGt                                                                                                                                              <  2:568506/250‑1 (MQ=255)
                                                                                                                                                              cAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                                                                                                                                                                                                        >  5:912510/1‑128 (MQ=255)
                                                                                                                                                              cAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCAtttt                                                                                                                                                                                                                                                                        <  6:912510/128‑1 (MQ=255)
                                                                                                                                                                           ttCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTc                                                                                                                                                                                                                                            >  4:544644/1‑143 (MQ=255)
                                                                                                                                                                             ccTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCg                                                                                                                                                                                                                                     <  5:685273/148‑1 (MQ=255)
                                                                                                                                                                               tGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACaa                                                                                                                                                                                                                                        <  2:544646/143‑1 (MQ=255)
                                                                                                                                                                                gCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACa                                                                                                                                                                                                                                  >  5:466119/1‑148 (MQ=255)
                                                                                                                                                                                 cAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACAt                                                                                                                                                                                                                                 <  5:138653/148‑1 (MQ=255)
                                                                                                                                                                                   aCCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATgg                                                                                                                                                                                                                               >  6:813086/1‑148 (MQ=255)
                                                                                                                                                                                         tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACa                                                                                                                                                                                                                                                      >  4:814158/1‑119 (MQ=255)
                                                                                                                                                                                         tGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTc                                                                                                                                                                                                        >  1:157395/1‑165 (MQ=255)
                                                                                                                                                                                             tttGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                                                                                                                                                                                    <  3:157394/165‑1 (MQ=255)
                                                                                                                                                                                                gTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGc                                                                                  <  2:424533/275‑1 (MQ=255)
                                                                                                                                                                                                 tAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAgc                                                                                                                                                                                                                <  7:733082/149‑1 (MQ=255)
                                                                                                                                                                                                 tAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAg                                                                                                                                                                                                                 >  6:346097/1‑148 (MQ=255)
                                                                                                                                                                                                  aGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGCTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGCTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCCTAGCCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGCCTTCCTTTACAAAACCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTCAGGcc                                                                                 >  2:293854/1‑274 (MQ=255)
                                                                                                                                                                                                   gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTc                                                                                                                                                                                                                                                              >  7:30448/1‑101 (MQ=255)
                                                                                                                                                                                                   gCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTc                                                                                                                                                                                                                                                              <  8:30448/101‑1 (MQ=255)
                                                                                                                                                                                                     ccGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAACGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCACCATCGCCATTAAGGCCAGtt                                                                             >  1:89630/1‑275 (MQ=255)
                                                                                                                                                                                                           aaGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGa                                                                                                                                                                                                      <  8:319218/149‑1 (MQ=255)
                                                                                                                                                                                                             ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                                                                                                                                                                                    >  7:119634/1‑149 (MQ=255)
                                                                                                                                                                                                             ggCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGata                                                                                                                                                                                                    <  8:656406/149‑1 (MQ=255)
                                                                                                                                                                                                                  tCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCt                                                                                                                                                                                                <  5:619693/148‑1 (MQ=255)
                                                                                                                                                                                                                      gTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAg                                                                                                                                                                                           <  8:159929/149‑1 (MQ=255)
                                                                                                                                                                                                                        cacaTCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGc                                                                                                                                                                                          >  6:225045/1‑148 (MQ=255)
                                                                                                                                                                                                                                gCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                                                                                                                                                                                >  8:697159/1‑149 (MQ=255)
                                                                                                                                                                                                                                gCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCt                                                                                                                                                                                <  7:619718/149‑1 (MQ=255)
                                                                                                                                                                                                                                 cATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGc                                                 <  2:654906/275‑1 (MQ=255)
                                                                                                                                                                                                                                      ataaaaCGTACGTGGTTAAGTTTGGAAGGCCGGCGAGTTGAACAGTCGGGCATTTTTTTACGGTTCGCGCAACAACCGAACTTGACAGTAGACATGGCCCGTGCAGCCAAGCGCAGAGTCGATAGGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGgaa                                            <  4:664498/273‑1 (MQ=255)
                                                                                                                                                                                                                                      aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                                                                                                                                                                         <  5:739415/119‑1 (MQ=255)
                                                                                                                                                                                                                                      aCAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGt                                                                                                                                                                                                         >  6:739415/1‑119 (MQ=255)
                                                                                                                                                                                                                                        aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGc                                                                                                                                                                                                                        >  7:229390/1‑102 (MQ=255)
                                                                                                                                                                                                                                        aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGc                                                                                                                                                                                                                        <  8:229390/102‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa                                                                                                                                                                                             <  7:647873/129‑1 (MQ=255)
                                                                                                                                                                                                                                        aaaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaa                                                                                                                                                                                             >  8:647873/1‑129 (MQ=255)
                                                                                                                                                                                                                                         aaaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAgc                                                                                                                                                                        >  5:49407/1‑148 (MQ=255)
                                                                                                                                                                                                                                          aaCGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                                                                                                                                                  >  4:551375/1‑169 (MQ=255)
                                                                                                                                                                                                                                              tacgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGt                                                                                                                                              <  2:551377/169‑1 (MQ=255)
                                                                                                                                                                                                                                               acgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                                                                                                                                                  >  2:747751/1‑164 (MQ=255)
                                                                                                                                                                                                                                                cgtTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCa                                                                                                                                                                                                                       >  2:240267/1‑95 (MQ=255)
                                                                                                                                                                                                                                                     ttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt                                                                                                                                                           >  8:280810/1‑149 (MQ=255)
                                                                                                                                                                                                                                                     ttAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGACGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCt                                                                                                                                                           >  7:366055/1‑149 (MQ=255)
                                                                                                                                                                                                                                                          gTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTaa                                                                                                                                                       >  5:391602/1‑148 (MQ=255)
                                                                                                                                                                                                                                                           ttCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGt                                                                                                                                           <  4:100517/159‑1 (MQ=255)
                                                                                                                                                                                                                                                              gAATGCCGGCAAGTTGAACACCCGGGCATTTTTTTAGAGTTCGCGCAACAATCGCACTTCCCAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGc                    <  4:498489/275‑1 (MQ=255)
                                                                                                                                                                                                                                                              gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                                                                                                                                                  >  7:7049/1‑149 (MQ=255)
                                                                                                                                                                                                                                                              gAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                                                                                                                                                  >  8:516134/1‑149 (MQ=255)
                                                                                                                                                                                                                                                               aaTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCTATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGc                    >  2:337095/1‑274 (MQ=255)
                                                                                                                                                                                                                                                               aaTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAg                                                                                                                                                  >  5:451195/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGc                                                                                                                                                 <  5:813086/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                aTGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGGTCGCGCAACATCCGAACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATTGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCaa                  <  2:452232/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCg                                                                                                                                                <  5:498074/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                 tGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGCCGTTTCCAGATAGCAGCGTTTGAAACCCa                                                                                                                 >  2:89838/1‑179 (MQ=255)
                                                                                                                                                                                                                                                                   cccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGt                                                                                                                                              <  5:641258/148‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaaa                                                                                                                                                                                            >  7:836456/1‑102 (MQ=255)
                                                                                                                                                                                                                                                                    ccGCGTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCaaa                                                                                                                                                                                            <  8:836456/102‑1 (MQ=255)
                                                                                                                                                                                                                                                                    ccGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg                                                                                                                                            >  8:816091/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                     cgctAGTTGAACACGCGGGCATGTTTGTAGAGTTCGCGCAACATCCGCACTGGCCAAGCGACATGGCCCGTGCAGCAAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATCAAGCGCAATGGCTTCCTTTAAGAAAGCGGTCGTATCCAGATAGGAGCGTTTGAAACCCAGCTCGCGAGCCTGCTCCAGCGCCATAAAGGCCAGTTTTTTTGCCAGCCATTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGGTCGCAAATATc             <  4:274776/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgTGTTCAACTCGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCg                                                                                                                                            >  5:555849/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgGGTTGAACACGCTGGCATTTTTTTGTAGTTCGCGCAACTACCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAAGGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCAtctc                                                                                                             <  4:89837/179‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgAGTTGAACACGCGTGCATGTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCTCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATc             <  1:91724/275‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgAGTTGAACACGCGGGCCTTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTc                                                                                                                                                                               <  5:261302/113‑1 (MQ=255)
                                                                                                                                                                                                                                                                     cgcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTc                                                                                                                                                                               >  6:261302/1‑113 (MQ=255)
                                                                                                                                                                                                                                                                     cacgAGGTGAACACTCGGGCATTTTTTTANNTTTCCCGCAACCTCCGCCCTTCACAATCGANAAGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTANAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATc             <  3:89629/273‑1 (MQ=255)
                                                                                                                                                                                                                                                                      gcgAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGtt                                                                                                                                          >  8:419892/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                        gAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTc                                                                                                                                        <  8:119634/149‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATa                                                                                                                                                                            <  4:199543/112‑1 (MQ=255)
                                                                                                                                                                                                                                                                         aGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGttt                                                                                                                                         <  1:199821/147‑1 (MQ=255)
                                                                                                                                                                                                                                                                            tGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATACGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAg                                                                                                                                     >  6:591378/1‑148 (MQ=255)
                                                                                                                                                                                                                                                                              aaCACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATa                                                                                                                                  >  8:157817/1‑149 (MQ=255)
                                                                                                                                                                                                                                                                                cacGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGCTAACAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATCCATCTTTTTCAGTTCGAAAATGTCCGACCCACTcc  >  4:78536/1‑275 (MQ=255)
                                                                                                                                                                                                                                                                                 acGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAg                                                                                                                                                                         >  7:159590/1‑107 (MQ=255)
                                                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                                  
CCCAGGTTGAACAACTGATGACGCTGGCAGAGAAATTTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAACAAAACGTACGTTGTTAATGTTCGAATGCCCGCGAGTTGAACACGCGGGCATTTTTTTAGAGTTCGCGCAACATCCGCACTTCACAATCGACATGGCCCGTGCAGCCAAGCGCATAGTCGATATGCTCAAAGCCCAAATGCTCATAAAGCGCAATGGCTTCCTTTAAAAAAGCGGTCGTTTCCAGATAGCAGCGTTTGAAACCCATCTCTCGCGCCTGCTCCATCGCCATTAAGGCCAGTTTTTTTGCCAGCCCTTTGCCGCGGATAGCCGGAAGAAAATACATCTTTTGCAGTTCGCAAATATCCGACTCACTCC  >  NC_000913/4478747‑4479293

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