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 F48 I1 R1 58 71.3 3054194 98.3% 3002272 146.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 336,125 +T coding (210/237 nt) PP_0278 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC336,1241.T96.1% 166.2 / 0.7 51coding (209/237 nt)PP_0278hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (23/26);  total (24/27)
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.78e-01

ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATT  >  NC_002947_CJ‑RC/335980‑336267
                                                                                                                                                  |                                                                                                                                                
aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAttac                                                                                                                                              <  7:126910/149‑4 (MQ=255)
aCCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGAttac                                                                                                                                              <  8:356068/149‑4 (MQ=255)
     tCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCCGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCtgtg                                                                                                                                         <  6:28871/149‑1 (MQ=255)
      cccTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGc                                                                                                                                        <  4:153564/149‑1 (MQ=255)
        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGa                                                                                                                                      <  2:58887/149‑1 (MQ=255)
        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGa                                                                                                                                      >  8:181315/1‑149 (MQ=255)
        cTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGa                                                                                                                                      >  7:207323/1‑149 (MQ=255)
                 cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCa                                                                                                                              <  7:143336/148‑1 (MQ=255)
                 cTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTAGATCCTGATGGGCTACATCAAGCAACTGGTGCGCGAGGGCGATATGGAGGACGACCATAAAATACTGGGTGCAGCCGATAACCTGTGCGACATGGca                                                                                                                               >  1:144343/1‑146 (MQ=255)
                    cgcgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAgg                                                                                                                           <  7:273944/148‑1 (MQ=255)
                      cgTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGgcgc                                                                                                                        <  3:334770/149‑1 (MQ=255)
                       gTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCAAAGGCGCt                                                                                                                       >  1:274342/1‑149 (MQ=255)
                         agagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtga                                                                                                                     >  5:34088/1‑149 (MQ=255)
                         agagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtga                                                                                                                     >  5:64547/1‑149 (MQ=255)
                          gagCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCtgat                                                                                                                    <  3:222369/149‑1 (MQ=255)
                             cccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCCAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGAGGaa                                                                                                                 >  1:336381/1‑149 (MQ=255)
                             cccGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGaa                                                                                                                 >  1:325034/1‑149 (MQ=255)
                                    tCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                          <  6:11717/149‑1 (MQ=255)
                                      cccGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATc                                                                                                          >  5:23792/1‑147 (MQ=255)
                                        cGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCACCCGATTACCTGTCAGACATGCCCAAGACGCTAATAAACAATATCaag                                                                                                       >  4:259556/1‑148 (MQ=255)
                                              cgcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGcc                                                                                                 >  8:385414/1‑148 (MQ=255)
                                               gcgGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGcc                                                                                                 >  4:185399/1‑147 (MQ=255)
                                                    cTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCga                                                                                                        >  5:375398/1‑135 (MQ=255)
                                                            aTGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCAt                                                                                  >  3:23913/1‑149 (MQ=255)
                                                               tctcGATCCTGATGGGCTACATCAAGCACCTGGTTCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGTGCATCCTATTACCTGTGCGACATGGCAAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTa                                                                                <  7:37579/148‑1 (MQ=255)
                                                               tctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAg                                                                               >  7:112449/1‑149 (MQ=255)
                                                                ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                              >  3:255210/1‑149 (MQ=255)
                                                                ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                              >  1:121009/1‑149 (MQ=255)
                                                                ctcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                              >  5:278736/1‑149 (MQ=255)
                                                                 tcGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAgg                                                                              >  7:9993/1‑148 (MQ=255)
                                                                   gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGc                                                                             <  5:134163/147‑1 (MQ=255)
                                                                   gATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCa                                                                            <  7:385414/148‑1 (MQ=255)
                                                                    aTCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAtt                                                                          >  3:311826/1‑149 (MQ=255)
                                                                        tGATTGGTTACATCAATCAACTGGTGCGGGAGGGCGATATGGTGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAc                                                                       <  3:259556/148‑1 (MQ=255)
                                                                        tGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGAcg                                                                      >  7:379402/1‑149 (MQ=255)
                                                                               cTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGCGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCAt                                                               >  1:303860/1‑149 (MQ=255)
                                                                                         cACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGc                                                     <  6:215620/149‑1 (MQ=255)
                                                                                         cACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGc                                                     <  2:124118/149‑1 (MQ=255)
                                                                                          aCCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCg                                                    >  7:297491/1‑149 (MQ=255)
                                                                                           ccTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCGAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCgg                                                   <  4:226910/149‑1 (MQ=255)
                                                                                                   gcgAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGCCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCGTGCGt                                           <  8:144751/149‑1 (MQ=255)
                                                                                                         gCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCGAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGaa                                     <  4:109219/149‑1 (MQ=255)
                                                                                                         gCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGaa                                     <  7:181315/149‑1 (MQ=255)
                                                                                                             tatGGAGGATGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATgg                                  <  2:332272/148‑1 (MQ=255)
                                                                                                             tatGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCAAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATggg                                 <  2:274342/149‑1 (MQ=255)
                                                                                                               tGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCt                               <  2:119236/149‑1 (MQ=255)
                                                                                                                       gacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCCAGCGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagca                       <  2:373467/149‑1 (MQ=255)
                                                                                                                       gacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagca                       <  4:23913/149‑1 (MQ=255)
                                                                                                                       gacCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagca                       <  8:388436/149‑1 (MQ=255)
                                                                                                                       gacCAATAATTCCTGCGCCCAGCCGCTTACCTGTGCGACATGGCCTAGGCGCTGATGAACTATATCGAGATCGCCTAGGGTAAGGCTCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAagca                       <  2:365155/149‑1 (MQ=255)
                                                                                                                                      ggggCAGCCTATTACCGGTGCCACAGGCCCACGGCGCTGAGGAACGAAATCGAGATCCCCACGGGTACGGCGCATTAGCCATTGCCGCGGTCATTTGGGGGCGCGGCCTTGCGGCGCTAATGGGCTGCAAAGCAGCCCCGGCGATTTCa        <  8:298479/149‑1 (MQ=255)
                                                                                                                                        ggCAGCCTATTACCTGTGCCACATGGCCAAGGCGCTTATCAACGATATCGAGATCGCCAAGGGTAAGGCCCATTAGGCATTGCCGCGGTCATTTGGGGGAGCGGCCTTGCGTCGCTAATGGGCTGAAAAGCCGCCCCGGCGATTTCaga      <  8:379402/149‑1 (MQ=255)
                                                                                                                                            gCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGAtt  >  1:219461/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATT  >  NC_002947_CJ‑RC/335980‑336267

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
AACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTG  >  NC_002947_CJ‑RC/335979‑336268
                                                                                                                                                   |                                                                                                                                                
AACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACC                                                                                                                                              <  NB501016:282:H2V7GAFX2:4:11511:15740:12943/151‑1 (MQ=60)
        CCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGAC                                                                                                                                      >  NB501016:282:H2V7GAFX2:4:21402:20691:15262/1‑151 (MQ=60)
                 CCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAA                                                                                                                              <  NB501016:282:H2V7GAFX2:4:11602:24812:14512/150‑1 (MQ=60)
                 CCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTAGATCCTGATGGGCTACATCAAGCAACTGGTGCGCGAGGGCGATATGGAGGACGACCATAAAATACTGGGTGCAGCCGATAACCTGTGCGACATGGCAA                                                                                                                               >  NB501016:282:H2V7GAFX2:1:11302:16939:17509/1‑149 (MQ=60)
                    ACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGC                                                                                                                           <  NB501016:282:H2V7GAFX2:4:21502:1967:12342/150‑1 (MQ=60)
                      GCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCT                                                                                                                        <  NB501016:282:H2V7GAFX2:2:21306:9623:12165/151‑1 (MQ=60)
                       CGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTG                                                                                                                       >  NB501016:282:H2V7GAFX2:4:21404:21297:2883/1‑151 (MQ=60)
                       CGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCAAAGGCGCTG                                                                                                                       >  NB501016:282:H2V7GAFX2:1:21203:3185:10054/1‑151 (MQ=60)
                         TAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGAT                                                                                                                     >  NB501016:282:H2V7GAFX2:3:11501:3156:14578/1‑151 (MQ=60)
                         TAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGAT                                                                                                                     >  NB501016:282:H2V7GAFX2:3:11407:5297:14462/1‑151 (MQ=60)
                          CGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATG                                                                                                                    <  NB501016:282:H2V7GAFX2:2:21106:13875:12110/151‑1 (MQ=60)
                             GCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCCAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGAGGAAC                                                                                                                 >  NB501016:282:H2V7GAFX2:1:21304:24324:2070/1‑151 (MQ=60)
                             GCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAAC                                                                                                                 >  NB501016:282:H2V7GAFX2:1:21301:11319:3950/1‑151 (MQ=60)
                                CGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGAT                                                                                                              >  NB501016:282:H2V7GAFX2:2:21109:3854:7432/1‑151 (MQ=60)
                                      CCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCG                                                                                                          >  NB501016:282:H2V7GAFX2:3:11405:23078:13696/1‑149 (MQ=60)
                                                    CCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAG                                                                                                        >  NB501016:282:H2V7GAFX2:3:21611:24912:17031/1‑137 (MQ=60)
                                                            CATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATT                                                                                  >  NB501016:282:H2V7GAFX2:2:11105:20704:13262/1‑151 (MQ=60)
                                                               GTCTCGATCCTGATGGGCTACATCAAGCACCTGGTTCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGTGCATCCTATTACCTGTGCGACATGGCAAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAG                                                                                <  NB501016:282:H2V7GAFX2:4:11407:13990:13542/150‑1 (MQ=60)
                                                               GTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGG                                                                               >  NB501016:282:H2V7GAFX2:4:11509:18657:1759/1‑151 (MQ=60)
                                                                TCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGC                                                                              >  NB501016:282:H2V7GAFX2:1:11210:8051:11699/1‑151 (MQ=60)
                                                                TCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGC                                                                              >  NB501016:282:H2V7GAFX2:3:21505:18033:17768/1‑151 (MQ=60)
                                                                TCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGC                                                                              >  NB501016:282:H2V7GAFX2:2:21201:10741:12480/1‑151 (MQ=60)
                                                                 CTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGC                                                                              >  NB501016:282:H2V7GAFX2:4:11402:3823:16297/1‑150 (MQ=60)
                                                                   CGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCAT                                                                            <  NB501016:282:H2V7GAFX2:4:21611:9100:7413/150‑1 (MQ=60)
                                                                   CGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCA                                                                             <  NB501016:282:H2V7GAFX2:3:11602:23132:5932/149‑1 (MQ=60)
                                                                   CGATCATGATGGTCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGTGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATT                                                                           <  NB501016:282:H2V7GAFX2:1:11209:17494:16559/151‑1 (MQ=60)
                                                                    GATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTG                                                                          >  NB501016:282:H2V7GAFX2:2:21302:14200:3940/1‑151 (MQ=60)
                                                                        ttgattggtTACATCAATCAACTGGTGCGGGAGGGCGATATGGTGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACG                                                                       <  NB501016:282:H2V7GAFX2:2:21203:3461:5271/141‑1 (MQ=60)
                                                                        CTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGC                                                                      >  NB501016:282:H2V7GAFX2:4:21610:3498:4553/1‑151 (MQ=60)
                                                                               GCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGCGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATT                                                               >  NB501016:282:H2V7GAFX2:1:21209:6970:3347/1‑151 (MQ=60)
                                                                                          CACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGG                                                    >  NB501016:282:H2V7GAFX2:4:21506:22083:17166/1‑151 (MQ=60)
                                                                                                         GGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAAT                                     <  NB501016:282:H2V7GAFX2:4:11609:14435:18841/151‑1 (MQ=60)
                                                                                                                            ATAAATTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGggggggggggggggggggggggggggggacgacc                  >  NB501016:282:H2V7GAFX2:3:11507:6328:20248/1‑117 (MQ=60)
                                                                                                                                 TTCCTGGGGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGggggggggggggggggggggggggcgggacgccccgcgc             >  NB501016:282:H2V7GAFX2:1:21206:11612:6399/1‑112 (MQ=60)
                                                                                                                                        GGGCAGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGggggggggggggggcagggggggggggacggccccggcccgggga       >  NB501016:282:H2V7GAFX2:1:21202:25885:15216/1‑105 (MQ=60)
                                                                                                                                            AGCCGATTACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTG  >  NB501016:282:H2V7GAFX2:1:21104:24746:19713/1‑151 (MQ=60)
                                                                                                                                                   |                                                                                                                                                
AACCCATCCCTCAATCTCCTACGCGTAGAGCCCGGCGTCCCCGTCGACGCGGCCTACGAGCATGTCTCGATCCTGATGGGCTACATCAAGCACCTGGTGCGCGAGGGCGATATGGAGGACGACCATAAATTCCTGGGGGCAGCCGATACCTGTGCGACATGGCCAAGGCGCTGATGAACGATATCGAGATCGCCAAGGGTAAGGCGCATTAGGCATTGACGCGGTCATTTGTGGGAGCGGCCTTGCGTCGCGAATGGGCTGCAAAGCAGCCCCGGCGATTTCAGAGATTG  >  NC_002947_CJ‑RC/335979‑336268

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