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
A3 F23 I1 R1 59 87.9 3762788 98.7% 3713871 145.1

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 4,888,761 (A)6→5 coding (117/1776 nt) gcl → glyoxylate carboligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,888,7560A.97.5% 361.9 / 3.3 81coding (112/1776 nt)gclglyoxylate carboligase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (42/37);  total (43/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

TATACAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCAC  >  NC_002947_CJ‑RC/4888613‑4888898
                                                                                                                                               |                                                                                                                                              
tataCAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGaaaa                                                                                                                                             <  1:121478/147‑1 (MQ=255)
tataCAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGaaaa                                                                                                                                             <  5:13379/147‑1 (MQ=255)
 ataCAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGaaaaag                                                                                                                                           >  1:119632/1‑147 (MQ=255)
  taCAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGaaaaagt                                                                                                                                          <  2:79305/148‑3 (MQ=255)
    cAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGaaaaagtcg                                                                                                                                        >  3:175767/1‑144 (MQ=255)
     aaTTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCgg                                                                                                                                      <  7:470689/148‑1 (MQ=255)
      aTTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGt                                                                                                                                     >  1:127874/1‑148 (MQ=255)
       ttACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTg                                                                                                                                    >  6:57145/1‑148 (MQ=255)
           aaaaTCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCAt                                                                                                                                >  4:42735/1‑148 (MQ=255)
            aaaTCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATc                                                                                                                               <  7:191724/148‑1 (MQ=255)
              aTCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTACTGGTCATCCGCCGAGAAGGTGTAGAAACCGCGTTAGGCATCCCGGGGCCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGa                                                                                                                             >  5:87102/1‑148 (MQ=255)
               tCGATGAGGCACAATCCTTGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTTAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGAt                                                                                                                            <  5:5653/148‑1 (MQ=255)
               tCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGAt                                                                                                                            >  4:160864/1‑148 (MQ=255)
                   tGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACg                                                                                                                        >  4:359333/1‑148 (MQ=255)
                      ggCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTcc                                                                                                                     <  2:206172/148‑1 (MQ=255)
                          caAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGAt                                                                                                                            <  1:330535/137‑1 (MQ=255)
                          caAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGAt                                                                                                                            >  2:330535/1‑137 (MQ=255)
                          caAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGAATCGATCACGTCctcgc                                                                                                                 >  2:348027/1‑148 (MQ=255)
                           aaaCCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCctcgct                                                                                                                <  1:111614/148‑1 (MQ=255)
                           aaaCCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCctcgct                                                                                                                <  1:104915/148‑1 (MQ=255)
                              ccATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGt                                                                                                             <  3:36899/148‑1 (MQ=255)
                                   aGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGt                                                                                                        >  8:335184/1‑148 (MQ=255)
                                     cAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTTGGTGGCATCGATCACGTCCTCGCTCGTCACGTc                                                                                                       >  1:433744/1‑147 (MQ=255)
                                        aaTGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAg                                                                                                   <  1:437709/148‑1 (MQ=255)
                                                 aaTCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCg                                                                                           >  2:350494/1‑147 (MQ=255)
                                                 aaTCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCg                                                                                           >  5:73443/1‑147 (MQ=255)
                                                  aTCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGca                                                                                         <  7:217284/148‑1 (MQ=255)
                                                    cGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGcaca                                                                                       >  5:406464/1‑148 (MQ=255)
                                                     gATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACAt                                                                                      <  3:160864/148‑1 (MQ=255)
                                                         cAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGc                                                                                   <  4:102433/147‑1 (MQ=255)
                                                           gCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGcc                                                                                  >  4:450798/1‑146 (MQ=255)
                                                           gCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGcc                                                                                  >  2:68758/1‑146 (MQ=255)
                                                             cGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAgg                                                                              <  4:4161/148‑1 (MQ=255)
                                                             cGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAgg                                                                              <  4:175767/148‑1 (MQ=255)
                                                             cGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAgg                                                                              <  4:298786/148‑1 (MQ=255)
                                                              gTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAggg                                                                             <  8:311815/148‑1 (MQ=255)
                                                               ttCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGc                                                                            >  3:119679/1‑148 (MQ=255)
                                                                    gTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTacac                                                                       >  1:121104/1‑148 (MQ=255)
                                                                    gTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTac                                                                         <  7:291847/146‑1 (MQ=255)
                                                                     tCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACAcc                                                                      >  2:321886/1‑148 (MQ=255)
                                                                        tGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCGCGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCcgc                                                                   >  3:312143/1‑148 (MQ=255)
                                                                          cgcCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCgcgc                                                                 >  7:215865/1‑148 (MQ=255)
                                                                            ccGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGcc                                                                >  4:86838/1‑147 (MQ=255)
                                                                                     tgtAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGc                                                       >  8:227790/1‑147 (MQ=255)
                                                                                     tgtAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCa                                                      <  1:350494/148‑1 (MQ=255)
                                                                                      gtAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCaa                                                     >  4:53876/1‑148 (MQ=255)
                                                                                       tAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGcaca                                                                                       >  6:380382/1‑113 (MQ=255)
                                                                                       tAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGcaca                                                                                       <  5:380382/113‑1 (MQ=255)
                                                                                        aGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACa                                                   >  2:43231/1‑148 (MQ=255)
                                                                                          aTACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATc                                                 >  7:251448/1‑148 (MQ=255)
                                                                                          aTACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATc                                                 >  4:330243/1‑148 (MQ=255)
                                                                                             ccGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCgg                                               >  4:156117/1‑148 (MQ=255)
                                                                                              cgcgTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGgtg                                             <  7:388961/148‑1 (MQ=255)
                                                                                                cgTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGgtgtg                                           <  6:194809/148‑1 (MQ=255)
                                                                                                cgTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGgtgtg                                           <  3:330243/148‑1 (MQ=255)
                                                                                                 gTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGgtgtg                                           <  3:156117/148‑1 (MQ=255)
                                                                                                 gTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGgtgtgt                                          >  5:370513/1‑148 (MQ=255)
                                                                                                  ttCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGtgtgtg                                         >  7:153056/1‑148 (MQ=255)
                                                                                                  ttCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGtgtgtg                                         >  2:338449/1‑148 (MQ=255)
                                                                                                        aTCCCGGTGGCTGCCATCAAACCTTTGTAATCGGCGCTGAAAAAGTCGGTGGCATCGATCTCGTCCTCGCTCGTGACGTCGAAGGTGCCTCGCACATGGCCGATGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCTTCgg                                   <  8:180748/148‑1 (MQ=255)
                                                                                                        aTCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCgg                                   <  1:348027/148‑1 (MQ=255)
                                                                                                        aTCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCgg                                   <  1:2096/148‑1 (MQ=255)
                                                                                                         tCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGc                                  <  3:86838/148‑1 (MQ=255)
                                                                                                            cGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACt                               >  8:376498/1‑148 (MQ=255)
                                                                                                             gggggCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTTGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACtt                              <  2:433744/148‑1 (MQ=255)
                                                                                                                ggCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACTCGGGCAACATCGGTGTGTGCATCGGCACTTCCg                           <  8:67473/148‑1 (MQ=255)
                                                                                                                  cTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGc                         >  6:204045/1‑148 (MQ=255)
                                                                                                                  cTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCAGc                         <  8:9164/148‑1 (MQ=255)
                                                                                                                   tGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGcc                        >  8:425476/1‑148 (MQ=255)
                                                                                                                    gCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGACGGCTACACCCCCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGcc                        >  2:171650/1‑147 (MQ=255)
                                                                                                                     ccATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGccc                       <  5:57145/147‑1 (MQ=255)
                                                                                                                     ccATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCt                      >  4:355538/1‑148 (MQ=255)
                                                                                                                       aTCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGc                    <  7:250884/148‑1 (MQ=255)
                                                                                                                         aaaCCCGTTGTATTCGGCCCTGAAAAAGTCGGTGTCATCGATCACGTCCTCGCACGTCACTTCGAAGGTGCCTCGCACATGGCCGAGGGCTAAACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCAGGCCCTGcag                  >  3:343617/2‑146 (MQ=255)
                                                                                                                            cccGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCa               >  6:47032/1‑148 (MQ=255)
                                                                                                                            cccGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCa               >  7:269678/1‑148 (MQ=255)
                                                                                                                            cacGTTGTATTCGGCACTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATAGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCa               >  6:87581/3‑148 (MQ=255)
                                                                                                                               gTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCg            <  3:355538/148‑1 (MQ=255)
                                                                                                                                   tATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACAt        <  2:119632/148‑1 (MQ=255)
                                                                                                                                    aTTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCg                  >  5:346517/1‑137 (MQ=255)
                                                                                                                                    aTTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCg                  <  6:346517/137‑1 (MQ=255)
                                                                                                                                    aTTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGCCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATg       >  1:418081/1‑148 (MQ=255)
                                                                                                                                     ttCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATgg      <  5:47032/148‑1 (MQ=255)
                                                                                                                                       cGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGt     >  6:321636/1‑147 (MQ=255)
                                                                                                                                        ggCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTc    <  1:322435/147‑1 (MQ=255)
                                                                                                                                         gCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCAc  <  2:121104/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                              
TATACAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCAC  >  NC_002947_CJ‑RC/4888613‑4888898

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTATACAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCACCGGCCTGTACA  >  NC_002947_CJ‑RC/4888612‑4888909
                                                                                                                                                     |                                                                                                                                                    
GTATACAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAG                                                                                                                                                      <  NB501016:282:H2V7GAFX2:1:11206:16635:10565/150‑1 (MQ=60)
GTATACAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAG                                                                                                                                                      <  NB501016:282:H2V7GAFX2:3:11403:5002:2072/150‑1 (MQ=60)
   TACAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGG                                                                                                                                                 >  NB501016:282:H2V7GAFX2:2:11303:6993:18614/1‑151 (MQ=60)
     CAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTG                                                                                                                                               <  NB501016:282:H2V7GAFX2:4:21610:14444:19951/151‑1 (MQ=60)
     CAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTG                                                                                                                                               >  NB501016:282:H2V7GAFX2:1:11207:14776:10068/1‑151 (MQ=60)
            AAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGA                                                                                                                                        <  NB501016:282:H2V7GAFX2:4:11605:4768:5643/151‑1 (MQ=60)
             AAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTACTGGTCATCCGCCGAGAAGGTGTAGAAACCGCGTTAGGCATCCCGGGGCCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGAT                                                                                                                                       >  NB501016:282:H2V7GAFX2:3:11502:11754:9899/1‑151 (MQ=60)
               ATCGATGAGGCACAATCCTTGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTTAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCA                                                                                                                                     <  NB501016:282:H2V7GAFX2:3:11401:22354:17700/151‑1 (MQ=60)
                         CACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCA                                                                                                                                     <  NB501016:282:H2V7GAFX2:1:21202:12275:8667/141‑1 (MQ=60)
                           CAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCG                                                                                                                         <  NB501016:282:H2V7GAFX2:1:11204:1090:1338/151‑1 (MQ=60)
                           CAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCG                                                                                                                         <  NB501016:282:H2V7GAFX2:1:11205:26814:1333/151‑1 (MQ=60)
                              ACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCC                                                                                                                      <  NB501016:282:H2V7GAFX2:2:11106:3817:15772/151‑1 (MQ=60)
                                    AGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTTGGTGGCATCGATCACGTCCTCGCTCGTCACGTCG                                                                                                                 >  NB501016:282:H2V7GAFX2:1:21307:8675:5545/1‑150 (MQ=60)
                                        AAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGT                                                                                                            <  NB501016:282:H2V7GAFX2:1:21307:4648:17268/151‑1 (MQ=60)
                                                GCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGC                                                                                                     >  NB501016:282:H2V7GAFX2:3:11412:10076:9613/1‑150 (MQ=60)
                                                  AATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACA                                                                                                  <  NB501016:282:H2V7GAFX2:4:11609:7645:4549/151‑1 (MQ=60)
                                                   ATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACAT                                                                                                 >  NB501016:282:H2V7GAFX2:3:21602:6568:18121/1‑151 (MQ=60)
                                                     CGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGG                                                                                               <  NB501016:282:H2V7GAFX2:2:11301:22189:14672/151‑1 (MQ=60)
                                                              CGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCT                                                                                      >  NB501016:282:H2V7GAFX2:2:11207:2118:10335/1‑151 (MQ=60)
                                                                   TGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACC                                                                                 >  NB501016:282:H2V7GAFX2:1:11206:16888:9477/1‑151 (MQ=60)
                                                                    GGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACAC                                                                                  <  NB501016:282:H2V7GAFX2:4:21408:24539:5765/149‑1 (MQ=60)
                                                                       CATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCGCGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCC                                                                             >  NB501016:282:H2V7GAFX2:2:21201:20077:10840/1‑151 (MQ=60)
                                                                         TGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCC                                                                           >  NB501016:282:H2V7GAFX2:4:11608:12329:19643/1‑151 (MQ=60)
                                                                                     GTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAAC                                                               <  NB501016:282:H2V7GAFX2:1:21205:10554:8500/151‑1 (MQ=60)
                                                                                      TGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATG                                                                                                <  NB501016:282:H2V7GAFX2:3:21511:8406:1990/117‑1 (MQ=60)
                                                                                         AGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCG                                                           >  NB501016:282:H2V7GAFX2:4:21402:24867:8504/1‑151 (MQ=60)
                                                                                              CCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTG                                                      <  NB501016:282:H2V7GAFX2:4:21510:2998:8758/151‑1 (MQ=60)
                                                                                                GCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTG                                                    >  NB501016:282:H2V7GAFX2:3:21509:1913:10703/1‑151 (MQ=60)
                                                                                                GCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTG                                                    <  NB501016:282:H2V7GAFX2:2:21205:12874:14333/151‑1 (MQ=60)
                                                                                                 CGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTG                                                    <  NB501016:282:H2V7GAFX2:2:11212:21831:20263/151‑1 (MQ=60)
                                                                                                 CGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGC                                                   >  NB501016:282:H2V7GAFX2:4:11511:12367:6617/1‑151 (MQ=60)
                                                                                                        CATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCA                                            <  NB501016:282:H2V7GAFX2:1:21205:25489:1271/151‑1 (MQ=60)
                                                                                                        CATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCA                                            <  NB501016:282:H2V7GAFX2:1:11101:9555:7016/151‑1 (MQ=60)
                                                                                                         ATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCAC                                           <  NB501016:282:H2V7GAFX2:2:11202:19862:10664/151‑1 (MQ=60)
                                                                                                                     GCCATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTG                                <  NB501016:282:H2V7GAFX2:3:11409:9694:18126/150‑1 (MQ=60)
                                                                                                                       CATCAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCG                             <  NB501016:282:H2V7GAFX2:4:21402:24962:6911/151‑1 (MQ=60)
                                                                                                                        ATAAACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGTCATCGATCACGTCCTCGCACGTCACTTCGAAGGTGCCTCGCACATGGCCGAGGGCTAAACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCAGGCCCTGCAGG                            >  NB501016:282:H2V7GAFX2:2:21207:6967:15383/1‑151 (MQ=60)
                                                                                                                           AACCCGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCAC                         >  NB501016:282:H2V7GAFX2:4:21405:8608:1995/1‑151 (MQ=60)
                                                                                                                               CGTTGTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGAC                     <  NB501016:282:H2V7GAFX2:2:21209:8051:12420/151‑1 (MQ=60)
                                                                                                                                   GTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGC                           >  NB501016:282:H2V7GAFX2:3:21505:14363:14953/1‑141 (MQ=60)
                                                                                                                                   GTATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGCCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGG                 >  NB501016:282:H2V7GAFX2:1:21304:23767:11528/1‑151 (MQ=60)
                                                                                                                                     ATTCGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTC               <  NB501016:282:H2V7GAFX2:3:11408:11631:6692/151‑1 (MQ=60)
                                                                                                                                        CGGCCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCAC             <  NB501016:282:H2V7GAFX2:1:21201:9113:3228/150‑1 (MQ=60)
                                                                                                                                           CCCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCACCGGC         >  NB501016:282:H2V7GAFX2:3:21404:26181:8874/1‑151 (MQ=60)
                                                                                                                                            CCTGAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGACGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGCTCACCGCCC        >  NB501016:282:H2V7GAFX2:1:11212:15440:18025/1‑151 (MQ=60)
                                                                                                                                               tgAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCACCGGCCTG      >  NB501016:282:H2V7GAFX2:4:21512:24260:4883/3‑151 (MQ=60)
                                                                                                                                                  AAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCACCGGCCTGTAC   <  NB501016:282:H2V7GAFX2:4:21603:19602:19062/151‑1 (MQ=60)
                                                                                                                                                   AAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCACCGGCCTGTACA  <  NB501016:282:H2V7GAFX2:2:11107:2816:13886/151‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                    
GTATACAATTACAAAATCGATGAGGCACAAACCATGAGCAAAATGAGAGCAATCGATGCAGCCGTTCTGGTCATGCGCCGTGAAGGTGTAGATACCGCGTTCGGCATCCCGGGGGCTGCCATCAACCCGTTGTATTCGGCCCTGAAAAAAGTCGGTGGCATCGATCACGTCCTCGCTCGTCACGTCGAAGGTGCCTCGCACATGGCCGAGGGCTACACCCGCGCCAACCCGGGCAACATCGGTGTGTGCATCGGCACTTCCGGCCCTGCCGGCACCGACATGGTCACCGGCCTGTACA  >  NC_002947_CJ‑RC/4888612‑4888909

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: