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
A2 F43 I1 R1 23 34.4 1194698 96.8% 1156467 139.5

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,9290CT100.0% 177.3 / NA 50E234E (GAG→GAArhlBATP‑dependent RNA helicase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (24/26);  total (24/26)

ACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCA  >  NC_000913/3964785‑3965074
                                                                                                                                                |                                                                                                                                                 
acGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtctt                                                                                                                                               <  6:98950/149‑1 (MQ=255)
acGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTtctt                                                                                                                                               >  5:98950/1‑149 (MQ=255)
  gATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTTGTTAGAAGGGTAGAAAAGTtcttct                                                                                                                                             >  5:46544/1‑149 (MQ=255)
       ccATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttc                                                                                                                                                 >  1:25905/1‑140 (MQ=255)
       ccATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGttc                                                                                                                                                 <  2:25905/140‑1 (MQ=255)
        cATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATa                                                                                                                                       >  5:32117/1‑149 (MQ=255)
        cATCTGCTGCCAGGTGGCCCCAGATCACTTCACAACGGGCTTTGGTGTTGGCCAAAATAATCACACTGTCTAGCCACTCTTCTTCGAACAGCGTTAGCACCAAACCCATTTTTTCTTCGTTAGAAAGGTAGAAAAGTTCTTCTTTAATa                                                                                                                                       >  8:67502/1‑149 (MQ=255)
         aTCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATAc                                                                                                                                      >  4:70138/1‑149 (MQ=255)
                 ccAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGccc                                                                                                                              >  1:167445/1‑149 (MQ=255)
                    ggTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGtt                                                                                                                           <  5:143677/149‑1 (MQ=255)
                      gggCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGtttt                                                                                                                         <  6:32117/148‑1 (MQ=255)
                       ggCACCAGATCTCTTCACAACTGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTGTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGtttt                                                                                                                         <  7:67502/148‑1 (MQ=255)
                          cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCt                                                                                                                       >  1:23775/1‑147 (MQ=255)
                          cccAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCt                                                                                                                       <  2:23775/147‑1 (MQ=255)
                                  tctTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTc                                                                                                             >  1:49842/1‑149 (MQ=255)
                                    ttCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCgg                                                                                                                <  2:79872/144‑1 (MQ=255)
                                    ttCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCa                                                                                                           >  6:150936/1‑149 (MQ=255)
                                    ttCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCa                                                                                                           >  1:31227/1‑149 (MQ=255)
                                    ttCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCa                                                                                                           >  7:13384/1‑149 (MQ=255)
                                         aaCGGTGTTTGGTGTTGGCGAAAATACTCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGtttt                                                                                                                         <  3:91707/130‑1 (MQ=255)
                                         aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGtttt                                                                                                                         >  4:91707/1‑130 (MQ=255)
                                         aaCGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGCCCGTTTTCTGTTCCGGTTCCACTTCaa                                                                                                     <  4:34589/149‑1 (MQ=255)
                                           cGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAat                                                                                                    <  4:106488/149‑1 (MQ=255)
                                                tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTTGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATAtt                                                                                               <  6:46544/149‑1 (MQ=255)
                                                tttGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATAtt                                                                                               >  1:1633/1‑149 (MQ=255)
                                                  tGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTc                                                                                              <  3:70138/148‑1 (MQ=255)
                                                    gtgtTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCgg                                                                                            >  1:66269/1‑148 (MQ=255)
                                                          gCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGt                                                                                      >  2:79100/1‑148 (MQ=255)
                                                            gAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCa                                                                                                      >  6:66543/1‑130 (MQ=255)
                                                            gAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCa                                                                                                      <  5:66543/130‑1 (MQ=255)
                                                                      gCTCTGTCTGGCCACTCCTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGaa                                                                         <  7:23487/149‑1 (MQ=255)
                                                                       ctctGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAAc                                                                        <  1:79100/149‑1 (MQ=255)
                                                                       ctctGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAAc                                                                        <  2:167445/149‑1 (MQ=255)
                                                                          tgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCgg                                                                                                                <  8:113099/106‑1 (MQ=255)
                                                                          tgtctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCgg                                                                                                                >  7:113099/1‑106 (MQ=255)
                                                                            tctgGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGc                                                                      <  2:31227/146‑1 (MQ=255)
                                                                                gCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGtt                                                                 <  1:120117/147‑1 (MQ=255)
                                                                                 ccACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCAc                                                              <  5:150936/149‑1 (MQ=255)
                                                                                 ccACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCAc                                                              <  2:49842/149‑1 (MQ=255)
                                                                                     tcttctTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGTGTAGAAAAGTTCTTCTTTAATACGGTGGCACGTTTTCTGTTCCGGTTCCACTTCAATATATTATGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAc                                                          >  7:28593/1‑149 (MQ=255)
                                                                                         cttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTAcc                                                         <  2:66269/146‑1 (MQ=255)
                                                                                          ttcGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCgg                                                      >  2:17917/1‑148 (MQ=255)
                                                                                                 aGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGaa                                                <  5:68979/147‑1 (MQ=255)
                                                                                                 aGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGaa                                                >  6:68979/1‑147 (MQ=255)
                                                                                                  gCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGc                                             <  8:13384/149‑1 (MQ=255)
                                                                                                          agcaAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCg                                       >  3:65531/1‑147 (MQ=255)
                                                                                                             aaaCGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGaa                                  <  4:65531/149‑1 (MQ=255)
                                                                                                                 gCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTAAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAg                               >  2:163883/1‑148 (MQ=255)
                                                                                                                       tttCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGtt                         <  1:163883/148‑1 (MQ=255)
                                                                                                                                     gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGc                                                                      <  1:158062/89‑1 (MQ=255)
                                                                                                                                     gggTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGc                                                                      >  2:158062/1‑89 (MQ=255)
                                                                                                                                             aaGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCa  >  1:22439/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
ACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCA  >  NC_000913/3964785‑3965074

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ 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
CACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCAT  >  NC_000913/3964784‑3965075
                                                                                                                                                 |                                                                                                                                                  
CACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTC                                                                                                                                               >  NB501016:228:HFK2FAFXY:3:21410:4076:11952/1‑151 (MQ=60)
  CGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTTGTTAGAAGGGTAGAAAAGTTCTTCTT                                                                                                                                             >  NB501016:228:HFK2FAFXY:3:11510:3867:12180/1‑151 (MQ=60)
  CGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGT                                                                                                                                                    <  NB501016:228:HFK2FAFXY:1:21212:7810:18434/144‑1 (MQ=60)
       ACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCT                                                                                                                                                 >  NB501016:228:HFK2FAFXY:1:11111:24408:10043/1‑142 (MQ=60)
        CCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATAC                                                                                                                                       >  NB501016:228:HFK2FAFXY:3:11503:2644:18937/1‑151 (MQ=60)
                 GCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGCAAAGTTCTTCTTTAATACGGTGGCCCG                                                                                                                              >  NB501016:228:HFK2FAFXY:1:21311:5656:11700/1‑151 (MQ=60)
                    AGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTT                                                                                                                           <  NB501016:228:HFK2FAFXY:3:21607:5811:15592/151‑1 (MQ=60)
                       TGGCACCAGATCTCTTCACAACTGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTGTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTC                                                                                                                         <  NB501016:228:HFK2FAFXY:4:21401:26113:2291/150‑1 (MQ=60)
                          CCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTG                                                                                                                       >  NB501016:228:HFK2FAFXY:1:11110:13239:12707/1‑149 (MQ=60)
                                  CTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCC                                                                                                             >  NB501016:228:HFK2FAFXY:1:11209:6296:9980/1‑151 (MQ=60)
                                    CTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCAC                                                                                                           >  NB501016:228:HFK2FAFXY:1:11201:12044:14246/1‑151 (MQ=60)
                                    CTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCAC                                                                                                           >  NB501016:228:HFK2FAFXY:4:11407:18127:15819/1‑151 (MQ=60)
                                         CAACGGTGTTTGGTGTTGGCGAAAATACTCGCGCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTC                                                                                                                         <  NB501016:228:HFK2FAFXY:2:21112:5789:18553/132‑1 (MQ=60)
                                                GTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTC                                                                                               >  NB501016:228:HFK2FAFXY:1:11101:19221:13734/1‑151 (MQ=60)
                                                  TTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCG                                                                                              <  NB501016:228:HFK2FAFXY:2:11312:24341:16252/150‑1 (MQ=60)
                                                    GGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGC                                                                                            >  NB501016:228:HFK2FAFXY:1:11304:8122:9386/1‑150 (MQ=60)
                                                            CGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAA                                                                                                      <  NB501016:228:HFK2FAFXY:3:11607:3999:8777/132‑1 (MQ=60)
                                                                      CGCTCTGTCTGGCCACTCCTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAAC                                                                         <  NB501016:228:HFK2FAFXY:4:11412:12041:8604/151‑1 (MQ=60)
                                                                       GCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACG                                                                        <  NB501016:228:HFK2FAFXY:1:11309:2059:15990/151‑1 (MQ=60)
                                                                          CTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGT                                                                                                                >  NB501016:228:HFK2FAFXY:4:21601:7220:7156/1‑108 (MQ=60)
                                                                                GGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTC                                                                 <  NB501016:228:HFK2FAFXY:1:21203:25002:4348/149‑1 (MQ=60)
                                                                                 GCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACG                                                              <  NB501016:228:HFK2FAFXY:3:21611:11267:6483/151‑1 (MQ=60)
                                                                                     CTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGTGTAGAAAAGTTCTTCTTTAATACGGTGGCACGTTTTCTGTTCCGGTTCCACTTCAATATATTATGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACC                                                          >  NB501016:228:HFK2FAFXY:4:11503:1311:3247/1‑151 (MQ=60)
                                                                                                 CAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAA                                                <  NB501016:228:HFK2FAFXY:3:11608:26221:10422/149‑1 (MQ=60)
                                                                                                          CAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGG                                       >  NB501016:228:HFK2FAFXY:2:11310:25086:14255/1‑149 (MQ=60)
                                                                                                                       TTTTCTTCGTTAGAAGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTG                         <  NB501016:228:HFK2FAFXY:1:21310:20377:1519/150‑1 (MQ=60)
                                                                                                                                     AGGGTAGAAAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCC                                                                      <  NB501016:228:HFK2FAFXY:1:21307:23333:10389/91‑1 (MQ=60)
                                                                                                                                     AGGGTAGAAAAGTTCTTCTTAAATACGGTGGCCCGTTTTCTGTTCC                                                                                                                   <  NB501016:228:HFK2FAFXY:4:11507:2970:7465/46‑1 (MQ=60)
                                                                                                                                             AAAGTTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCAT  >  NB501016:228:HFK2FAFXY:1:11110:24454:1698/1‑151 (MQ=60)
                                                                                                                                                 |                                                                                                                                                  
CACGATGACCATCTGCTGCCAGGTGGCCCCAGATCTCTTCACAACGGTGTTTGGTGTTGGCGAAAATAATCGCTCTGTCTGGCCACTCTTCTTCGATCAGCGTTTGCAGCAAACGCATTTTTTCTTCGTTAGAAGGGTAGAAAAGCTCTTCTTTAATACGGTGGCCCGTTTTCTGTTCCGGTTCCACTTCAATATATTCGGCATTGTTCATCTGCTCGAACGCCAGTTCACGTACCCGGTACGAAAGCGTGGCGGAGAACAGCATGTTGAGGCGCTGGTTTGCAGGCGGCAT  >  NC_000913/3964784‑3965075

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