Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F87 I1 R1 22 84.3 3290658 97.6% 3211682 138.4

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,964,929 C→T E234E (GAG→GAA rhlB ← ATP‑dependent RNA helicase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,964,9290CT99.5% 650.4 / NA 191E234E (GAG→GAArhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (98/92);  minor base A (0/1);  total (98/93)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 193 total aligned reads displayed.

ACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCA  >  NC_000913/3964783‑3965074
                                                                                                                                                  |                                                                                                                                                 
acacGATGACCATCTGCTGCCAGTGGTCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGGCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttc                                                                                                                                                 <  4:245933/149‑1 (MQ=255)
acacGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttc                                                                                                                                                 >  5:141908/1‑149 (MQ=255)
    gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtcttct                                                                                                                                             >  8:404214/1‑149 (MQ=255)
    gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtcttct                                                                                                                                             >  1:22211/1‑149 (MQ=255)
    gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtcttct                                                                                                                                             >  8:165801/1‑149 (MQ=255)
    gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtcttct                                                                                                                                             <  2:408556/149‑1 (MQ=255)
    gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtcttct                                                                                                                                             >  4:181970/1‑149 (MQ=255)
     gaagagaTCTGGGGCCACCTGGCAGCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                            >  7:355542/7‑149 (MQ=255)
      tGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                            >  7:411369/1‑148 (MQ=255)
        aCCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTaa                                                                                                                                         <  4:112836/149‑1 (MQ=255)
        aCCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTaa                                                                                                                                         <  5:277050/149‑1 (MQ=255)
          cATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATa                                                                                                                                       >  3:147720/1‑149 (MQ=255)
           aTCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTCGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTATCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATAc                                                                                                                                      <  6:167985/149‑1 (MQ=255)
              tgctgcCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAAt                                                                                                                                        <  7:149246/144‑1 (MQ=255)
                ctgcCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGt                                                                                                                                   <  7:404214/147‑1 (MQ=255)
                ctgcCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTg                                                                                                                                  <  3:371889/148‑1 (MQ=255)
                      ggTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGtt                                                                                                                           <  5:240564/149‑1 (MQ=255)
                       gTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGttt                                                                                                                          >  8:72832/1‑149 (MQ=255)
                           ccccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGt                                                                                                                            >  2:280687/1‑143 (MQ=255)
                           ccccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTg                                                                                                                      >  6:330941/1‑149 (MQ=255)
                           ccccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTg                                                                                                                      >  7:181729/1‑149 (MQ=255)
                            cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtctt                                                                                                                                               >  8:105062/1‑123 (MQ=255)
                            cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                            >  6:260584/1‑126 (MQ=255)
                            cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGt                                                                                                                     >  3:46478/1‑149 (MQ=255)
                            cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACCGTGGCCCGTTTTCTGt                                                                                                                     >  5:15573/1‑149 (MQ=255)
                                      ttCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTcc                                                                                                            >  1:141366/1‑148 (MQ=255)
                                        cacaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCAc                                                                                                          <  4:45735/148‑1 (MQ=255)
                                         acaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTcc                                                                                                            <  5:184749/145‑1 (MQ=255)
                                         acaACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTcc                                                                                                            >  6:184749/1‑145 (MQ=255)
                                           aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTc                                                                                                       >  7:208099/1‑148 (MQ=255)
                                           aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTc                                                                                                       >  3:305855/1‑148 (MQ=255)
                                            aCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGTCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCaa                                                                                                     <  5:202517/149‑1 (MQ=255)
                                            aCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCaa                                                                                                     <  8:220188/149‑1 (MQ=255)
                                             cGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTATTCTTCGATCAGAGTTTGCAGCAAACGCATTTTTTCTTCGATAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTACGGTTCCCCTTGaa                                                                                                     >  6:198342/1‑148 (MQ=255)
                                              ggTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGGCCGTTTTCTGTTCCGGTTCCACTTCAata                                                                                                   >  5:329941/1‑149 (MQ=255)
                                              ggTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAata                                                                                                   <  4:147720/149‑1 (MQ=255)
                                               gtgtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAata                                                                                                   >  3:17741/1‑148 (MQ=255)
                                               gtgtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                    >  7:58736/1‑147 (MQ=255)
                                               gtgtTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                    <  8:58736/147‑1 (MQ=255)
                                                  tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCAc                                                                                                          >  5:199268/1‑138 (MQ=255)
                                                  tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATAtt                                                                                               >  7:140298/1‑149 (MQ=255)
                                                  ttgggTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATAtt                                                                                               >  2:392197/4‑149 (MQ=255)
                                                      gtgtTGGCGAAAATCATCGCTCTGGCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGATCCACTTCAATATATTCgg                                                                                            <  2:133121/148‑1 (MQ=255)
                                                       tgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCa                                                                                          >  4:380124/1‑149 (MQ=255)
                                                        gtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAt                                                                                         >  2:35823/1‑149 (MQ=255)
                                                             cGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGtt                                                                                     >  1:419620/1‑148 (MQ=255)
                                                              gAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAAt                                                                                                                                        >  5:331094/1‑96 (MQ=255)
                                                               aaaaTAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCAt                                                                                  <  7:72832/149‑1 (MQ=255)
                                                                     atCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTc                                                                            <  5:311011/149‑1 (MQ=255)
                                                                     atCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTc                                                                            <  1:40993/149‑1 (MQ=255)
                                                                        gCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATAtt                                                                                               <  5:341468/127‑1 (MQ=255)
                                                                        gCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGaa                                                                         >  3:237424/1‑149 (MQ=255)
                                                                         ctctGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTATTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAAc                                                                        <  7:381310/149‑1 (MQ=255)
                                                                            tgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCa                                                                                          <  4:90531/128‑1 (MQ=255)
                                                                            tgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCa                                                                                          >  3:90531/1‑128 (MQ=255)
                                                                             gtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTc                                                                                              >  5:277299/1‑123 (MQ=255)
                                                                              tctgGCCTCTATTCTTAGATAAGAGTTTGCAGCTAAAGCATTTTTTCTTCGTTAGAAGGTTAGAAATGTTCTTCTTTAATACGGTGGACCGTTTTCTGTTCCGGTTCCACTTGATTTTATTCGGCATTGTTCATCTGCTGGAACGACa                                                                    <  6:115533/148‑1 (MQ=255)
                                                                              tctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCa                                                                    <  6:335876/148‑1 (MQ=255)
                                                                                  gCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTc                                                                            >  2:328703/1‑136 (MQ=255)
                                                                                  gCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTc                                                                            <  1:328703/136‑1 (MQ=255)
                                                                                  gCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCa                                                               >  1:81212/1‑149 (MQ=255)
                                                                                   ccACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCAc                                                              <  2:81212/149‑1 (MQ=255)
                                                                                    cACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGtt                                                                                                                           >  6:20135/1‑87 (MQ=255)
                                                                                    cACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCt                                                                                <  4:199355/130‑1 (MQ=255)
                                                                                    cACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCt                                                                                >  3:199355/1‑130 (MQ=255)
                                                                                          tcttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAGCCg                                                       >  1:357980/1‑149 (MQ=255)
                                                                                           cttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCgg                                                      >  6:260676/1‑149 (MQ=255)
                                                                                           cttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTGCTTCTTTAATACGGTTGCACGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                         <  8:355542/146‑1 (MQ=255)
                                                                                            ttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCgg                                                      >  2:162093/1‑148 (MQ=255)
                                                                                              cGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGt                                                     <  5:42266/147‑1 (MQ=255)
                                                                                                       tttGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTg                                          <  1:162093/149‑1 (MQ=255)
                                                                                                          gcagcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGc                                        <  4:305855/148‑1 (MQ=255)
                                                                                                          gcagcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATAATCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGc                                        <  5:260676/148‑1 (MQ=255)
                                                                                                           cagcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCgg                                      <  2:150389/149‑1 (MQ=255)
                                                                                                           cagcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCgg                                      <  8:383583/149‑1 (MQ=255)
                                                                                                            agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATc                                                                                 <  5:254307/105‑1 (MQ=255)
                                                                                                            agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGga                                     <  1:126229/149‑1 (MQ=255)
                                                                                                              caAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTcttctt                                                                                                                                            <  6:380125/44‑1 (MQ=255)
                                                                                                               aaaCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGaa                                  >  2:182933/1‑149 (MQ=255)
                                                                                                                 aCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGaaa                                               >  6:76315/1‑134 (MQ=255)
                                                                                                                 aCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGaa                                  >  1:273027/1‑147 (MQ=255)
                                                                                                                 aCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACa                                <  3:380124/149‑1 (MQ=255)
                                                                                                                     aTTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGc                                                                      >  8:196050/1‑107 (MQ=255)
                                                                                                                      ttttttCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGGTCCGGTTCCACTTCAATATATTAGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATg                           >  5:9109/1‑149 (MQ=255)
                                                                                                                         tttCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGtt                         <  3:186150/148‑1 (MQ=255)
                                                                                                                           tcttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGa                       <  1:306528/148‑1 (MQ=255)
                                                                                                                            cttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATc                                                                                 >  2:416940/1‑89 (MQ=255)
                                                                                                                            cttcGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAg                      <  5:170660/148‑1 (MQ=255)
                                                                                                                                     aaGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGa                       <  7:309164/138‑1 (MQ=255)
                                                                                                                                     aaGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGa                       >  8:309164/1‑138 (MQ=255)
                                                                                                                                      aGGGTAGAAAAGTTCTTCTTTAATACTGTGGCCCGTTTTCTTGTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGAATGTTGAGGCGCTGGttgg           >  2:335890/1‑147 (MQ=255)
                                                                                                                                          tAGAAAAGTTGTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGCTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACTTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAgg       <  6:141908/149‑1 (MQ=255)
                                                                                                                                          tAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAgg       >  8:372178/1‑149 (MQ=255)
                                                                                                                                            gAACAGTTCTGCTTTAATACGGTGGCCCGTTTTCTGGTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCTCGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcg     <  2:347561/149‑1 (MQ=255)
                                                                                                                                            gAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggc      <  4:34016/148‑1 (MQ=255)
                                                                                                                                             gaaaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg    >  3:13488/2‑149 (MQ=255)
                                                                                                                                              aaaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAggcgg    >  6:67248/1‑148 (MQ=255)
                                                                                                                                               aaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAg                      >  2:44923/1‑129 (MQ=255)
                                                                                                                                               aaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAg                      <  1:44923/129‑1 (MQ=255)
                                                                                                                                               aaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCa  >  4:48769/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCA  >  NC_000913/3964783‑3965074

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCAT  >  NC_000913/3964782‑3965075
                                                                                                                                                   |                                                                                                                                                  
GACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCT                                                                                                                                                 >  NB501016:230:HFVJLAFXY:3:11601:12180:18517/1‑151 (MQ=60)
  CACGATGACCATCTGCGGCCAGGTGGCCCCAGATCTCTTCACACCGGTGTTTGGTGTTGGCGAAAATAATCGCGCTGTCTGGCCACTCTTCTTCGAGCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTC                                                                                                                                               <  NB501016:230:HFVJLAFXY:1:21210:25761:2370/151‑1 (MQ=60)
  AACGATGACCATCTGCTGCCAGGTGGCCCCAGTTCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTC                                                                                                                                               <  NB501016:230:HFVJLAFXY:3:11602:13176:16688/151‑1 (MQ=60)
    CGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTT                                                                                                                                             >  NB501016:230:HFVJLAFXY:1:11104:3349:9627/1‑151 (MQ=60)
    CGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTT                                                                                                                                             <  NB501016:230:HFVJLAFXY:4:11604:22916:2208/151‑1 (MQ=60)
    CGATGAACATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTT                                                                                                                                             >  NB501016:230:HFVJLAFXY:2:11112:17434:5547/1‑151 (MQ=60)
     tgaagagatctggggccacctggcagCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                            >  NB501016:230:HFVJLAFXY:4:21512:6607:11032/27‑151 (MQ=60)
     GATGAGCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                            <  NB501016:230:HFVJLAFXY:3:21407:21271:8083/151‑1 (MQ=60)
     GATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                            <  NB501016:230:HFVJLAFXY:3:11510:4798:16229/151‑1 (MQ=60)
      ATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                            >  NB501016:230:HFVJLAFXY:4:21610:24453:6914/1‑150 (MQ=60)
        GACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAAT                                                                                                                                         <  NB501016:230:HFVJLAFXY:3:21412:11472:18873/151‑1 (MQ=60)
         ACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTTGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTTCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATA                                                                                                                                        <  NB501016:230:HFVJLAFXY:4:21610:24002:8245/151‑1 (MQ=60)
          CCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATAC                                                                                                                                       >  NB501016:230:HFVJLAFXY:2:11305:12667:1608/1‑151 (MQ=60)
              CTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATA                                                                                                                                        <  NB501016:230:HFVJLAFXY:4:11601:20449:5873/146‑1 (MQ=60)
                GCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGG                                                                                                                                  <  NB501016:230:HFVJLAFXY:2:11305:4032:3934/150‑1 (MQ=60)
                GCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGG                                                                                                                                  <  NB501016:230:HFVJLAFXY:2:21310:15771:14439/150‑1 (MQ=60)
                GCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTG                                                                                                                                   <  NB501016:230:HFVJLAFXY:4:21608:2893:20383/149‑1 (MQ=60)
                      AGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTT                                                                                                                           <  NB501016:230:HFVJLAFXY:2:11210:15683:17675/151‑1 (MQ=60)
                      AGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTT                                                                                                                           <  NB501016:230:HFVJLAFXY:3:21406:12388:13772/151‑1 (MQ=60)
                          GGCCCGAGATCTCTTAACAACGTTGTTTGGTGTTGGTTAATATACTCGATCTGTCTGGCCACTCTTATTCGATGAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTG                                                                                                                       <  NB501016:230:HFVJLAFXY:3:21511:24881:4933/151‑1 (MQ=60)
                           GCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGT                                                                                                                      >  NB501016:230:HFVJLAFXY:4:11606:8930:17440/1‑151 (MQ=60)
                           GCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTT                                                                                                                            <  NB501016:230:HFVJLAFXY:1:21112:1258:4837/145‑1 (MQ=60)
                            CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTT                                                                                                                     >  NB501016:230:HFVJLAFXY:2:11109:12506:16405/1‑151 (MQ=60)
                            CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACCGTGGCCCGTTTTCTGTT                                                                                                                     >  NB501016:230:HFVJLAFXY:3:11403:6569:14681/1‑151 (MQ=60)
                            CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTC                                                                                                                                               <  NB501016:230:HFVJLAFXY:4:11505:20045:15190/125‑1 (MQ=60)
                            CCCCAGATCTCTTCACAACGGTGTTTGGTGGTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                            <  NB501016:230:HFVJLAFXY:3:21410:1798:3512/128‑1 (MQ=60)
                             CCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTC                                                                                                                    >  NB501016:230:HFVJLAFXY:3:21512:13164:6682/1‑151 (MQ=60)
                                      CTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTACGGTTACA                                                                                                            >  NB501016:230:HFVJLAFXY:1:11212:18348:18119/1‑150 (MQ=60)
                                         CACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCA                                                                                                            <  NB501016:230:HFVJLAFXY:3:11609:3851:3272/147‑1 (MQ=60)
                                           CAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCA                                                                                                       >  NB501016:230:HFVJLAFXY:2:21210:23613:14707/1‑150 (MQ=60)
                                           CAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTATGTTCCGGTTCCACTTCA                                                                                                       >  NB501016:230:HFVJLAFXY:4:11611:6134:10116/1‑150 (MQ=60)
                                            ATCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGTCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAAA                                                                                                     <  NB501016:230:HFVJLAFXY:3:11612:6404:3346/151‑1 (MQ=60)
                                            AACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAAT                                                                                                     <  NB501016:230:HFVJLAFXY:2:21108:18943:7405/151‑1 (MQ=60)
                                             ACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATA                                                                                                    <  NB501016:230:HFVJLAFXY:4:21411:20663:17197/151‑1 (MQ=60)
                                              CGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGGCCGTTTTCTGTTCCGGTTCCACTTCAATAT                                                                                                   >  NB501016:230:HFVJLAFXY:3:21511:25687:5121/1‑151 (MQ=60)
                                               GGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATAT                                                                                                   >  NB501016:230:HFVJLAFXY:2:11104:6164:7494/1‑150 (MQ=60)
                                               GGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATA                                                                                                    >  NB501016:230:HFVJLAFXY:4:11410:6699:5129/1‑149 (MQ=60)
                                                  GTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTC                                                                                               >  NB501016:230:HFVJLAFXY:4:11511:19137:15837/1‑151 (MQ=60)
                                                  GTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACT                                                                                                          >  NB501016:230:HFVJLAFXY:3:11611:24918:11927/1‑140 (MQ=60)
                                                           tctcgataatatTCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTATAAGGGTAGAAAAGTTCTTCTTTAATACGGTGTCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGT                                                                                       <  NB501016:230:HFVJLAFXY:3:21610:25600:13974/138‑1 (MQ=60)
                                                             GCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTC                                                                                     >  NB501016:230:HFVJLAFXY:1:21312:7002:12721/1‑150 (MQ=60)
                                                              CGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCAT                                                                                   >  NB501016:230:HFVJLAFXY:4:21605:12532:7673/1‑151 (MQ=60)
                                                              CGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATA                                                                                                                                        >  NB501016:230:HFVJLAFXY:3:21511:10581:8953/1‑98 (MQ=60)
                                                               GAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATC                                                                                  <  NB501016:230:HFVJLAFXY:4:11412:15645:6391/151‑1 (MQ=60)
                                                               GAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATC                                                                                  <  NB501016:230:HFVJLAFXY:3:21504:5931:19560/151‑1 (MQ=60)
                                                                     AATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCG                                                                            <  NB501016:230:HFVJLAFXY:3:21507:22436:15990/151‑1 (MQ=60)
                                                                     AATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCG                                                                            <  NB501016:230:HFVJLAFXY:1:11107:23294:9015/151‑1 (MQ=60)
                                                                     AATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCG                                                                            <  NB501016:230:HFVJLAFXY:1:21307:2679:17965/151‑1 (MQ=60)
                                                                      ATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGA                                                                           <  NB501016:230:HFVJLAFXY:2:21204:14453:7566/151‑1 (MQ=60)
                                                                       TCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGA                                                                           >  NB501016:230:HFVJLAFXY:2:11107:23575:9475/1‑150 (MQ=60)
                                                                        CGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAAC                                                                         >  NB501016:230:HFVJLAFXY:2:21109:4793:8871/1‑151 (MQ=60)
                                                                        CGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTC                                                                                               <  NB501016:230:HFVJLAFXY:3:21601:12504:6018/129‑1 (MQ=60)
                                                                         GCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTATTACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACG                                                                        <  NB501016:230:HFVJLAFXY:4:21604:20251:19024/151‑1 (MQ=60)
                                                                            CTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCAT                                                                                          >  NB501016:230:HFVJLAFXY:2:11206:10541:7758/1‑130 (MQ=60)
                                                                             TGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCG                                                                                              >  NB501016:230:HFVJLAFXY:3:21412:6158:19834/1‑125 (MQ=60)
                                                                              GTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCC                                                                      <  NB501016:230:HFVJLAFXY:1:11101:15532:9299/148‑1 (MQ=60)
                                                                               TCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCAC                                                                                                           <  NB501016:230:HFVJLAFXY:3:11404:11824:7529/110‑1 (MQ=60)
                                                                                CTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTC                                                                 >  NB501016:230:HFVJLAFXY:4:11512:20471:3570/1‑151 (MQ=60)
                                                                                  GGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCAC                                                               >  NB501016:230:HFVJLAFXY:1:11202:23349:4088/1‑151 (MQ=60)
                                                                                  GGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCG                                                                            <  NB501016:230:HFVJLAFXY:1:21208:20961:14232/138‑1 (MQ=60)
                                                                                  GGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTT                                                                                      <  NB501016:230:HFVJLAFXY:4:21601:18492:5627/128‑1 (MQ=60)
                                                                                  GGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGGTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTT                                                                  <  NB501016:230:HFVJLAFXY:1:21312:6071:19320/148‑1 (MQ=60)
                                                                                    CCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACG                                                              <  NB501016:230:HFVJLAFXY:4:21510:6783:18932/150‑1 (MQ=60)
                                                                                    CCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTG                                                                                >  NB501016:230:HFVJLAFXY:2:21102:14354:2301/1‑132 (MQ=60)
                                                                                    CCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTT                                                                                                                           <  NB501016:230:HFVJLAFXY:3:11404:24774:11287/89‑1 (MQ=60)
                                                                                       CTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCGTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTTCC                                                          >  NB501016:230:HFVJLAFXY:3:11510:1380:4823/1‑151 (MQ=60)
                                                                                         CTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACC                                                          <  NB501016:230:HFVJLAFXY:1:11102:9958:2942/149‑1 (MQ=60)
                                                                                          TTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGATCGAACGCCAGTTCACGTAGCCGG                                                       >  NB501016:230:HFVJLAFXY:1:21301:1612:19009/1‑151 (MQ=60)
                                                                                           TCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGT                                                      >  NB501016:230:HFVJLAFXY:1:11302:21498:9977/1‑151 (MQ=60)
                                                                                            CTTCGATCGGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCAGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGT                                                      >  NB501016:230:HFVJLAFXY:1:11311:18669:6945/1‑150 (MQ=60)
                                                                                            CTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTTATACGGTGGCGCGTTTTCGGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATGTGCTCGATCGCCTGTTCACGTtcgcgtt                                                      >  NB501016:230:HFVJLAFXY:3:11501:26205:1308/1‑143 (MQ=60)
                                                                                            CTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCAGCTAGAACGCCAGTTCACGTACCCGGT                                                      >  NB501016:230:HFVJLAFXY:1:21212:15531:2714/1‑150 (MQ=60)
                                                                                              TCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTA                                                     <  NB501016:230:HFVJLAFXY:3:11408:7552:8676/149‑1 (MQ=60)
                                                                                               CGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACG                                                   >  NB501016:230:HFVJLAFXY:2:11303:17419:18162/1‑150 (MQ=60)
                                                                                                   gctcttcccgatctAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAG                                               <  NB501016:230:HFVJLAFXY:3:11502:23361:6920/136‑1 (MQ=60)
                                                                                                       GTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGG                                          <  NB501016:230:HFVJLAFXY:1:11304:19199:10505/151‑1 (MQ=60)
                                                                                                          gacagcgtcatcctggatcaTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGC                                         <  NB501016:230:HFVJLAFXY:2:21211:13114:18615/129‑1 (MQ=60)
                                                                                                          TGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCG                                        <  NB501016:230:HFVJLAFXY:2:11108:4554:16476/150‑1 (MQ=60)
                                                                                                          TGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATAATCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCG                                        <  NB501016:230:HFVJLAFXY:3:21410:19688:3849/150‑1 (MQ=60)
                                                                                                            CAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAG                                     <  NB501016:230:HFVJLAFXY:1:11210:9618:5300/151‑1 (MQ=60)
                                                                                                            CAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGG                                       >  NB501016:230:HFVJLAFXY:2:21310:11968:15356/1‑149 (MQ=60)
                                                                                                            CAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCT                                                                                 <  NB501016:230:HFVJLAFXY:3:21409:11301:1790/107‑1 (MQ=60)
                                                                                                              GCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTT                                                                                                                                            >  NB501016:230:HFVJLAFXY:3:21607:19105:15365/1‑46 (MQ=60)
                                                                                                                 AACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAG                                <  NB501016:230:HFVJLAFXY:2:21312:11323:3463/151‑1 (MQ=60)
                                                                                                                 AACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAAC                                  >  NB501016:230:HFVJLAFXY:1:21110:22567:17879/1‑149 (MQ=60)
                                                                                                                  ACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGC                               <  NB501016:230:HFVJLAFXY:4:11503:6551:11610/151‑1 (MQ=60)
                                                                                                                   CGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAG                                                                    <  NB501016:230:HFVJLAFXY:2:21101:10736:5985/113‑1 (MQ=60)
                                                                                                                     CATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCC                                                                      <  NB501016:230:HFVJLAFXY:4:11609:12687:8463/109‑1 (MQ=60)
                                                                                                                      ATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGGTCCGGTTCCACTTCAATATATTAGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGT                           >  NB501016:230:HFVJLAFXY:3:11402:4829:12356/1‑151 (MQ=60)
                                                                                                                       TTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGAGGAGAACAGCATGTT                          >  NB501016:230:HFVJLAFXY:3:11608:9267:13641/1‑151 (MQ=60)
                                                                                                                         TTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTG                         <  NB501016:230:HFVJLAFXY:2:11311:9281:14709/150‑1 (MQ=60)
                                                                                                                           TTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAG                       <  NB501016:230:HFVJLAFXY:1:21204:25856:14900/150‑1 (MQ=60)
                                                                                                                            TCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGG                      <  NB501016:230:HFVJLAFXY:3:11606:4467:15296/150‑1 (MQ=60)
                                                                                                                            TCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCT                                                                                 <  NB501016:230:HFVJLAFXY:1:21312:8591:4390/91‑1 (MQ=60)
                                                                                                                             CTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCG                    >  NB501016:230:HFVJLAFXY:3:21411:5775:6396/1‑151 (MQ=60)
                                                                                                                             CTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCG                    >  NB501016:230:HFVJLAFXY:3:21403:10820:16589/1‑151 (MQ=60)
                                                                                                                                     GAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTG            <  NB501016:230:HFVJLAFXY:3:21511:8340:8482/151‑1 (MQ=60)
                                                                                                                                     GAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAG                       <  NB501016:230:HFVJLAFXY:4:21504:14581:9519/140‑1 (MQ=60)
                                                                                                                                       AGGGTAGAAAAGTGCTGCTTTAATACGGTGGCCCGTTTTCTGGTCCGGTTCCTCTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGCGAACAGCATGTTGAGGCGCTGGTTTGCA          <  NB501016:230:HFVJLAFXY:1:11206:24943:3418/151‑1 (MQ=60)
                                                                                                                                            CGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGGTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGG     <  NB501016:230:HFVJLAFXY:4:11402:2632:1129/151‑1 (MQ=60)
                                                                                                                                             GGAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGC    >  NB501016:230:HFVJLAFXY:2:11103:11496:11054/1‑151 (MQ=60)
                                                                                                                                              AAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGG     <  NB501016:230:HFVJLAFXY:3:11604:11730:16912/149‑1 (MQ=60)
                                                                                                                                               AAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGG                      <  NB501016:230:HFVJLAFXY:1:11108:17101:2045/131‑1 (MQ=60)
                                                                                                                                               AAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGAATTGTTCATCTGCTCGAACGCCAGTTCACGTTCCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCAT  >  NB501016:230:HFVJLAFXY:1:21207:22170:19475/1‑151 (MQ=60)
                                                                                                                                                   |                                                                                                                                                  
GACACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCAT  >  NC_000913/3964782‑3965075

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: