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
A4 F27 I1 R1 56 84.2 3463790 98.6% 3415296 146.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 5,677,568 G→C intergenic (+115/+25) 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,5680GC97.6% 287.5 / ‑3.9 82intergenic (+115/+25)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base C (45/35);  total (47/35)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAA  >  NC_002947_CJ‑RC/5677427‑5677703
                                                                                                                                                |                                                                                                                                       
ccAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTAGTGATAGGGCCGCAAAGCCGcccccta                                                                                                                                  >  4:230794/1‑144 (MQ=255)
 cAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                   <  2:324221/148‑1 (MQ=255)
 cAGCGCGGCTGACCCGCAAGGGCTGATAGAAGTCCCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                   <  2:202556/148‑1 (MQ=37)
  aGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCgg                                                                                                                                  >  6:15014/1‑148 (MQ=255)
     gctgctGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGtt                                                                                                                                >  6:424893/1‑147 (MQ=255)
        gctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTCa                                                                                                                            >  4:307693/1‑148 (MQ=255)
           gACCCGCAAGGGCTGATAGAAGTACCCAGGGGGCTGCTTCGCAGCCCTATCACTACACAAGGCGGCTCCTGCAGAGATTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGGGGCCCCCTCGGTTTTCAGCt                                                                                                                         <  2:322095/148‑1 (MQ=37)
            aCCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTg                                                                                                                        >  6:75150/1‑148 (MQ=255)
               cGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCa                                                                                                                     <  4:340523/148‑1 (MQ=255)
                gCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGTATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCa                                                                                                                     >  7:430022/1‑147 (MQ=37)
                gCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGGTAGGGCCGCAAAGCGGCCCCCGCGGGTTTCAGCTggga                                                                                                                     >  3:90958/1‑145 (MQ=37)
                gCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCa                                                                                                                     >  4:124814/1‑147 (MQ=255)
                      gCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTa                                                                                                              <  5:93494/148‑1 (MQ=255)
                       cTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTAc                                                                                                             <  6:302264/148‑1 (MQ=255)
                           tAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCg                                                                                                                                   <  3:111361/122‑1 (MQ=255)
                                    cgcgggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCacac                                                                                                <  6:128333/145‑1 (MQ=255)
                                    cgcgggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTCCAGAGAATTTCTTCATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCacac                                                                                                <  1:222417/145‑1 (MQ=37)
                                     gAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACAcc                                                                                               >  1:407307/1‑148 (MQ=255)
                                      aGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCt                                                                                              >  2:414724/1‑148 (MQ=255)
                                      aGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCAGAAATTACTTATGATTCACAcc                                                                                               >  6:56066/1‑147 (MQ=37)
                                       gggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                             <  5:117901/148‑1 (MQ=255)
                                          ggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                          <  4:134997/148‑1 (MQ=255)
                                           gctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGac                                                                                         <  5:125057/148‑1 (MQ=255)
                                           gctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGac                                                                                         <  8:371471/148‑1 (MQ=255)
                                           gctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAAATTCTTGATCTCTACAGGAGCGGCCTTGAGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGATTCACACCTTCGac                                                                                         >  8:353518/1‑148 (MQ=37)
                                             tgctTCGCTGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTTTTCTCTACAGGAGCGGCCTTGTGTCTTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                       <  5:270745/148‑1 (MQ=37)
                                               ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCacac                                                                                                >  7:207880/1‑137 (MQ=255)
                                                ttCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACAcc                                                                                               <  8:207880/137‑1 (MQ=255)
                                                  cGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGata                                                                                  >  1:9706/1‑148 (MQ=255)
                                                     aGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGTAGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtcc                                                                               <  3:307693/148‑1 (MQ=38)
                                                     aGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACATCTGCTTCACACCTTCGACACTGatat                                                                                 >  4:219009/1‑146 (MQ=38)
                                                        ccTATCACGACACAAGGCCGCTCCTGCAGAGAAATTCTTGATCTCTACAGGAGCGGCATTGTGTCGTGATAGGGCCTCAAAGCGGCGCCCTCGGTTTTCACCTGGCAGAAATTACTTCTGCTTCACACCATCTACACTGATAtcaac                                                                             >  4:352954/1‑144 (MQ=37)
                                                         cTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggt                                                                           >  2:281947/1‑148 (MQ=255)
                                                               cGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                      <  5:424893/147‑1 (MQ=255)
                                                               cGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAggg                                                                     >  7:259387/1‑148 (MQ=255)
                                                                  cacaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCt                                                                  >  1:253380/1‑148 (MQ=255)
                                                                     aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                >  1:252362/1‑147 (MQ=255)
                                                                     aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                >  3:117395/1‑147 (MQ=255)
                                                                     aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                >  1:408852/1‑147 (MQ=255)
                                                                     aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCGTCGCCACTGATATCCAGGTCCAGGGTCTGc                                                                >  7:178412/1‑147 (MQ=39)
                                                                         ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                           <  7:257865/148‑1 (MQ=255)
                                                                         ccGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                           <  8:356012/148‑1 (MQ=255)
                                                                          cGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                          >  2:40888/1‑148 (MQ=255)
                                                                           gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTa                                                          >  4:272671/1‑147 (MQ=255)
                                                                             tCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCa                                                       >  2:310601/1‑148 (MQ=255)
                                                                               cTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCAcc                                                     >  3:70468/1‑148 (MQ=255)
                                                                                  cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                  >  8:202885/1‑148 (MQ=255)
                                                                                  cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGGATCACACCGTCGACACAGATATCCACGTCCAGGGTCTGCGAAGTAGCACCCgg                                                  >  7:421024/1‑148 (MQ=38)
                                                                                   agagAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCggg                                                 <  2:408852/148‑1 (MQ=255)
                                                                                   agagAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCCGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                  <  4:88388/147‑1 (MQ=255)
                                                                                    gagaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGc                                                <  3:272671/148‑1 (MQ=255)
                                                                                        atatCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCATAAAATACATCTTCAACACACATTCCACACTTATATCCACGTACACCGTCTGCGAAGAAGCACCCGGGCCtt                                             >  8:400831/4‑147 (MQ=19)
                                                                                          ttCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGa                                          <  1:414724/148‑1 (MQ=255)
                                                                                            cTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATg                                        >  2:77334/1‑148 (MQ=255)
                                                                                             ttGATCTCTACCGGAGGGGCTTTGTGGGTTTATTGGGCGGCAACGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATg                                        <  1:218665/147‑1 (MQ=255)
                                                                                                aTCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGcc                                      >  4:217420/1‑146 (MQ=255)
                                                                                                 tctcTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                   >  1:416280/1‑148 (MQ=255)
                                                                                                 tctcAACAGCAGCGGCCTTGTGCCGTGATTGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                   <  3:13162/148‑1 (MQ=255)
                                                                                                  ctctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                   <  4:70468/147‑1 (MQ=255)
                                                                                                   tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGt                                 >  3:400337/1‑148 (MQ=255)
                                                                                                    ctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                >  6:213268/1‑148 (MQ=255)
                                                                                                    ctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTATGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGAACCCGGGCGTTTGATGCCGAAGTc                                >  5:52631/1‑148 (MQ=38)
                                                                                                      ccAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAACGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGaa                                   <  5:131008/142‑1 (MQ=255)
                                                                                                      aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCTGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAg                                  <  2:383478/144‑1 (MQ=255)
                                                                                                      aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCg                               >  6:393620/1‑147 (MQ=255)
                                                                                                        aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgt                              <  5:396443/146‑1 (MQ=255)
                                                                                                         ggAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttca                           <  3:124814/148‑1 (MQ=255)
                                                                                                             cGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggtt                       <  5:15014/148‑1 (MQ=255)
                                                                                                               gCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttca                     >  6:128383/1‑148 (MQ=255)
                                                                                                                  ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGGTGGCAGAAATTACTTCTGCTTGACACCTTCGACACTGATATCCAGGTACAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTGAGGTTCAGCg                  >  5:238949/1‑148 (MQ=255)
                                                                                                                       tCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg             >  1:307190/1‑148 (MQ=255)
                                                                                                                         gTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggt           >  2:427605/1‑148 (MQ=255)
                                                                                                                          tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtggtg          >  7:11412/1‑148 (MQ=255)
                                                                                                                              aGGGCCGCAAAGCGGCCCCGTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCg       <  2:219559/147‑1 (MQ=255)
                                                                                                                                ggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt     >  2:62155/1‑147 (MQ=255)
                                                                                                                                ggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt     >  4:381099/1‑147 (MQ=255)
                                                                                                                                ggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt     >  6:241395/1‑147 (MQ=255)
                                                                                                                                ggCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt     >  6:236023/1‑147 (MQ=255)
                                                                                                                                 gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt     >  8:61079/1‑146 (MQ=255)
                                                                                                                                  ccGCAAAGTGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa  <  1:102710/148‑1 (MQ=255)
                                                                                                                                  ccGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa  <  5:219863/148‑1 (MQ=255)
                                                                                                                                  ccGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAAAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa  <  2:275174/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                       
CCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAA  >  NC_002947_CJ‑RC/5677427‑5677703

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
GCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCC  >  NC_002947_CJ‑RC/5677425‑5677718
                                                                                                                                                 |                                                                                                                                                      
GCCCAGCGCTGCTGACCCGCAAGGGCTTATAGAACTACCGAGGGGGATGCTTCGCAGCCATATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCG                                                                                                                                                  <  NB501016:282:H2V7GAFX2:4:21409:11420:1365/151‑1 (MQ=60)
GCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCG                                                                                                                                                  >  NB501016:282:H2V7GAFX2:3:11508:9348:7418/1‑151 (MQ=60)
            gcacccgcaagggctgaaccacgacccgcGGGGGCTGCTGCGCACCCCTACCACGACACACGCCCGCTCCTGCAGATAATTTCTTCACCTCACCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGG                                                                                                                                      <  NB501016:282:H2V7GAFX2:2:11106:18823:11154/122‑1 (MQ=39)
                 CGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGTATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCAGA                                                                                                                                  >  NB501016:282:H2V7GAFX2:4:21610:1977:5949/1‑150 (MQ=60)
                      ggtcttatacactccccgtggggtcggcttcccccccaaaccccgacccacgcccgctcctcactataatttcttcatctaacccgtcgggcCCTTGTGTCGTTATTGGGCCGAAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTAC                                                                                                                            <  NB501016:282:H2V7GAFX2:2:21305:6077:19942/59‑1 (MQ=26)
                      GGGCTTATAGACGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACACGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTAC                                                                                                                            <  NB501016:282:H2V7GAFX2:3:11504:26582:20189/151‑1 (MQ=60)
                           GATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGG                                                                                                                                                 <  NB501016:282:H2V7GAFX2:2:11207:19452:11486/125‑1 (MQ=60)
                              tgaagtaccgagggggctctttccctgccctatcacgacacagggccccaccggcatattttttcttgctgtctcctgcAGCGGCCTTGTGTCTTGATTGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCT                                                                                                                     <  NB501016:282:H2V7GAFX2:2:11305:22472:7888/71‑1 (MQ=46)
                                    cccgcGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTCCAGAGAATTTCTTCATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACC                                                                                                              <  NB501016:282:H2V7GAFX2:1:11312:17063:17536/146‑1 (MQ=60)
                                      CGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTT                                                                                                            >  NB501016:282:H2V7GAFX2:1:21308:20211:18518/1‑151 (MQ=60)
                                       GAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTC                                                                                                           <  NB501016:282:H2V7GAFX2:3:11508:7639:18162/151‑1 (MQ=60)
                                           GGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACA                                                                                                       <  NB501016:282:H2V7GAFX2:3:11510:9874:1579/151‑1 (MQ=60)
                                             GCTGCTTCGCTGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTTTTCTCTACAGGAGCGGCCTTGTGTCTTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACT                                                                                                     <  NB501016:282:H2V7GAFX2:3:21410:11958:2037/151‑1 (MQ=60)
                                                GCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCT                                                                                                             >  NB501016:282:H2V7GAFX2:4:11610:18653:3091/1‑140 (MQ=60)
                                                   TCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATC                                                                                               >  NB501016:282:H2V7GAFX2:1:11102:3418:9016/1‑151 (MQ=60)
                                                     caAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGTAGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCA                                                                                             <  NB501016:282:H2V7GAFX2:2:21205:18571:14437/149‑1 (MQ=60)
                                                               CACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGG                                                                                    <  NB501016:282:H2V7GAFX2:3:21611:16901:15555/150‑1 (MQ=60)
                                                                ACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTC                                                                                  >  NB501016:282:H2V7GAFX2:4:21406:11558:10982/1‑151 (MQ=60)
                                                                   ACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGC                                                                               >  NB501016:282:H2V7GAFX2:1:21105:25136:15620/1‑151 (MQ=60)
                                                                      CAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGA                                                                             >  NB501016:282:H2V7GAFX2:2:11208:17852:11026/1‑150 (MQ=60)
                                                                      CAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGA                                                                             >  NB501016:282:H2V7GAFX2:1:21309:22257:4991/1‑150 (MQ=60)
                                                                      CAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGA                                                                             >  NB501016:282:H2V7GAFX2:1:21105:24408:12586/1‑150 (MQ=60)
                                                                      CAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCGTCGCCACTGATATCCAGGTCCAGGGTCTGCGA                                                                             >  NB501016:282:H2V7GAFX2:4:11605:15355:5423/1‑150 (MQ=60)
                                                                         GGCCGCTCCTGCAGAGAATTTCTTGATCTCTACTGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTA                                                                         <  NB501016:282:H2V7GAFX2:4:21406:17942:6234/151‑1 (MQ=60)
                                                                            CGCTCATGCAGATAATTTCTTTATCTCTACAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAG                                                                                            <  NB501016:282:H2V7GAFX2:3:11409:18952:12434/129‑1 (MQ=60)
                                                                                CCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACAATGATATCCAGGTCCAGGGTCTGCGAAGTATCACCCC                                                                  >  NB501016:282:H2V7GAFX2:2:11112:17092:16748/1‑151 (MQ=60)
                                                                                   GCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGGATCACACCGTCGACACAGATATCCACGTCCAGGGTCTGCGAAGTAGCACCCGGGC                                                               >  NB501016:282:H2V7GAFX2:4:21608:11795:15687/1‑151 (MQ=60)
                                                                                    AAGAGAATTTCTTGATCTCTACAGGTGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCC                                                              <  NB501016:282:H2V7GAFX2:2:21109:14586:16787/151‑1 (MQ=60)
                                                                                          ATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAT                                                        <  NB501016:282:H2V7GAFX2:1:21310:23841:5951/151‑1 (MQ=60)
                                                                                                 CATCTCAACAGCAGCGGCCTTGTGCCGTGATTGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAG                                                 <  NB501016:282:H2V7GAFX2:2:11103:10268:5345/151‑1 (MQ=60)
                                                                                                  ATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGT                                                >  NB501016:282:H2V7GAFX2:1:21310:4020:11278/1‑151 (MQ=60)
                                                                                                    CTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCG                                              >  NB501016:282:H2V7GAFX2:2:21309:26366:9459/1‑151 (MQ=60)
                                                                                                     TCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTATGCTTCACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGAACCCGGGCGTTTGATGCCGAAGTCGT                                             >  NB501016:282:H2V7GAFX2:3:11410:21730:5275/1‑151 (MQ=60)
                                                                                                      CTCCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAACGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAG                                                 <  NB501016:282:H2V7GAFX2:3:11510:25480:20199/146‑1 (MQ=60)
                                                                                                        ACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTT                                            <  NB501016:282:H2V7GAFX2:3:21606:13018:18575/149‑1 (MQ=60)
                                                                                                         CAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAG                                         <  NB501016:282:H2V7GAFX2:2:11209:1386:14573/151‑1 (MQ=60)
                                                                                                             AGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTC                                     <  NB501016:282:H2V7GAFX2:3:11403:24091:12480/151‑1 (MQ=60)
                                                                                                                   CTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGGTGGCAGAAATTACTTCTGCTTGACACCTTCGACACTGATATCCAGGTACAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTGAGGTTCAGCGTG                               >  NB501016:282:H2V7GAFX2:3:21404:17343:4049/1‑151 (MQ=60)
                                                                                                                        GTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGT                          >  NB501016:282:H2V7GAFX2:1:21202:11738:19770/1‑151 (MQ=60)
                                                                                                                           GTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGC                       >  NB501016:282:H2V7GAFX2:4:11402:15474:16939/1‑151 (MQ=60)
                                                                                                                                  GGCCGCAAAGTGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAG                <  NB501016:282:H2V7GAFX2:1:11204:1232:15320/151‑1 (MQ=60)
                                                                                                                                  GGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAG                <  NB501016:282:H2V7GAFX2:3:11612:7550:11905/151‑1 (MQ=60)
                                                                                                                                               GCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCTCCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCG    <  NB501016:282:H2V7GAFX2:3:11611:8336:11096/151‑1 (MQ=60)
                                                                                                                                                 CCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCC  >  NB501016:282:H2V7GAFX2:1:21204:19352:1234/1‑151 (MQ=60)
                                                                                                                                                 |                                                                                                                                                      
GCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCC  >  NC_002947_CJ‑RC/5677425‑5677718

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: