Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A2 F12 I1 R1 58 76.7 3263832 98.3% 3208346 146.5

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 4,768,005 T→A intergenic (+112/‑70) fpvA → / → tpaK TonB‑dependent outer membrane ferripyoverdine receptor FpvA/probable terephthalate transporter, MFS superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,768,0050TA75.5% 131.9 / 45.4 53intergenic (+112/‑70)fpvA/tpaKTonB‑dependent outer membrane ferripyoverdine receptor FpvA/probable terephthalate transporter, MFS superfamily protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (14/26);  minor base . (7/6);  total (21/32)
Fisher's exact test for biased strand distribution p-value = 3.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

CAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAC  >  NC_002947_CJ‑RC/4767860‑4768140
                                                                                                                                                 |                                                                                                                                       
cAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATaa                                                                                                                                      <  5:225399/149‑1 (MQ=255)
  aCTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAt                                                                                                                                     <  8:35323/148‑1 (MQ=255)
           tCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGAc                                                                                                                           <  4:131619/149‑1 (MQ=255)
               cGCCAAGTACAGCATCAGATAGCCTGACGATGAACGCAAAAAACCTCACCCAGGTGCGGTTTTTTGAATTCACTATTAGTCACAATGGACGTAACGTATGAATACAAAAGACAATTAATCATCGGCAAGTATAATGTTATCAGACCt                                                                                                                         >  8:178120/1‑147 (MQ=255)
                 ccAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGa                                                                                                                      >  7:243405/1‑148 (MQ=255)
                     gTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGATTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATGACAAATAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                 >  2:324234/1‑149 (MQ=255)
                      tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                 <  8:106007/148‑1 (MQ=255)
                      tACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATtgt                                                                                                                 <  7:409571/148‑1 (MQ=255)
                        cAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGa                                                                                                               <  4:32932/148‑1 (MQ=255)
                         aggttCTGTTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGa                                                                                                               <  8:120810/144‑1 (MQ=21)
                          gCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGc                                                                                                             <  8:388253/148‑1 (MQ=255)
                                 aTCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATaa                                                                                                       <  1:143948/147‑1 (MQ=255)
                                    gCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAAtt                                                                                                  <  1:303987/149‑1 (MQ=255)
                                      cTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCt                                                                                                >  4:104223/1‑149 (MQ=255)
                                         aCGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                             <  3:90639/149‑1 (MQ=255)
                                          cGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTa                                                                                              <  7:381644/147‑1 (MQ=255)
                                           gTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCtt                                                                                               >  4:128054/1‑145 (MQ=255)
                                              gAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGAt                                                                                          >  7:306168/1‑147 (MQ=255)
                                                 cGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTaa                                                                                       <  7:178120/147‑1 (MQ=255)
                                                 cGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCa                                                                     >  6:318824/1‑149 (MQ=12)
                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                             >  1:379498/1‑140 (MQ=255)
                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTaa                                                                                             <  2:379498/140‑1 (MQ=255)
                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATc                                                                      >  7:253636/1‑147 (MQ=12)
                                                  gCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCat                                                                    <  1:206889/149‑1 (MQ=12)
                                                    aaaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTca                                                                                   >  4:88209/1‑148 (MQ=255)
                                                    aaaaaaCCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATAATCGGCTCGTATAATGTGATCAGACATGGAATTGTGAGCGGATAAAAATTCTTAAGATTAACTca                                                                                   >  4:272024/1‑148 (MQ=17)
                                                          ccGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCt                                                             >  1:102649/1‑148 (MQ=12)
                                                          ccGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCt                                                             >  1:35775/1‑148 (MQ=12)
                                                             cACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTa                                                           >  8:159465/1‑147 (MQ=12)
                                                             cACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAg                                                          <  8:193119/148‑1 (MQ=12)
                                                                ccAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTACGGAGGTCAGGTATGATTACGATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATTTCATATGGCCTTAGCAc                                                       >  4:320574/1‑148 (MQ=9)
                                                                   ggTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATc                                                   <  1:33563/149‑1 (MQ=12)
                                                                    gTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCatat                                                                  >  3:207447/1‑149 (MQ=255)
                                                                     tGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATg                                                                 >  4:62624/1‑149 (MQ=255)
                                                                                gAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACt                                       <  5:278129/148‑1 (MQ=12)
                                                                                  aTTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGcc                                    <  5:271518/149‑1 (MQ=12)
                                                                                    tCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATc                                                   <  8:222664/148‑1 (MQ=255)
                                                                                    tCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCATTAGCACCATCGCGCGTGACACTGCCg                                   >  5:90328/1‑148 (MQ=9)
                                                                                      aCTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATcgc                                                 <  4:276861/148‑1 (MQ=255)
                                                                                      aCTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGT‑‑‑‑‑‑‑‑AT‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGAt                                <  2:136241/149‑1 (MQ=12)
                                                                                          gtagtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGt                                             >  6:359795/1‑148 (MQ=255)
                                                                                             gtCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGaca                                         >  2:229858/1‑149 (MQ=255)
                                                                                                  tATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTg                                      <  5:406250/147‑1 (MQ=255)
                                                                                                       aGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGAtt                               <  5:126557/149‑1 (MQ=255)
                                                                                                          tCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCATTAGCACCATCGCGCGTGACACTGCCGGATTTCa                            <  2:305256/149‑1 (MQ=255)
                                                                                                                    ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgc                   <  2:156811/148‑1 (MQ=255)
                                                                                                                    ttACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTcgc                   <  3:62624/148‑1 (MQ=255)
                                                                                                                      aCTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCgg                 <  4:113890/148‑1 (MQ=255)
                                                                                                                          ttgaCAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGTGGCCAGt            <  6:57516/146‑1 (MQ=255)
                                                                                                                               aaTTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGt        <  8:243405/148‑1 (MQ=255)
                                                                                                                               aaTTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTc       <  4:207447/149‑1 (MQ=255)
                                                                                                                                    atcatcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAc  >  2:388762/1‑149 (MQ=255)
                                                                                                                                     tcatcGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAc  <  7:362588/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                       
CAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATAC  >  NC_002947_CJ‑RC/4767860‑4768140

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
GCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACATCGTG  >  NC_002947_CJ‑RC/4767859‑4768152
                                                                                                                                                  |                                                                                                                                                   
GCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAAT                                                                                                                                                 <  NB501016:282:H2V7GAFX2:3:21404:26727:4901/151‑1 (MQ=60)
                 GCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAAT                                                                                                                                >  NB501016:282:H2V7GAFX2:4:21406:2756:9503/1‑151 (MQ=60)
                     AGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTG                                                                                                                            <  NB501016:282:H2V7GAFX2:4:21611:21634:5079/151‑1 (MQ=60)
                         CAGCTTCTGATCGCCTTACGTTGAACGAAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAG                                                                                                        <  NB501016:282:H2V7GAFX2:2:11307:3269:17714/151‑1 (MQ=56)
                         CAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAG                                                                                                        <  NB501016:282:H2V7GAFX2:3:21508:1326:10264/151‑1 (MQ=60)
                              TCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATT                                                                                                     <  NB501016:282:H2V7GAFX2:1:21203:16314:12366/149‑1 (MQ=60)
                                 GATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAAC                                                                                                                  <  NB501016:282:H2V7GAFX2:1:11301:12920:4463/149‑1 (MQ=60)
                                    CGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTC                                                                                                             <  NB501016:282:H2V7GAFX2:1:21205:8796:11927/151‑1 (MQ=60)
                                    CGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACA                                                                                             <  NB501016:282:H2V7GAFX2:1:21202:8048:12793/151‑1 (MQ=57)
                                         GACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                                                                        <  NB501016:282:H2V7GAFX2:2:11205:5300:4642/151‑1 (MQ=60)
                                         GACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAA                                                                                                         <  NB501016:282:H2V7GAFX2:4:21606:7238:4750/150‑1 (MQ=60)
                                            GTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGAT                                                                                     <  NB501016:282:H2V7GAFX2:4:21409:2142:10963/151‑1 (MQ=45)
                                              TGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTA                                                                                                    >  NB501016:282:H2V7GAFX2:4:21505:11291:1791/1‑150 (MQ=52)
                                                AACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCATGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAAC                                                                                                  <  NB501016:282:H2V7GAFX2:4:11607:16581:2776/150‑1 (MQ=46)
                                                  CGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAG                                                                                                        >  NB501016:282:H2V7GAFX2:1:21308:21985:3282/1‑142 (MQ=45)
                                                  CGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATA                                                                               <  NB501016:282:H2V7GAFX2:1:11312:12697:12187/151‑1 (MQ=35)
                                                  CGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCAT                                                                                >  NB501016:282:H2V7GAFX2:4:21408:7060:4257/1‑150 (MQ=35)
                                                          ACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTT                                                                        >  NB501016:282:H2V7GAFX2:1:11106:7243:14561/1‑150 (MQ=30)
                                                          ACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTT                                                                        >  NB501016:282:H2V7GAFX2:1:11205:3258:18847/1‑150 (MQ=30)
                                                                   AGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCG                                                              <  NB501016:282:H2V7GAFX2:1:11106:16643:7971/151‑1 (MQ=32)
                                                                    GGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATT                                                                             >  NB501016:282:H2V7GAFX2:2:11312:20809:14394/1‑151 (MQ=5)
                                                                                TGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTG                                                  <  NB501016:282:H2V7GAFX2:3:21502:12796:3286/150‑1 (MQ=30)
                                                                                  AATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCG                                               <  NB501016:282:H2V7GAFX2:3:21412:11342:18858/151‑1 (MQ=32)
                                                                                    TTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCATTAGCACCATCGCGCGTGACACTGCCGG                                              >  NB501016:282:H2V7GAFX2:3:11505:7108:4391/1‑150 (MQ=22)
                                                                                              GTCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGAT                                   >  NB501016:282:H2V7GAFX2:4:11402:4790:4954/1‑151 (MQ=43)
                                                                                               TCTCTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATT                                  >  NB501016:282:H2V7GAFX2:4:11510:15215:11948/1‑151 (MQ=45)
                                                                                                  CTATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGC                                                 <  NB501016:282:H2V7GAFX2:3:21612:14001:15358/149‑1 (MQ=26)
                                                                                                    ATGGAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGCGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG                              <  NB501016:282:H2V7GAFX2:2:11201:25744:20208/150‑1 (MQ=50)
                                                                                                       GAGGTCAGGTATGATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTT                                          <  NB501016:282:H2V7GAFX2:3:11511:7225:8626/151‑1 (MQ=40)
                                                                                                              GGTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTC                   <  NB501016:282:H2V7GAFX2:4:21602:20814:9925/151‑1 (MQ=60)
                                                                                                               GTATGATTACTATTGACAATTAATCATCGGCTCGTAT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAATTGTGAGCGGATAAGAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCG                  <  NB501016:282:H2V7GAFX2:2:21201:5719:15219/151‑1 (MQ=60)
                                                                                                                    ATTACTATTGACAATTAATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCG                              <  NB501016:282:H2V7GAFX2:2:11112:6246:5504/150‑1 (MQ=60)
                                                                                                                                    AATCATCGGCTCGTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACC             <  NB501016:282:H2V7GAFX2:4:21602:21092:15308/151‑1 (MQ=60)
                                                                                                                                                GTATAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCTCGGCCAGTGAGTCGATACCAGTACATCGTG  >  NB501016:282:H2V7GAFX2:3:11610:14874:3614/1‑150 (MQ=60)
                                                                                                                                                  |                                                                                                                                                   
GCAACTTCACCATCACCGCCAAGTACAGCTTCTGATCGCCTGACGTTGAACGCAAAAAACCGCACCCAGGTGCGGTTTTTTGAATTCACTAGTAGTCTCTATGGAGGTCAGGTATGATTACTAATTACAATTAATCATCGGCTCGTTTAATGTGATCAGACCTGGAATTGTGAGCGGATAACAATTCTTAAGATTAACTCACACAGGAGATATCATATGTCCTTAGCACCATCGCGCGTGACACTGCCGGATTTCATCGATTCGCGGCCAGTGAGTCGATACCAGTACATCGTG  >  NC_002947_CJ‑RC/4767859‑4768152

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: