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 F38 I1 R1 54 64.3 2734350 98.3% 2687866 146.9

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 5,677,566 +G intergenic (+113/+27) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC5,677,5651.G98.0% 165.2 / ‑3.1 49intergenic (+112/+28)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (26/22);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

GAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAA  >  NC_002947_CJ‑RC/5677421‑5677703
                                                                                                                                                 |                                                                                                                                           
gAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTATTGATCTCTACAGGAGCGGCCTTGTGTCCTGATAGGGCCGCAAAgcggcc                                                                                                                                         >  1:99649/1‑146 (MQ=255)
 aGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcggcc                                                                                                                                         <  6:293822/148‑4 (MQ=255)
    gcCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAgcggccccct                                                                                                                                     <  8:220503/149‑8 (MQ=255)
        aGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGt                                                                                                                                 <  3:124844/149‑1 (MQ=255)
             tgctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGtttt                                                                                                                              >  1:151122/1‑147 (MQ=255)
               ctgACCCGCAAGGGCTGATAGAAGTACCGAGTGGGCTGCTTCGCAGCCCTATCACGACAAAAGGCAGCTCCTGCAGAGAATTTCTTGATCTCGACAGGAGCGGCCTTGTGTCGTGATAGTGCCGCAAAGCGGCCCCCTCGGTTTTCAGc                                                                                                                          >  1:35582/1‑149 (MQ=38)
               ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGc                                                                                                                          <  8:246933/149‑1 (MQ=255)
               ctgACCCGCAAGGCCTGATACAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACCCAAGGCCGCTCCTGCAGAGAATTTCTTCATCTCTCCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAg                                                                                                                           <  7:285307/148‑1 (MQ=25)
                  aCCCTCACGGGCTCATAGAAGTACCCCGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACCGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTg                                                                                                                        <  4:192448/148‑1 (MQ=25)
                    ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTACTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTGAGCTgg                                                                                                                       >  3:15434/1‑147 (MQ=37)
                     cGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTGTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTGAGCTgg                                                                                                                       >  3:42227/1‑146 (MQ=37)
                          gggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAAtt                                                                                                               <  8:79743/149‑1 (MQ=255)
                                    acgtcccGAGCGGGCTGCCTCCCAGCCCTATCACCCCACACCGCCGCCCCTGCAGAGAATTACTTCACCTCTCCAGCAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                                                                                                     <  4:15434/144‑1 (MQ=25)
                                      gTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTca                                                                                                   <  7:341976/149‑1 (MQ=255)
                                      gTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTCAGCTGGCAGAAATTACTTGTGCTTca                                                                                                   >  2:11559/1‑149 (MQ=37)
                                             gcgggCTGCTGCGCCGCCCTCCCCCGCCACAAGGCCCCGCCTCCACCCAATTGCTTGCTCCCAACCGGCGCGGCGTTGTGTCGGGATAGGGCCGCAAAGCGCCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                            <  4:175241/147‑1 (MQ=255)
                                               gggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCGCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                          <  8:296368/149‑1 (MQ=255)
                                                     ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGc                                                                                                       >  8:112040/1‑130 (MQ=255)
                                                     ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGc                                                                                                       <  7:112040/130‑1 (MQ=255)
                                                                 aTCACGACACAAGGCCGCTCCTGTAGAGATCAAGAAATTCTCTGCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtc                                                                          >  2:48388/1‑147 (MQ=14)
                                                                   cACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                      <  4:206171/149‑1 (MQ=255)
                                                                    aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                      <  6:289815/148‑1 (MQ=255)
                                                                       acacAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCt                                                                  >  1:96816/1‑149 (MQ=255)
                                                                         acacGGCCGCTCTTGCTGACAATTTCTTTACCTCTACTGTAGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                <  5:183738/149‑1 (MQ=255)
                                                                         acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTAGTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                >  2:285081/1‑149 (MQ=255)
                                                                         acaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                <  7:211198/149‑1 (MQ=255)
                                                                           aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCg                                                               >  2:153080/1‑148 (MQ=255)
                                                                             ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                            >  5:238146/1‑149 (MQ=255)
                                                                             ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                            >  4:161019/1‑149 (MQ=255)
                                                                             ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                            <  2:151122/149‑1 (MQ=255)
                                                                              gCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                           <  6:238146/149‑1 (MQ=255)
                                                                                cgatCCTGCAGATAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAg                                                         <  3:133450/146‑1 (MQ=39)
                                                                                     cTGCAGAGAATTTCTTGATCTCTACAGGCGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCAc                                                      <  4:233169/147‑1 (MQ=255)
                                                                                       gCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                  >  3:93666/1‑149 (MQ=255)
                                                                                         agagAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGc                                                >  2:247729/1‑149 (MQ=255)
                                                                                                        ctctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCAGAAATTACTTGTGGTTCACACCTACGACACTGATATCCAGGTCCAGGGTCTGCGAAGTTGCACCAGGGCCGTTGATGGCGAAGt                                 >  4:64837/1‑149 (MQ=37)
                                                                                                              aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttca                           <  1:92627/149‑1 (MQ=255)
                                                                                                                gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                         >  3:115540/1‑149 (MQ=255)
                                                                                                                gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCACCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                         <  3:57991/149‑1 (MQ=255)
                                                                                                                       cTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCg                  <  1:48388/149‑1 (MQ=255)
                                                                                                                        ttGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGCTTAACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgt                 >  3:104482/1‑149 (MQ=255)
                                                                                                                           tgtCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggt              >  7:137657/1‑149 (MQ=255)
                                                                                                                             tCGTGAAAGGGCCGCAAAGCGGCCCCCTCGGTTATCAGCTGGCAGAAATTACTTCAGCTTCACACCTTCGACACTCATATGCAGGACCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGAAGCCGAAGTCGATCAGGTTCAGCgtggtgg            >  8:35499/1‑149 (MQ=255)
                                                                                                                                tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg         >  6:285486/1‑149 (MQ=255)
                                                                                                                                tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg         >  6:290669/1‑149 (MQ=255)
                                                                                                                                tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg         >  6:170634/1‑149 (MQ=255)
                                                                                                                                tGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTAAGCTGGCAGAAATTACTTCTGCTTCACACATTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCATTTGATGCCGAAGTCGTTCAGGTTCAGCGtggtggtgg         >  1:35670/1‑149 (MQ=255)
                                                                                                                                 gATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGc        >  7:44910/1‑149 (MQ=255)
                                                                                                                                  aTAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCg       >  5:10257/1‑149 (MQ=255)
                                                                                                                                       gCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa  <  1:57487/149‑1 (MQ=255)
                                                                                                                                         cGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtgg               >  6:207287/1‑134 (MQ=255)
                                                                                                                                         cGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtgg               <  5:207287/134‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
GAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAA  >  NC_002947_CJ‑RC/5677421‑5677703

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
TGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAG  >  NC_002947_CJ‑RC/5677420‑5677704
                                                                                                                                                   |                                                                                                                                          
TGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTATTGATCTCTACAGGAGCGGCCTTGTGTCCTGATAGGGCCGCAAAGCGGCCC                                                                                                                                       >  NB501016:282:H2V7GAFX2:1:11208:23750:16227/1‑151 (MQ=60)
        CAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTT                                                                                                                                 <  NB501016:282:H2V7GAFX2:2:11303:22298:5392/151‑1 (MQ=60)
             CTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGggcccctcgggttta                                                                                                                             >  NB501016:282:H2V7GAFX2:1:11307:13293:16713/1‑134 (MQ=60)
               GCTGACCCGCAAGGGCTGATAGAAGTACCGAGTGGGCTGCTTCGCAGCCCTATCACGACAAAAGGCAGCTCCTGCAGAGAATTTCTTGATCTCGACAGGAGCGGCCTTGTGTCGTGATAGTGCCGCAAAGCGGCCCCCTCGGTTTTCAGCT                                                                                                                          >  NB501016:282:H2V7GAFX2:1:11107:2381:17498/1‑151 (MQ=60)
               GCTGACCCGCAAGGCCTGATACAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACCCAAGGCCGCTCCTGCAGAGAATTTCTTCATCTCTCCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGC                                                                                                                           <  NB501016:282:H2V7GAFX2:4:21510:20431:12324/150‑1 (MQ=60)
                    CCCGCAAGGGCTGATAGAAGTACCGCGGGGGCTGCTTCGCAGCCCTATCACGGCACAAGGCCGCTCCGGCAGAGAATTTGTTGATCTGTACAGGGGCGGCCTTGTGTGGTGTTAGGGCCGGAAAGCGgccccctcgggtttgatgtggg                                                                                                                      >  NB501016:282:H2V7GAFX2:2:11202:3368:4836/1‑127 (MQ=60)
                    CCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTACTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTGAGCTGGC                                                                                                                       >  NB501016:282:H2V7GAFX2:2:11104:5838:7774/1‑149 (MQ=60)
                     CCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTGTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGTTTTGAGCTGGC                                                                                                                       >  NB501016:282:H2V7GAFX2:2:11110:14398:1769/1‑148 (MQ=60)
                        cacgggcttttagattTACCGCGGGGGGTGCTTCGCAGCCCTATCACGACACATGCCCGCCCCTCCAGATAATTTCTTGATCTCTACTGCCGGGGCCTTGTGTCGTGATTGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAAT                                                                                                                 <  NB501016:282:H2V7GAFX2:4:11603:15825:7240/135‑1 (MQ=46)
                         acgggcttatccaccacccgcgggggcggctacgcccccctaccaccccacacggccgcccctccacacaatttcttgaactaaaccgcAGCGGCCTTGTGGCGTGATAGGGCCGAAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATT                                                                                                                <  NB501016:282:H2V7GAFX2:2:21102:18480:5400/62‑1 (MQ=47)
                                      AGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCAC                                                                                                   <  NB501016:282:H2V7GAFX2:4:21610:12993:12854/151‑1 (MQ=60)
                                                     GCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCT                                                                                                       <  NB501016:282:H2V7GAFX2:4:11511:20178:6752/132‑1 (MQ=60)
                                                               CCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCGCGGGTTTCCGCTGGCAGAAATTACTTGTGCTTGACACCTTCGACACTGATATCCAGGTCC                                                                          >  NB501016:282:H2V7GAFX2:2:21107:5046:14399/1‑151 (MQ=60)
                                                                       GACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGGTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTG                                                                  >  NB501016:282:H2V7GAFX2:1:11208:9157:4907/1‑151 (MQ=60)
                                                                        ACACAAGGCCGCTCCTGCAGATAATTTCTTGATCTCTACTGGAGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGC                                                                 <  NB501016:282:H2V7GAFX2:1:21304:19240:2249/151‑1 (MQ=60)
                                                                         CACACGGCCGCTCTTGCTGACAATTTCTTTACCTCTACTGTAGCGGCCTTGTGTCGTTATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCG                                                                <  NB501016:282:H2V7GAFX2:3:21402:20865:20365/151‑1 (MQ=60)
                                                                         CACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCG                                                                <  NB501016:282:H2V7GAFX2:4:21407:18552:16365/151‑1 (MQ=60)
                                                                             AGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGT                                                            >  NB501016:282:H2V7GAFX2:3:21503:10946:3449/1‑151 (MQ=60)
                                                                                CCGATCCTGCAGATAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGC                                                         <  NB501016:282:H2V7GAFX2:2:11304:19783:18596/151‑1 (MQ=60)
                                                                                       TGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGG                                                  >  NB501016:282:H2V7GAFX2:2:11208:17563:16756/1‑151 (MQ=60)
                                                                                                              CAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAG                           <  NB501016:282:H2V7GAFX2:1:11207:25681:7436/151‑1 (MQ=60)
                                                                                                                GGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGT                         >  NB501016:282:H2V7GAFX2:2:11301:11040:8284/1‑151 (MQ=60)
                                                                                                                GGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGGCACCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGT                         <  NB501016:282:H2V7GAFX2:2:11201:17588:8924/151‑1 (MQ=60)
                                                                                                                       CCTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGT                  <  NB501016:282:H2V7GAFX2:1:11110:11070:8989/151‑1 (MQ=60)
                                                                                                                        CTTGTGTCGTGATAGGGCCGCAAAGCGGCCCCCACGGTTTTCAGCTGGCAGAAATTACTTCTGCTTAACACCTTCGACACTGATATACAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTG                 >  NB501016:282:H2V7GAFX2:2:11211:11937:2794/1‑151 (MQ=60)
                                                                                                                           GTGTCGTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTG              >  NB501016:282:H2V7GAFX2:4:11604:19430:14063/1‑151 (MQ=60)
                                                                                                                                GTGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTAAGCTGGCAGAAATTACTTCTGCTTCACACATTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCATTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGC         >  NB501016:282:H2V7GAFX2:1:11107:26408:17898/1‑151 (MQ=60)
                                                                                                                                 TGATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCG        >  NB501016:282:H2V7GAFX2:4:11409:3889:15946/1‑151 (MQ=60)
                                                                                                                                  GATAGGGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGT       >  NB501016:282:H2V7GAFX2:3:11403:23148:4372/1‑151 (MQ=60)
                                                                                                                                       GGCCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAG  <  NB501016:282:H2V7GAFX2:1:11112:20780:4893/151‑1 (MQ=60)
                                                                                                                                         CCGCAAAGCGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGT               <  NB501016:282:H2V7GAFX2:3:21408:23067:6907/136‑1 (MQ=60)
                                                                                                                                                   |                                                                                                                                          
TGAGCGCCCAGCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAG  >  NC_002947_CJ‑RC/5677420‑5677704

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