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 F58 I1 R1 60 66.0 2910360 98.3% 2860883 145.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 5,502,776 A→C intergenic (+19/‑116) cobK → / → PP_4833 precorrin‑6x reductase/SURF1 domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC5,502,7760AC98.1% 185.2 / ‑3.2 52intergenic (+19/‑116)cobK/PP_4833precorrin‑6x reductase/SURF1 domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (26/25);  minor base . (1/0);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

CCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACACAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCG  >  NC_002947_CJ‑RC/5502635‑5502919
                                                                                                                                             |                                                                                                                                               
ccGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAgg                                                                                                                                           >  6:167521/1‑148 (MQ=255)
   aaGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCAc                                                                                                                                        >  6:254247/1‑148 (MQ=255)
    aGCTAGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACt                                                                                                                                       <  4:349915/148‑1 (MQ=255)
          aaGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTAc                                                                                                                                  >  4:128788/1‑147 (MQ=255)
             tGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAgg                                                                                                                              <  4:75697/148‑1 (MQ=255)
                cgcgGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGAtt                                                                                                                           >  5:216619/1‑148 (MQ=255)
                  cgGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTc                                                                                                                         <  1:331194/148‑1 (MQ=255)
                  cgGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCa                                                                                                                        <  7:110949/149‑1 (MQ=255)
                     gAGNTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGa                                                                                                                      >  5:222428/1‑148 (MQ=255)
                     gAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCGAAAGGCACTCCTACAAGGATTTCAGa                                                                                                                      >  5:217785/1‑148 (MQ=255)
                                 ccGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTccaa                                                                                                          <  5:254247/148‑1 (MQ=255)
                                  cGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTccaac                                                                                                         >  5:785/1‑148 (MQ=255)
                                  cGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTccaac                                                                                                         <  3:22382/148‑1 (MQ=255)
                                       cTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCaaccaacc                                                                                                    >  4:149681/1‑148 (MQ=255)
                                        ttGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCaaccaacc                                                                                                    >  8:226790/1‑147 (MQ=255)
                                           tCCTGAAACGTGCGACGCTGCCGGTGGTGGTCAGGGCATTCGTCTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGGACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAg                                                                                                 <  6:299678/147‑1 (MQ=255)
                                              tGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGc                                                                                             >  2:154630/1‑148 (MQ=255)
                                               gAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGcc                                                                                            >  3:16204/1‑148 (MQ=255)
                                                 aaCGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCg                                                                                          >  6:81628/1‑148 (MQ=255)
                                                  aCGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGt                                                                                         <  3:311436/148‑1 (MQ=255)
                                                     tCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGcc                                                                                     <  7:226790/149‑1 (MQ=255)
                                                              gCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATctgct                                                                             >  4:148702/1‑148 (MQ=255)
                                                                   aGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCt                                                                         >  5:238672/1‑147 (MQ=255)
                                                                      tGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTaaa                                                                     <  1:233500/148‑1 (MQ=255)
                                                                         aCAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTc                                                                  >  3:186300/1‑147 (MQ=255)
                                                                          cAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTc                                                                  >  5:324046/1‑147 (MQ=255)
                                                                          cAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTc                                                                  >  3:90486/1‑147 (MQ=255)
                                                                            gggCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACg                                                              <  8:287115/149‑1 (MQ=255)
                                                                              gCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGc                                                             >  5:123246/1‑148 (MQ=255)
                                                                                          cTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGa                                                 <  5:166974/148‑1 (MQ=255)
                                                                                             aCGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAgg                                              <  6:337953/148‑1 (MQ=255)
                                                                                                 aTCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGttt                                           >  5:197521/1‑147 (MQ=255)
                                                                                                  tCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCAc                                        >  8:158665/1‑149 (MQ=255)
                                                                                                     ggCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCa                                     >  8:62604/1‑149 (MQ=255)
                                                                                                      gcaaaaGCGATCGCGATTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCATTGGCCCGTTGCCTATCTGCTGCCTTTAACTCGACGCCCCTCCCCCGGAAGGTTTCAc                                        <  2:171996/140‑1 (MQ=255)
                                                                                                      gCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTc                                          <  3:210932/143‑1 (MQ=255)
                                                                                                      gCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCa                                         <  3:147963/144‑1 (MQ=255)
                                                                                                       ccACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTc                                          <  6:253674/142‑1 (MQ=255)
                                                                                                       ccACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTc                                          >  5:253674/1‑142 (MQ=255)
                                                                                                        cacaGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCatat                                  >  7:258449/1‑149 (MQ=255)
                                                                                                          ccgcgcTCGCCCTTTCGGCGCTGCCCTTTCACAACCCAACGGCACTCCTACAATTATTTCAGCCTTTTCTTCAACCCACCCCCTTGTCCCGTGTCCGATCTGCTGCCTTAAACTCGCCGCCCCTCCCCCGGACGGTTTCACCCatattg                                <  8:148163/147‑3 (MQ=255)
                                                                                                               cTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAAt                            <  1:154630/148‑1 (MQ=255)
                                                                                                                  gCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCaa                        >  8:258434/1‑149 (MQ=255)
                                                                                                                  gCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCa                         <  3:128788/148‑1 (MQ=255)
                                                                                                                         aGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAAtggtgg                  <  7:345585/148‑1 (MQ=255)
                                                                                                                          ggCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAAGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGa                 >  5:160209/1‑148 (MQ=255)
                                                                                                                           gcgcTGCCCTTTCACAACCGAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGAt                <  6:222428/148‑1 (MQ=255)
                                                                                                                           gcgcTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGAt                <  6:217785/148‑1 (MQ=255)
                                                                                                                            cgcTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTTCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATct              <  7:282072/149‑1 (MQ=255)
                                                                                                                              cTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATctgg            <  7:158665/149‑1 (MQ=255)
                                                                                                                                    tttCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGt       >  4:187800/1‑148 (MQ=255)
                                                                                                                                     ttCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTc      >  5:238635/1‑148 (MQ=255)
                                                                                                                                        acaacccaaAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCtt    >  6:328238/1‑147 (MQ=255)
                                                                                                                                        acaacccaaAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCtt    <  5:328238/147‑1 (MQ=255)
                                                                                                                                         caacccaaAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCg  <  5:167521/148‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                               
CCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACACAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCG  >  NC_002947_CJ‑RC/5502635‑5502919

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
CGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACACAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGCATTGC  >  NC_002947_CJ‑RC/5502627‑5502926
                                                                                                                                                     |                                                                                                                                                      
CGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCC                                                                                                                                                       >  NB501016:282:H2V7GAFX2:4:11510:18175:9674/1‑151 (MQ=60)
                      GGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTT                                                                                                                                 >  NB501016:282:H2V7GAFX2:3:21407:22123:9028/1‑151 (MQ=60)
                         GCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAG                                                                                                                              <  NB501016:282:H2V7GAFX2:4:11510:13187:1165/151‑1 (MQ=60)
                         GCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAG                                                                                                                              <  NB501016:282:H2V7GAFX2:1:21306:25144:2371/151‑1 (MQ=60)
                           GGGAGNTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGAC                                                                                                                            >  NB501016:282:H2V7GAFX2:3:21408:1063:13549/1‑151 (MQ=60)
                           GGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCGAAAGGCACTCCTACAAGGATTTCAGAC                                                                                                                            >  NB501016:282:H2V7GAFX2:3:21407:26604:13548/1‑151 (MQ=60)
                                        GCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACC                                                                                                               <  NB501016:282:H2V7GAFX2:3:21503:24862:4199/151‑1 (MQ=60)
                                        GCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACC                                                                                                               >  NB501016:282:H2V7GAFX2:3:11401:17524:4111/1‑151 (MQ=60)
                                         CCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCA                                                                                                              <  NB501016:282:H2V7GAFX2:2:11105:2046:11522/151‑1 (MQ=60)
                                               CTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCC                                                                                                          >  NB501016:282:H2V7GAFX2:3:21603:19108:10854/1‑149 (MQ=60)
                                                     CTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCC                                                                                                  >  NB501016:282:H2V7GAFX2:2:11104:7175:7081/1‑151 (MQ=60)
                                                         AACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTG                                                                                              <  NB501016:282:H2V7GAFX2:2:21305:3806:7949/151‑1 (MQ=60)
                                                            GTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCG                                                                                           <  NB501016:282:H2V7GAFX2:4:21408:20312:13567/151‑1 (MQ=60)
                                                                         GGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTT                                                                               >  NB501016:282:H2V7GAFX2:3:21412:18021:2272/1‑150 (MQ=60)
                                                                             GTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACT                                                                          <  NB501016:282:H2V7GAFX2:1:21109:6972:16728/151‑1 (MQ=60)
                                                                               GGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCG                                                                        >  NB501016:282:H2V7GAFX2:2:21101:17450:8885/1‑150 (MQ=60)
                                                                                GACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCG                                                                        >  NB501016:282:H2V7GAFX2:3:21604:12888:15796/1‑150 (MQ=60)
                                                                                GACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCG                                                                        >  NB501016:282:H2V7GAFX2:2:11207:23517:3299/1‑150 (MQ=60)
                                                                                    GGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCC                                                                   >  NB501016:282:H2V7GAFX2:3:11601:3403:3959/1‑151 (MQ=60)
                                                                                                 ACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAG                                                      <  NB501016:282:H2V7GAFX2:3:11609:17221:5783/151‑1 (MQ=60)
                                                                                                       GGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGATGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTC                                                 >  NB501016:282:H2V7GAFX2:3:21403:5278:7283/1‑150 (MQ=60)
                                                                                                             GGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACC                                              <  NB501016:282:H2V7GAFX2:2:11306:7060:3728/147‑1 (MQ=60)
                                                                                                             GGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCAC                                               <  NB501016:282:H2V7GAFX2:2:21106:1746:14182/146‑1 (MQ=60)
                                                                                                             GGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCAC                                               >  NB501016:282:H2V7GAFX2:3:21503:10206:2009/1‑146 (MQ=60)
                                                                                                               CCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATG                                        >  NB501016:282:H2V7GAFX2:4:21502:9627:13774/1‑151 (MQ=60)
                                                                                                                      GCTCGCCCTTTAGGCGCTGCCATTTCACAACCCAAAGGCACTCATACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGC                                 <  NB501016:282:H2V7GAFX2:1:11306:12798:4508/151‑1 (MQ=60)
                                                                                                                         CGCCCTTTAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAAT                              <  NB501016:282:H2V7GAFX2:2:11302:9502:11233/151‑1 (MQ=60)
                                                                                                                                TAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGA                        <  NB501016:282:H2V7GAFX2:4:21607:16621:12168/150‑1 (MQ=60)
                                                                                                                                TAGGCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAAGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGAT                       >  NB501016:282:H2V7GAFX2:3:11607:22254:20142/1‑151 (MQ=60)
                                                                                                                                   GCGCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTTCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTG                    <  NB501016:282:H2V7GAFX2:4:21507:4365:5471/151‑1 (MQ=60)
                                                                                                                                     GCTGCCCTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGC                  <  NB501016:282:H2V7GAFX2:4:11607:9641:8456/151‑1 (MQ=60)
                                                                                                                                           CTTTCACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCT            >  NB501016:282:H2V7GAFX2:3:21412:23397:2105/1‑151 (MQ=60)
                                                                                                                                               CACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCG         <  NB501016:282:H2V7GAFX2:3:21605:24719:11703/150‑1 (MQ=60)
                                                                                                                                                ACAACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGC       <  NB501016:282:H2V7GAFX2:3:11609:21736:7770/151‑1 (MQ=60)
                                                                                                                                                  AACCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGCAT     >  NB501016:282:H2V7GAFX2:2:11312:4904:18006/1‑151 (MQ=60)
                                                                                                                                                    CCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGCATTG   >  NB501016:282:H2V7GAFX2:1:21212:2498:3327/1‑151 (MQ=60)
                                                                                                                                                    CCCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGCATTG   >  NB501016:282:H2V7GAFX2:1:11204:11265:14308/1‑151 (MQ=60)
                                                                                                                                                     CCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGCATTGG  >  NB501016:282:H2V7GAFX2:1:11312:21301:11301/1‑151 (MQ=60)
                                                                                                                                                     CCAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGCATTGC  <  NB501016:282:H2V7GAFX2:1:21111:3876:10275/151‑1 (MQ=60)
                                                                                                                                                     |                                                                                                                                                      
CGACCGAGCCGAAGCTGGAAGTGGCGCGGGAGCTTGGCGTGCCGGTGCTTGTCCTGAAACGTCCGACGCTGCCGGAGGTGGACAGGGCATTCGACTCACTGACGGATCTGGCCACAGCGCTCGCCCTTTAGGCGCTGCCCTTTCACAACACAAAGGCACTCCTACAAGGATTTCAGACTTTTCTTCCAACCAACCCAGTGGCCCGTTGCCGATCTGCTGCCTTAAACTCGACGCCCCTCCCCCGGAAGGTTTCACCCATATGGAAATGCAATGGTGGATCTGGCTGGTCTTCGGCATTGC  >  NC_002947_CJ‑RC/5502627‑5502926

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