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 F24 I1 R1 58 74.9 3151412 98.7% 3110443 147.7

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947_CJ‑RC 2,958,523 C→T A428V (GCA→GTA)  PP_2589 → aldehyde dehydrogenase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947_CJ‑RC2,958,5230CT98.6% 246.9 / NA 72A428V (GCA→GTA) PP_2589aldehyde dehydrogenase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (36/35);  minor base A (1/0);  total (37/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGC  >  NC_002947_CJ‑RC/2958380‑2958665
                                                                                                                                               |                                                                                                                                              
cATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGc                                                                                                                                             <  4:199722/147‑1 (MQ=255)
cATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCt                                                                                                                                           >  5:211374/1‑149 (MQ=255)
 aTCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTc                                                                                                                                          >  5:150533/1‑149 (MQ=255)
 aTCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTc                                                                                                                                          <  3:355291/149‑1 (MQ=255)
 aTCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTc                                                                                                                                          >  2:76689/1‑149 (MQ=255)
      gCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGgtgt                                                                                                                                     >  7:310817/1‑149 (MQ=255)
       ccGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGtgtg                                                                                                                                    >  6:88574/1‑149 (MQ=255)
        cGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTgg                                                                                                                                   <  2:172251/149‑1 (MQ=255)
          aCCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAg                                                                                                                                 >  7:75380/1‑149 (MQ=255)
              tttttGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGcgac                                                                                                                             >  2:25672/1‑149 (MQ=255)
                    ccGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACtt                                                                                                                        >  5:220366/1‑148 (MQ=255)
                     cGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTc                                                                                                                       >  4:305689/1‑148 (MQ=255)
                       gtgtCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAAc                                                                                                                    >  7:99217/1‑149 (MQ=255)
                       gtgtCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAAc                                                                                                                    >  3:94582/1‑149 (MQ=255)
                        tgtCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAAcc                                                                                                                   <  3:9043/149‑1 (MQ=255)
                           cGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAgg                                                                                                                >  6:66592/1‑149 (MQ=255)
                            gAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGt                                                                                                               >  7:232323/1‑149 (MQ=255)
                            gAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGt                                                                                                               >  4:288349/1‑149 (MQ=255)
                                 ccGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTTGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCGCGGTGTGGAGCGACGACTACAACCAGGTGCAcc                                                                                                          >  8:106100/1‑149 (MQ=255)
                                  cGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACTGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCg                                                                                                         <  7:398275/149‑1 (MQ=255)
                                    cATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCggg                                                                                                       >  7:260907/1‑149 (MQ=255)
                                      tGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCggg                                                                                                       >  8:402629/1‑147 (MQ=255)
                                         gCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCggg                                                                                                       >  6:145362/1‑144 (MQ=255)
                                         gCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCggg                                                                                                       <  5:145362/144‑1 (MQ=255)
                                          ccTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAgcgc                                                                                                 >  5:115911/1‑149 (MQ=255)
                                              gcgcgcGAGGAGGTGTTCGGCCCGGTACTGGTGGTCAGCGCCTTCGACACGGAAGAAGAAGTCGTGCGCCTGGCCAACGACAGTATCTATGGGCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGc                                                                                             <  4:268384/149‑1 (MQ=255)
                                                gcgcgAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGt                                                                                           <  2:91230/149‑1 (MQ=255)
                                                gcgcgAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGt                                                                                           <  1:250199/149‑1 (MQ=255)
                                                   cgAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGa                                                                                        <  3:213199/149‑1 (MQ=255)
                                                    gaggagGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGa                                                                                        >  3:184213/1‑148 (MQ=255)
                                                       gaggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGc                                                                                    >  6:326066/1‑149 (MQ=255)
                                                        aggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGcc                                                                                   <  8:99217/149‑1 (MQ=255)
                                                         ggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCg                                                                                  <  8:75380/149‑1 (MQ=255)
                                                         ggTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCg                                                                                  <  8:137606/149‑1 (MQ=255)
                                                           tgtTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCgg                                                                                 >  8:111674/1‑148 (MQ=255)
                                                             ttCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTAc                                                                              >  8:353105/1‑149 (MQ=255)
                                                             ttCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTAc                                                                              >  7:37972/1‑149 (MQ=255)
                                                                   ccGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATc                                                                        <  7:127637/149‑1 (MQ=255)
                                                                   ccGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATc                                                                        >  5:57003/1‑149 (MQ=255)
                                                                      gTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGt                                                                     <  6:57003/149‑1 (MQ=255)
                                                                              ggTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAg                                                             <  6:211374/149‑1 (MQ=255)
                                                                               gTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGt                                                            <  2:317293/149‑1 (MQ=255)
                                                                                tCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTc                                                           >  4:86304/1‑149 (MQ=255)
                                                                                   gcgcCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGAt                                                        <  8:108857/149‑1 (MQ=255)
                                                                                   gcgcCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGAt                                                        >  5:57956/1‑149 (MQ=255)
                                                                                   gcgcCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGAt                                                        <  5:198055/149‑1 (MQ=255)
                                                                                      ccTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATgcg                                                     >  6:335012/1‑149 (MQ=255)
                                                                                      ccTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATgcg                                                     <  8:232323/149‑1 (MQ=255)
                                                                                          cGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCt                                                   <  1:25672/147‑1 (MQ=255)
                                                                                            acacGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGAt                                               <  7:402629/149‑1 (MQ=255)
                                                                                              acGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGAtgt                                             >  5:139476/1‑149 (MQ=255)
                                                                                                ggAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCGTCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGAGGGGCTGGAtgtgt                                           >  2:356948/1‑148 (MQ=255)
                                                                                                ggAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCGCGGGGTGGAGCGACGACTGCAACCAGGTGCACCGGGTAGCGCGGGCGGTGAAAGCCGGTACGGGAGCGGTGCACACTGTCGATGCGCTGGGGGTGa                                           >  2:361911/1‑149 (MQ=255)
                                                                                                  aagaagaagCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACg                                         <  7:116769/149‑1 (MQ=255)
                                                                                                      agaagCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGTTGCACCGGGTAGCGCGGGCGTTGAAACCCGGTACCGAATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGc                                     >  1:74294/1‑149 (MQ=255)
                                                                                                           ccGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTc                                >  8:175121/1‑149 (MQ=255)
                                                                                                             gTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCgg                              <  5:41580/149‑1 (MQ=255)
                                                                                                             gTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCgg                              <  1:349765/149‑1 (MQ=255)
                                                                                                               gcgcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACGGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTg                            <  3:274107/149‑1 (MQ=255)
                                                                                                                cgcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTgg                           <  1:168657/149‑1 (MQ=255)
                                                                                                                 gcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTggg                          <  5:326066/149‑1 (MQ=255)
                                                                                                                 gcCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTggg                          <  7:353105/149‑1 (MQ=255)
                                                                                                                  ccTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTgggg                         <  4:344771/149‑1 (MQ=255)
                                                                                                                    ttgCCAACGACAGCATCTATGGCCTGGAAGCATCGGTGTGGAGCGAATACTTAAAACAGGTGCACCCGGTAGCGCGGGCTTTGAAAGCAGGTACCGTATCGCTGAACTCAGTCGATGCGCTAGATTTGACTGTTCAGTTCGGTggttg                        >  4:222085/3‑145 (MQ=255)
                                                                                                                      gCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGc                       <  1:87524/147‑1 (MQ=255)
                                                                                                                         aaCGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCa                  <  4:327202/149‑1 (MQ=255)
                                                                                                                               aGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCggg            <  4:184213/149‑1 (MQ=255)
                                                                                                                               aGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCggg            >  1:175262/1‑149 (MQ=255)
                                                                                                                                gCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCATTCGGTGGGGGCAAGCAGTCgggg           <  6:84036/149‑1 (MQ=255)
                                                                                                                                       tGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCATCCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGcc    >  8:191939/1‑149 (MQ=255)
                                                                                                                                       tGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGttt        <  4:151196/145‑1 (MQ=255)
                                                                                                                                         gCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCcgc  >  5:156650/1‑149 (MQ=255)
                                                                                                                                               |                                                                                                                                              
CATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGC  >  NC_002947_CJ‑RC/2958380‑2958665

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
ACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCG  >  NC_002947_CJ‑RC/2958379‑2958666
                                                                                                                                                |                                                                                                                                               
ACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTC                                                                                                                                           >  NB501016:282:H2V7GAFX2:3:21402:21108:14170/1‑151 (MQ=60)
 CATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCG                                                                                                                                          >  NB501016:282:H2V7GAFX2:3:11604:15652:3933/1‑151 (MQ=60)
 CATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCG                                                                                                                                          <  NB501016:282:H2V7GAFX2:2:21308:13736:9390/151‑1 (MQ=60)
      AGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTG                                                                                                                                     >  NB501016:282:H2V7GAFX2:4:21507:23477:20278/1‑151 (MQ=60)
          GACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGC                                                                                                                                 >  NB501016:282:H2V7GAFX2:4:11501:19747:18428/1‑151 (MQ=60)
                    GCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTC                                                                                                                        >  NB501016:282:H2V7GAFX2:3:21404:3472:12749/1‑150 (MQ=60)
                       GGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACC                                                                                                                    >  NB501016:282:H2V7GAFX2:4:11506:14378:3436/1‑151 (MQ=60)
                       GGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACC                                                                                                                    >  NB501016:282:H2V7GAFX2:2:11206:6217:9823/1‑151 (MQ=60)
                        GTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCA                                                                                                                   <  NB501016:282:H2V7GAFX2:2:11102:12337:14781/151‑1 (MQ=60)
                            CGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTG                                                                                                               >  NB501016:282:H2V7GAFX2:4:21406:12545:4690/1‑151 (MQ=60)
                                  CCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACTGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGG                                                                                                         <  NB501016:282:H2V7GAFX2:4:21612:25905:1205/151‑1 (MQ=60)
                                    GCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGT                                                                                                       >  NB501016:282:H2V7GAFX2:4:21411:6074:3435/1‑151 (MQ=60)
                                         AGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGT                                                                                                       <  NB501016:282:H2V7GAFX2:3:11603:23646:7251/146‑1 (MQ=60)
                                          GCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCG                                                                                                 >  NB501016:282:H2V7GAFX2:3:11510:17662:6427/1‑151 (MQ=60)
                                                CGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTT                                                                                           <  NB501016:282:H2V7GAFX2:1:21108:18957:19681/151‑1 (MQ=60)
                                                   GCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAA                                                                                        <  NB501016:282:H2V7GAFX2:2:21103:16182:4666/151‑1 (MQ=60)
                                                    CGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAA                                                                                        >  NB501016:282:H2V7GAFX2:2:11309:9604:20023/1‑150 (MQ=60)
                                                             GTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACC                                                                              >  NB501016:282:H2V7GAFX2:4:11407:9308:13158/1‑151 (MQ=60)
                                                                   CCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCG                                                                        <  NB501016:282:H2V7GAFX2:4:11511:12764:3807/151‑1 (MQ=60)
                                                                   CCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCG                                                                        >  NB501016:282:H2V7GAFX2:3:11411:8405:17356/1‑151 (MQ=60)
                                                                            GGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGATCAC                                                                >  NB501016:282:H2V7GAFX2:1:11308:26475:17089/1‑150 (MQ=60)
                                                                                   AGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATG                                                        >  NB501016:282:H2V7GAFX2:3:11412:2432:1748/1‑151 (MQ=60)
                                                                                   AGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATG                                                        <  NB501016:282:H2V7GAFX2:3:11612:17941:2173/151‑1 (MQ=60)
                                                                                          TCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTG                                                   <  NB501016:282:H2V7GAFX2:1:11105:5611:2591/149‑1 (MQ=60)
                                                                                            GACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATG                                               <  NB501016:282:H2V7GAFX2:4:21612:16610:16576/151‑1 (MQ=60)
                                                                                              CACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTG                                             >  NB501016:282:H2V7GAFX2:3:11602:18179:8054/1‑151 (MQ=60)
                                                                                                  GAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGG                                         <  NB501016:282:H2V7GAFX2:4:11509:17109:5488/151‑1 (MQ=60)
                                                                                                      AAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGTTGCACCGGGTAGCGCGGGCGTTGAAACCCGGTACCGAATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCC                                     >  NB501016:282:H2V7GAFX2:1:11201:17158:6191/1‑151 (MQ=60)
                                                                                                             CGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGT                              <  NB501016:282:H2V7GAFX2:3:11408:1755:17558/151‑1 (MQ=60)
                                                                                                             CGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGT                              <  NB501016:282:H2V7GAFX2:1:21304:2980:2744/151‑1 (MQ=60)
                                                                                                               GGCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACGGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGG                            <  NB501016:282:H2V7GAFX2:2:21204:6695:17923/151‑1 (MQ=60)
                                                                                                                GCGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGG                           <  NB501016:282:H2V7GAFX2:1:11306:3046:3314/151‑1 (MQ=60)
                                                                                                                 CGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGG                          <  NB501016:282:H2V7GAFX2:3:21512:4201:17947/151‑1 (MQ=60)
                                                                                                                 CGCCTGGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGG                          <  NB501016:282:H2V7GAFX2:4:21603:24440:13469/151‑1 (MQ=60)
                                                                                                                      GGCCAACGACAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCA                       <  NB501016:282:H2V7GAFX2:1:11203:24257:12324/149‑1 (MQ=60)
                                                                                                                               CAGCATCTATGGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGG            >  NB501016:282:H2V7GAFX2:1:11307:14389:7852/1‑151 (MQ=60)
                                                                                                                                         GGCCTGGTAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCG  >  NB501016:282:H2V7GAFX2:3:11605:21153:4103/1‑151 (MQ=60)
                                                                                                                                                |                                                                                                                                               
ACATCGAGCCGACCATTTTTGCCGGTGTCGAGGGCCGCATGAGCCTGGCGCGCGAGGAGGTGTTCGGCCCGGTACTGGCGGTCAGCGCCTTCGACACGGAAGAAGAAGCCGTGCGCCTGGCCAACGACAGCATCTATGGCCTGGCAGCCTCGGTGTGGAGCGACGACTTCAACCAGGTGCACCGGGTAGCGCGGGCGTTGAAAGCCGGTACCGTATCGGTGAACACAGTCGATGCGCTGGATGTGACGGTGCCGTTCGGTGGGGGCAAGCAGTCGGGGTTTGGCCGCG  >  NC_002947_CJ‑RC/2958379‑2958666

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: