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 F14 I0 R1 206 76.7 3255580 98.8% 3216513 146.3

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947_CJ‑RC 4,348,812 A→G 55.3% intergenic (‑70/+404) PP_3820 ← / ← galU group II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC4,348,8120AG55.3% ‑6.2 / 152.9 123intergenic (‑70/+404)PP_3820/galUgroup II intron‑encoding maturase/UTP‑glucose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (23/32);  new base G (31/37);  total (54/69)
Fisher's exact test for biased strand distribution p-value = 7.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

TCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCA  >  NC_002947_CJ‑RC/4348697‑4348953
                                                                                                                   |                                                                                                                                             
tCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAg                                                                                                                <  6:62245/147‑1 (MQ=11)
  tGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTgg                                                                                                             >  4:44523/1‑148 (MQ=17)
    cATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGATGCGCGGTAGCACAGTg                                                                                                              <  6:137561/145‑1 (MQ=17)
    cATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCaa                                                                                                          <  4:222564/149‑1 (MQ=17)
    cATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCa                                                                                                           <  7:239627/148‑1 (MQ=32)
      ttACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAg                                                                                                         >  5:115252/1‑148 (MQ=17)
      ttACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGc                                                                                                        <  4:117601/149‑1 (MQ=33)
       tACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGc                                                                                                        >  6:37575/1‑148 (MQ=33)
         cGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGt                                                                                                      <  3:187501/148‑1 (MQ=21)
         cGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                                                     >  3:108431/1‑149 (MQ=34)
          ggTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTa                                                                                                     >  3:344014/1‑148 (MQ=34)
           gTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAc                                                                                                    <  8:222777/148‑1 (MQ=21)
           gTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACa                                                                                                   <  1:373861/149‑1 (MQ=21)
               ccAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAt                                                                                                <  1:98965/148‑1 (MQ=21)
               ccAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtc                                                                                               <  3:217535/149‑1 (MQ=34)
               ccAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAt                                                                                                <  7:15659/148‑1 (MQ=34)
                cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtct                                                                                              >  3:12695/1‑149 (MQ=21)
                cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtc                                                                                               <  7:131227/148‑1 (MQ=21)
                cAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGAtc                                                                                               >  5:287518/1‑148 (MQ=34)
                  aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTcc                                                                                            >  4:67567/1‑149 (MQ=21)
                  aCGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGACCTTGGTGTCGCGGTTAACCGCCACTTACTACGGCTCCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTcc                                                                                            >  3:5627/1‑149 (MQ=255)
                    gagaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTccc                                                                                           <  8:140657/148‑1 (MQ=21)
                     agaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCa                                                                                          <  8:401332/148‑1 (MQ=34)
                      gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                        >  1:138357/1‑149 (MQ=34)
                      gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAg                                                                                         <  6:360208/148‑1 (MQ=34)
                      gaCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAg                                                                                         <  5:250677/148‑1 (MQ=34)
                       aCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCTCAGTGGCAAGCGTACAGATCTCCCAgg                                                                                        <  6:20252/148‑1 (MQ=21)
                       aCGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                        <  7:236886/148‑1 (MQ=34)
                           cGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTaa                                                                                    >  8:156568/1‑148 (MQ=34)
                                 cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                             <  2:401329/149‑1 (MQ=21)
                                 cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                             <  4:128/149‑1 (MQ=34)
                                 cTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                                                              <  5:226705/148‑1 (MQ=34)
                                  tACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                             >  2:134277/1‑148 (MQ=34)
                                   aCTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcgc                                                                             >  5:197942/1‑147 (MQ=34)
                                    cTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGa                                                                           <  8:17780/148‑1 (MQ=21)
                                    cTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAc                                                                          >  1:267192/1‑149 (MQ=34)
                                      ggcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAcc                                                                         >  8:288035/1‑148 (MQ=21)
                                       gcggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGAcc                                                                         <  6:287163/147‑1 (MQ=21)
                                        cggcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCtt                                                                       <  8:113612/148‑1 (MQ=21)
                                          gcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGTCGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                     <  7:227109/148‑1 (MQ=21)
                                          gcATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTcc                                                                     <  8:49616/148‑1 (MQ=34)
                                           cATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                                                              <  7:126071/138‑1 (MQ=21)
                                            aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                  >  2:111290/1‑149 (MQ=21)
                                            aTAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                  <  1:355797/149‑1 (MQ=21)
                                            aTAACCTGTTCGGCCCTTGGTGGCGCAGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGc                                                                  <  2:374284/149‑1 (MQ=21)
                                             tAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTcgcg                                                                              >  7:78261/1‑136 (MQ=21)
                                              aaCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCt                                                                 >  2:16078/1‑148 (MQ=34)
                                               aCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGTTAATTCGCGCGACCTTCCTGCtt                                                                <  6:160572/148‑1 (MQ=21)
                                               aCCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTa                                                               >  3:310283/1‑149 (MQ=21)
                                                ccTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                              <  2:194387/149‑1 (MQ=21)
                                                ccTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                              <  3:161334/149‑1 (MQ=34)
                                                ccTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAt                                                              >  4:379756/1‑149 (MQ=34)
                                                   gTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                           >  4:325186/1‑149 (MQ=21)
                                                    ttCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGc                                                            <  4:15872/147‑1 (MQ=21)
                                                    ttCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGcc                                                           >  8:404849/1‑148 (MQ=34)
                                                       ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTg                                                         <  5:384287/147‑1 (MQ=21)
                                                       ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTc                                                       <  1:112960/149‑1 (MQ=34)
                                                       ggCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTc                                                       >  4:333499/1‑149 (MQ=34)
                                                          ccTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTa                                                               <  6:58781/138‑1 (MQ=21)
                                                           cTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGa                                                    >  8:235088/1‑148 (MQ=21)
                                                            ttGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGAt                                                   >  5:409991/1‑148 (MQ=21)
                                                               gtggCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc                                               <  4:344014/149‑1 (MQ=34)
                                                               gtggCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTAc                                               <  3:336993/149‑1 (MQ=34)
                                                                  gcgcGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGt                                             <  7:299050/148‑1 (MQ=34)
                                                                   cgcgGTTAACCGCCACTTACTACGGCTCCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTca                                           <  4:5627/149‑1 (MQ=21)
                                                                     cgGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTcaca                                         <  2:246937/149‑1 (MQ=34)
                                                                     cgGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTca                                           >  7:50175/1‑147 (MQ=34)
                                                                        ttAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCg                                      >  3:217969/1‑149 (MQ=34)
                                                                           aCCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTc                                   >  2:201878/1‑149 (MQ=34)
                                                                            ccGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTc                                   <  8:222823/148‑1 (MQ=34)
                                                                             cGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTcc                                  >  5:144539/1‑148 (MQ=21)
                                                                              gCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCg                                 <  6:181123/148‑1 (MQ=21)
                                                                              gCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGt                                >  3:121087/1‑149 (MQ=21)
                                                                              gCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCg                                 >  6:245184/1‑148 (MQ=34)
                                                                                cACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGc                              <  4:310283/149‑1 (MQ=21)
                                                                                 aCTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCa                             <  1:134277/149‑1 (MQ=34)
                                                                                  cTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCa                             <  7:69580/148‑1 (MQ=21)
                                                                                  cTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGc                              >  7:356049/1‑147 (MQ=34)
                                                                                  cTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCa                             <  6:287518/148‑1 (MQ=34)
                                                                                      ctacGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCTGTAGCACAGTGGCAAGCGTACAGATCTCCCAGTGTAATTCGCGCGACCTTCCTGCTTATTCCTGTCGGATCTACGTCACAGCGTTCCGTGCATGt                          <  7:356755/147‑1 (MQ=11)
                                                                                       tacGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTAtt                       >  3:149030/1‑149 (MQ=21)
                                                                                         cGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTAtt                       <  7:347626/147‑1 (MQ=21)
                                                                                          ggCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTggg                    >  4:188399/1‑149 (MQ=21)
                                                                                          ggCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCACGGTAATTCGAGCGACCTTCCTGCTTAAGACTGTCGGATCTACGTCACAGCGTTCCGAGCAAGTATTggg                    >  2:368476/1‑149 (MQ=11)
                                                                                             ttCGGCTGACTTCTGAACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCt                  >  5:172489/1‑148 (MQ=21)
                                                                                                ggCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTgg              <  2:267192/149‑1 (MQ=34)
                                                                                                  cTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGc             >  5:311143/1‑148 (MQ=21)
                                                                                                     aCTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAt          >  5:255565/1‑148 (MQ=21)
                                                                                                     aCTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAt          <  4:108431/148‑1 (MQ=34)
                                                                                                     aCTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAt          >  5:349603/1‑148 (MQ=34)
                                                                                                        tCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtg       >  6:380324/1‑148 (MQ=21)
                                                                                                         cTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt      <  4:220719/148‑1 (MQ=21)
                                                                                                         cTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt      <  3:179939/148‑1 (MQ=21)
                                                                                                         cTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGAttt        <  8:189545/146‑1 (MQ=34)
                                                                                                         cTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt      <  4:217969/148‑1 (MQ=34)
                                                                                                         cTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtgt      <  1:333501/148‑1 (MQ=34)
                                                                                                         cTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTtg       <  6:349009/147‑1 (MQ=34)
                                                                                                            cACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGc    <  8:340908/147‑1 (MQ=21)
                                                                                                             aCGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCa  >  6:278162/1‑148 (MQ=21)
                                                                                                                   aTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAgg                                                                                        <  1:328734/56‑1 (MQ=255)
                                                                                                                   |                                                                                                                                             
TCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCCA  >  NC_002947_CJ‑RC/4348697‑4348953

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ 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
CGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCC  >  NC_002947_CJ‑RC/4348668‑4348952
                                                                                                                                                |                                                                                                                                            
CGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACCGCTTCGGCTGACTTCTGCACGCTCATCCCGT                                                                                                                                        >  NB501016:282:H2V7GAFX2:3:11404:15690:18901/1‑151 (MQ=0)
 GGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTC                                                                                                                                       <  NB501016:282:H2V7GAFX2:4:11612:20129:3200/151‑1 (MQ=0)
  GGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCG                                                                                                                                      >  NB501016:282:H2V7GAFX2:4:11611:14683:14519/1‑151 (MQ=0)
       ACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTT                                                                                                                                 <  NB501016:282:H2V7GAFX2:3:21501:15046:4193/151‑1 (MQ=0)
       ACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTT                                                                                                                                 <  NB501016:282:H2V7GAFX2:4:11402:6413:10815/151‑1 (MQ=0)
         AGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCG                                                                                                                               <  NB501016:282:H2V7GAFX2:4:11410:3894:6765/151‑1 (MQ=0)
             TCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACG                                                                                                                            <  NB501016:282:H2V7GAFX2:4:21611:18196:5469/150‑1 (MQ=0)
             TCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACCTCTGCACGCTCATCCCGTCGCCTTTCGACG                                                                                                                            <  NB501016:282:H2V7GAFX2:1:11311:12076:11712/150‑1 (MQ=0)
                   GTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGT                                                                                                                      >  NB501016:282:H2V7GAFX2:4:11609:5851:7761/1‑150 (MQ=0)
                     GGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGC                                                                                                                   >  NB501016:282:H2V7GAFX2:2:11203:5774:15694/1‑151 (MQ=0)
                      GGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGC                                                                                                                   >  NB501016:282:H2V7GAFX2:3:11507:13371:4203/1‑150 (MQ=0)
                                GCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAG                                                                                                        <  NB501016:282:H2V7GAFX2:4:21405:8580:18821/151‑1 (MQ=4)
                                 CATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGC                                                                                                       >  NB501016:282:H2V7GAFX2:3:11509:17247:5230/1‑151 (MQ=6)
                                     ACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTA                                                                                                    <  NB501016:282:H2V7GAFX2:2:11309:8768:19152/150‑1 (MQ=21)
                                     ACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAC                                                                                                   >  NB501016:282:H2V7GAFX2:2:11208:22481:9163/1‑151 (MQ=22)
                                      CGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTAC                                                                                                   >  NB501016:282:H2V7GAFX2:2:21304:3499:16782/1‑150 (MQ=22)
                                       GGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAG                                                                                                 <  NB501016:282:H2V7GAFX2:1:21307:17265:6445/151‑1 (MQ=21)
                                           ACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATC                                                                                              <  NB501016:282:H2V7GAFX2:1:11205:13043:9293/150‑1 (MQ=21)
                                           ACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCT                                                                                             <  NB501016:282:H2V7GAFX2:4:11403:12053:12504/151‑1 (MQ=22)
                                           ACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCT                                                                                             <  NB501016:282:H2V7GAFX2:2:21103:13169:3305/151‑1 (MQ=22)
                                           ACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCT                                                                                             >  NB501016:282:H2V7GAFX2:3:21503:25359:7318/1‑151 (MQ=22)
                                            CCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTC                                                                                            <  NB501016:282:H2V7GAFX2:4:11510:6378:18683/151‑1 (MQ=21)
                                            CCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTC                                                                                            >  NB501016:282:H2V7GAFX2:2:11103:4642:7627/1‑151 (MQ=21)
                                              AACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGACCTTGGTGTCGCGGTTAACCGCCACTTACTACGGCTCCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGAGGCAAGCGTACAGATCACCC                                                                                          >  NB501016:282:H2V7GAFX2:2:11102:6166:1634/1‑151 (MQ=18)
                                                  AGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGG                                                                                      <  NB501016:282:H2V7GAFX2:3:21408:20904:12076/151‑1 (MQ=22)
                                                  AGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGG                                                                                      >  NB501016:282:H2V7GAFX2:1:11212:8289:6930/1‑151 (MQ=22)
                                                   GACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGT                                                                                     <  NB501016:282:H2V7GAFX2:1:21209:11387:3678/151‑1 (MQ=21)
                                                   GACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGT                                                                                     <  NB501016:282:H2V7GAFX2:4:21405:10500:9958/151‑1 (MQ=22)
                                                             CCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCG                                                                           <  NB501016:282:H2V7GAFX2:3:21404:5465:3434/151‑1 (MQ=22)
                                                              CTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCG                                                                           >  NB501016:282:H2V7GAFX2:3:11610:22630:17737/1‑150 (MQ=22)
                                                                ACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACC                                                                        >  NB501016:282:H2V7GAFX2:1:21111:14128:5322/1‑151 (MQ=22)
                                                                      GGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGTCGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTT                                                                  <  NB501016:282:H2V7GAFX2:4:21403:5084:17216/151‑1 (MQ=20)
                                                                      GGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCG                                                                           <  NB501016:282:H2V7GAFX2:4:11510:23463:1487/142‑1 (MQ=21)
                                                                        CATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCT                                                                <  NB501016:282:H2V7GAFX2:1:21303:19799:12043/151‑1 (MQ=21)
                                                                        CATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCG                                                                           >  NB501016:282:H2V7GAFX2:4:11501:22288:16577/1‑140 (MQ=22)
                                                                           AACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAT                                                             >  NB501016:282:H2V7GAFX2:2:21210:11580:4111/1‑151 (MQ=21)
                                                                            ACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATG                                                            <  NB501016:282:H2V7GAFX2:2:11306:26571:11318/151‑1 (MQ=22)
                                                                            ACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATG                                                            <  NB501016:282:H2V7GAFX2:2:11305:1175:11323/151‑1 (MQ=22)
                                                                                   CGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTC                                                      <  NB501016:282:H2V7GAFX2:3:21608:9820:3741/150‑1 (MQ=21)
                                                                                   CGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCG                                                     <  NB501016:282:H2V7GAFX2:1:11207:8052:17593/151‑1 (MQ=22)
                                                                                     GCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATG                                                            >  NB501016:282:H2V7GAFX2:3:11411:25968:10775/1‑142 (MQ=21)
                                                                                       CCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATC                                                 >  NB501016:282:H2V7GAFX2:1:21310:21920:12901/1‑151 (MQ=21)
                                                                                       CCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATC                                                 >  NB501016:282:H2V7GAFX2:3:21612:13577:18171/1‑151 (MQ=21)
                                                                                           GGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACG                                             <  NB501016:282:H2V7GAFX2:2:21303:19482:10992/151‑1 (MQ=22)
                                                                                              GGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCA                                          <  NB501016:282:H2V7GAFX2:4:21503:3250:13767/151‑1 (MQ=22)
                                                                                                CGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCAC                                         >  NB501016:282:H2V7GAFX2:4:11409:16501:6851/1‑150 (MQ=22)
                                                                                                   GGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCG                                     <  NB501016:282:H2V7GAFX2:4:21501:15150:15144/151‑1 (MQ=21)
                                                                                                    GTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGT                                    >  NB501016:282:H2V7GAFX2:2:21103:7707:4958/1‑151 (MQ=22)
                                                                                                        ACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCG                                >  NB501016:282:H2V7GAFX2:3:11602:15871:3560/1‑151 (MQ=21)
                                                                                                          CGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTG                              >  NB501016:282:H2V7GAFX2:2:11210:13130:12725/1‑151 (MQ=21)
                                                                                                             CACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTAAGCCTGTCGGATCTACGTCACAGCGTTCCGTGCA                            >  NB501016:282:H2V7GAFX2:2:21202:1525:10676/1‑150 (MQ=20)
                                                                                                             CACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAA                           <  NB501016:282:H2V7GAFX2:1:11211:22195:12291/151‑1 (MQ=22)
                                                                                                             CACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCA                            >  NB501016:282:H2V7GAFX2:4:21601:1797:9175/1‑150 (MQ=22)
                                                                                                              ACTTACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAG                          <  NB501016:282:H2V7GAFX2:4:11412:5753:7805/151‑1 (MQ=21)
                                                                                                              ACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAG                          <  NB501016:282:H2V7GAFX2:4:11603:5389:9022/151‑1 (MQ=22)
                                                                                                                  ACTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCTGTAGCACAGTGGCAAGCGTACAGATCTCCCAGTGTAATTCGCGCGACCTTCCTGCTTATTCCTGTCGGATCTACGTCACAGCGTTCCGTGCATGTAT                       <  NB501016:282:H2V7GAFX2:4:21601:11022:11685/150‑1 (MQ=19)
                                                                                                                   CTACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTG                     >  NB501016:282:H2V7GAFX2:2:11303:20668:9161/1‑151 (MQ=21)
                                                                                                                     ACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGG                   <  NB501016:282:H2V7GAFX2:2:11209:16157:13441/151‑1 (MQ=21)
                                                                                                                     ACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGG                    <  NB501016:282:H2V7GAFX2:3:21505:6618:17117/150‑1 (MQ=21)
                                                                                                                     ACGGCTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGG                    <  NB501016:282:H2V7GAFX2:4:21511:19066:19728/150‑1 (MQ=21)
                                                                                                                        GCTTCGGCTGACTTCTGAACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTT                >  NB501016:282:H2V7GAFX2:3:11606:12449:13227/1‑151 (MQ=21)
                                                                                                                         CTTCGGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTT               >  NB501016:282:H2V7GAFX2:4:11408:6745:14786/1‑151 (MQ=21)
                                                                                                                             GGCTGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCG           >  NB501016:282:H2V7GAFX2:3:21507:5378:14581/1‑151 (MQ=21)
                                                                                                                                TGACTTCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATT        >  NB501016:282:H2V7GAFX2:3:21409:2506:10511/1‑151 (MQ=21)
                                                                                                                                TGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATT        >  NB501016:282:H2V7GAFX2:3:21602:6166:9200/1‑151 (MQ=22)
                                                                                                                                     TCTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTG    <  NB501016:282:H2V7GAFX2:2:11308:17214:13741/150‑1 (MQ=21)
                                                                                                                                     TCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTG    <  NB501016:282:H2V7GAFX2:1:21210:23911:19870/150‑1 (MQ=22)
                                                                                                                                     TCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTG    <  NB501016:282:H2V7GAFX2:2:11112:16673:6938/150‑1 (MQ=22)
                                                                                                                                      CTGCACGCTCGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCC  >  NB501016:282:H2V7GAFX2:3:21402:3112:11363/1‑151 (MQ=21)
                                                                                                                                               CGTCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTG                                                                  >  NB501016:282:H2V7GAFX2:4:21410:23224:13852/1‑78 (MQ=28)
                                                                                                                                               CATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGG                                                                                      <  NB501016:282:H2V7GAFX2:1:21210:22094:3955/58‑1 (MQ=33)
                                                                                                                                                |                                                                                                                                            
CGGGATGACAGTTTCCGCTGTGGGAAATTTCTGCATTACGGTTACCAACGAGACGGCGACGCCTACTGGCGGCATAACCTGTTCGGCCCTTGGTGGCGCGGTTAACCGCCACTTACTACGGCTTCGGCTGACTTCTGCACGCTCATCCCGTCGCCTTTCGACGCGCGGTAGCACAGTGGCAAGCGTACAGATCTCCCAGGGTAATTCGCGCGACCTTCCTGCTTATGCCTGTCGGATCTACGTCACAGCGTTCCGTGCAAGTATTGGGCTTTGGCGATTTGTGCC  >  NC_002947_CJ‑RC/4348668‑4348952

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