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
A5 F18 I1 R1 59 77.7 3515218 98.5% 3462489 145.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 5,677,566 +G intergenic (+113/+27) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC5,677,5651.G100.0% 156.0 / NA 44intergenic (+112/+28)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (0/0);  new base G (25/19);  total (25/19)

TGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677420‑5677707
                                                                                                                                                  |                                                                                                                                              
tGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcggc                                                                                                                                             <  6:223889/149‑3 (MQ=255)
       ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                       >  2:77707/1‑149 (MQ=255)
        ccgcgcGGCCGACCCGCAAGGGCTGATAGACGTACCGCGGGGGGTGCTTCGCAGCCCTAACACCACACAAGGCCGCTCCTGCCGAGAATTTCTTGATCTCTACAGGAGGGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCgg                                                                                                                                      <  6:25643/147‑1 (MQ=25)
               gctgACCCGCAAGGGCTGATAGTAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcg                                                                                                                                               <  4:37781/132‑1 (MQ=255)
               gctgACCCGCAAGGGCTGATAGTAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcg                                                                                                                                               >  3:37781/1‑132 (MQ=255)
                ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCGGCAGAGAATTTCTTGATGTCTACAGGAGCGGCCTTGTGTCGTGGTAGGGCCGCAAAGCGGCCCCCGGGGTTTTCAGc                                                                                                                              >  6:286068/1‑149 (MQ=37)
                   aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTg                                                                                                                            <  8:51760/148‑1 (MQ=255)
                     ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGATTTTCGTGATCTCTAGAGGAGCGGCCTTGTGTCGTGATTGGGCCGGCAAGCGGCGCGCGGGGGGTTCAGCTgg                                                                                                                           >  4:63798/1‑147 (MQ=25)
                     ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTgg                                                                                                                           >  5:182122/1‑147 (MQ=255)
                         aaGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGaaa                                                                                                                     <  4:175091/149‑1 (MQ=255)
                                      aGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTc                                                                                                        <  7:253753/149‑1 (MQ=255)
                                        tACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTcac                                                                                                      <  2:345902/149‑1 (MQ=255)
                                        tACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTca                                                                                                       >  8:381730/1‑148 (MQ=255)
                                         aCCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTcaca                                                                                                     >  5:104699/1‑149 (MQ=255)
                                            gAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCt                                                                                                  >  4:215818/1‑149 (MQ=255)
                                              tagggCTGCGAAGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                                >  5:85825/3‑149 (MQ=37)
                                                      ctTCGCAGCCCTATCACGACACAAGGCCGGTCGTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTg                                                                                         <  3:55833/148‑1 (MQ=38)
                                                      ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGa                                                                                        <  6:110405/149‑1 (MQ=255)
                                                                cTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTCTGTCGTGATAGGGCCCCAAAGCGGCCCCCACGGTTATCACCTGGCAGAAAATACTTCTGCTTCACACCTTCGACACAGATAaccaggag                                                                              >  1:44490/1‑147 (MQ=37)
                                                                        acacAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCt                                                                      >  6:299467/1‑149 (MQ=255)
                                                                          acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                    >  8:158382/1‑149 (MQ=255)
                                                                           caAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                    >  3:103457/1‑148 (MQ=255)
                                                                                 cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                >  4:121844/1‑146 (MQ=255)
                                                                                 cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                <  3:121844/146‑1 (MQ=255)
                                                                                         cAGATAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                      <  5:286068/148‑1 (MQ=255)
                                                                                              aaTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCttt                                                <  5:163289/149‑1 (MQ=255)
                                                                                                  tCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAAATACTTCTGCTTCACACCTTCGACACTGATATCCAGCTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAt                                             <  6:85825/148‑1 (MQ=39)
                                                                                                   cTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATg                                            <  1:77707/148‑1 (MQ=255)
                                                                                                    ttGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGcc                                          >  6:429610/1‑149 (MQ=255)
                                                                                                      gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCg                                         <  4:103457/148‑1 (MQ=255)
                                                                                                            tACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACACGGGCCTTTGATGCCGAAg                                      <  6:138073/145‑1 (MQ=255)
                                                                                                               aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttca                               <  3:162390/149‑1 (MQ=255)
                                                                                                                 gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                             >  4:366430/1‑149 (MQ=255)
                                                                                                                 gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                             >  3:392903/1‑149 (MQ=255)
                                                                                                                     ggCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttca                         <  5:144506/149‑1 (MQ=255)
                                                                                                                     ggCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTACGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggtt                           >  3:51034/1‑147 (MQ=255)
                                                                                                                      gCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttcag                        >  4:349684/1‑149 (MQ=255)
                                                                                                                            tgtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggt                  >  5:329102/1‑149 (MQ=255)
                                                                                                                             gtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg                 >  3:161620/1‑149 (MQ=255)
                                                                                                                                gTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg              <  4:392903/149‑1 (MQ=255)
                                                                                                                                  gATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGc            >  8:183650/1‑149 (MQ=255)
                                                                                                                                     aGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt         >  6:151726/1‑149 (MQ=255)
                                                                                                                                      gggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTGCAGGTTCAGCGTGGTGGTGGCGtt         >  2:212037/1‑148 (MQ=255)
                                                                                                                                          cGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc    >  1:61432/1‑149 (MQ=255)
                                                                                                                                            cAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc    >  2:2464/1‑147 (MQ=255)
                                                                                                                                            cAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCa  <  7:59797/149‑1 (MQ=255)
                                                                                                                                            cAAAGCGGCCCCCGCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATCTCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGGCGTTCAGGTTCAGCGGGGTGGGGGCGTTGAAGCCa  <  6:182122/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                              
TGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947_CJ‑RC/5677420‑5677707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ 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
TGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCG  >  NC_002947_CJ‑RC/5677434‑5677710
                                                                                                                                     |                                                                                                                                                
TGCTGACCCGCAAGGGCTGATAGTAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG                                                                                                                                                 >  NB501016:282:H2V7GAFX2:2:11107:2531:8433/1‑134 (MQ=60)
    CACCCGCAAGGGCTGATTGATGTACCGTGGGGGCTGCTTCGACGCCCTATCACGACACAAGGCCGCTCCTGCAGATAATTTCTTGATCTCTACCGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGC                                                                                                                             <  NB501016:282:H2V7GAFX2:1:21202:15341:7930/151‑1 (MQ=60)
      CCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGC                                                                                                                             >  NB501016:282:H2V7GAFX2:3:11606:12831:9258/1‑149 (MQ=60)
                      GAAGTACCGAGGGGGCTGCTTCGCAGCCCAATAACGACACAAGGCCGCTCCTGCAGAGGATTTCTTGTTCTCTACAGGAGCGGCCTTGTGTGGGGATAGGGCCGCAAAGCGgcccccgggggtttgaggggggagaaatttattatggttt                                                                                                          >  NB501016:282:H2V7GAFX2:1:11302:22190:16513/1‑111 (MQ=60)
                       AAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCA                                                                                                          <  NB501016:282:H2V7GAFX2:4:21404:7818:19257/151‑1 (MQ=60)
                          TACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACAC                                                                                                       >  NB501016:282:H2V7GAFX2:3:11506:6345:5918/1‑151 (MQ=60)
                               atagggctgcgaaGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTGTGCTTCACACCTTCG                                                                                                  >  NB501016:282:H2V7GAFX2:3:11503:7072:5524/14‑151 (MQ=60)
                                  ttcagacgtgggttcttcgcatcaacacatgccggcccctgctgatattttcttctttTCTACAGTAGCGGCTTTGTGTCTTGATTGGGCCGCAACGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTTCTTCTGCTTCACACCTTCGAC                                                                                                <  NB501016:282:H2V7GAFX2:4:21402:1733:14270/92‑1 (MQ=54)
                                       GCTTCGCAGCCCTATCACGACACAAGGCCGGTCGTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGA                                                                                           <  NB501016:282:H2V7GAFX2:2:11110:17994:8024/150‑1 (MQ=60)
                                                 CCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTCTGTCGTGATAGGGCCCCAAAGCGGCCCCCACGGTTATCACCTGGCAGAAAATACTTCTGCTTCACACCTTCGACACAGATAACCAGGagc                                                                                >  NB501016:282:H2V7GAFX2:1:11107:7624:15327/1‑148 (MQ=60)
                                                          ACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCCCAAAGCGGCCCCCTCGGTTTTGAGCTGGCAGAAATTACATCTGCTTgacacaatcgacaaagataaacatgtacacggtctgg                                                                       >  NB501016:282:H2V7GAFX2:1:11203:20542:7094/1‑114 (MQ=60)
                                                            ACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCG                                                                      >  NB501016:282:H2V7GAFX2:2:11206:1701:5082/1‑150 (MQ=60)
                                                                  CCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGT                                                                  <  NB501016:282:H2V7GAFX2:2:11209:17535:3989/148‑1 (MQ=60)
                                                                          GCAGATAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGG                                                        <  NB501016:282:H2V7GAFX2:3:21411:19860:15005/150‑1 (MQ=60)
                                                                               GAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTG                                                  <  NB501016:282:H2V7GAFX2:3:11603:26272:11509/151‑1 (MQ=60)
                                                                                    TCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGC                                              <  NB501016:282:H2V7GAFX2:1:11112:19785:17295/150‑1 (MQ=60)
                                                                                                CAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAG                                 <  NB501016:282:H2V7GAFX2:2:11303:15656:12449/151‑1 (MQ=60)
                                                                                                  GGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGT                               >  NB501016:282:H2V7GAFX2:2:21307:17491:17324/1‑151 (MQ=60)
                                                                                                      CGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAG                           <  NB501016:282:H2V7GAFX2:3:11512:4794:13294/151‑1 (MQ=60)
                                                                                                      CGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTACGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTC                             >  NB501016:282:H2V7GAFX2:2:11109:1770:11700/1‑149 (MQ=60)
                                                                                                             GTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTAAGCGTGGTG                    >  NB501016:282:H2V7GAFX2:3:21506:13905:20123/1‑151 (MQ=60)
                                                                                                              TGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGG                   >  NB501016:282:H2V7GAFX2:2:11303:24067:10183/1‑151 (MQ=60)
                                                                                                                           CCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCC      >  NB501016:282:H2V7GAFX2:1:11110:21190:7936/1‑151 (MQ=60)
                                                                                                                             GCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAG    <  NB501016:282:H2V7GAFX2:1:21305:4200:14964/151‑1 (MQ=60)
                                                                                                                             GCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAG    <  NB501016:282:H2V7GAFX2:4:11409:12488:19931/151‑1 (MQ=60)
                                                                                                                               AAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCG  <  NB501016:282:H2V7GAFX2:3:11601:13237:15772/151‑1 (MQ=60)
                                                                                                                               AAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCG  <  NB501016:282:H2V7GAFX2:1:11101:16230:8584/151‑1 (MQ=60)
                                                                                                                                     |                                                                                                                                                
TGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCG  >  NC_002947_CJ‑RC/5677434‑5677710

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