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
A1 F1 I2 R1 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,999,772 T→C H4R (CAC→CGC)  b1922 ← sigma factor 28

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,999,7720TC99.2% 391.9 / ‑5.5 127H4R (CAC→CGC) b1922sigma factor 28
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (40/86);  total (40/87)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 136 total aligned reads displayed.

AGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCT  >  NC_000913/1999682‑1999861
                                                                                          |                                                                                         
aGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAgcgt                                                                                         <  2:2630153/93‑1 (MQ=255)
 gATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGtt                                                                                        >  2:1724972/1‑93 (MQ=255)
  aTCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGttt                                                                                       <  1:3976355/93‑1 (MQ=255)
    cGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTAt                                                                                     >  1:1543503/1‑93 (MQ=255)
     gTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATc                                                                                    <  1:3131407/93‑1 (MQ=255)
     gTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGATTATc                                                                                    >  2:3389079/1‑93 (MQ=255)
        aaGTTCCATGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCCGGGCATAACGCTGCCACAGCGAGCGTTTATTCAt                                                                                 >  2:4374128/1‑93 (MQ=255)
        aaGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCCCCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCAt                                                                                 <  2:934171/93‑1 (MQ=255)
         aGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCAtt                                                                                <  2:2791033/93‑1 (MQ=255)
          gTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTa                                                                               >  2:3536032/1‑93 (MQ=255)
           ttCCACGCTCTCGGGCTGTCGAACCTGCAGGCGCTATGCTTCGTGACGCACCAGGGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTac                                                                              <  2:3289395/93‑1 (MQ=255)
              cACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACAcc                                                                           <  1:67715/93‑1 (MQ=255)
                cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGTCATAACGCTGCCACAGCGAGCGTTTATCCATTACACCtt                                                                         <  2:493824/93‑1 (MQ=255)
                cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCtt                                                                         <  1:1681493/93‑1 (MQ=255)
                cGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCtt                                                                         <  1:2046877/93‑1 (MQ=255)
                 gCTCGCGGGCAGTCGAACCTGCAGGCACAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTc                                                                        <  2:4142500/93‑1 (MQ=255)
                  cTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCa                                                                       <  1:833108/93‑1 (MQ=255)
                    cgcgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGc                                                                     <  1:4247389/93‑1 (MQ=255)
                      cgGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCgg                                                                   >  1:2407553/1‑93 (MQ=255)
                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCgga                                                                  >  1:353432/1‑92 (MQ=255)
                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  <  1:3414957/93‑1 (MQ=255)
                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  >  1:995785/1‑93 (MQ=255)
                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  <  2:3264186/93‑1 (MQ=255)
                       gggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGt                                                                  >  1:4503156/1‑93 (MQ=255)
                        ggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGta                                                                 <  2:2296801/93‑1 (MQ=255)
                        ggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGta                                                                 <  1:2145145/93‑1 (MQ=255)
                        ggCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGta                                                                 <  2:3378865/93‑1 (MQ=255)
                         gCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGtat                                                                >  2:2461244/1‑93 (MQ=255)
                          cAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGtata                                                               <  1:402898/93‑1 (MQ=255)
                           aGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATag                                                              >  1:1913426/1‑93 (MQ=255)
                            gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  2:2753812/93‑1 (MQ=255)
                            gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  1:2770865/93‑1 (MQ=255)
                            gTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATaga                                                             <  2:2680261/93‑1 (MQ=255)
                              cGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGt                                                           >  1:973298/1‑93 (MQ=255)
                              cGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGt                                                           <  2:1668284/93‑1 (MQ=255)
                                 aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATCACGCTGCCACAGCGAGCGTTTATACATTACACCTTCAGCGGTATAGAGTGaa                                                        >  1:666782/1‑93 (MQ=255)
                                 aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                                        >  1:3603444/1‑93 (MQ=255)
                                 aCCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGaa                                                        <  2:2547868/93‑1 (MQ=255)
                                   cTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAAtt                                                      >  2:915876/1‑93 (MQ=255)
                                    tGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                                     >  1:347867/1‑93 (MQ=255)
                                    tGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTc                                                     >  2:1009714/1‑93 (MQ=255)
                                     gCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCa                                                    <  1:2632528/93‑1 (MQ=255)
                                     gCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCa                                                    <  1:2609669/93‑1 (MQ=255)
                                      gaGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  2:1167187/92‑1 (MQ=255)
                                      cAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  2:1881100/93‑1 (MQ=255)
                                      cAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  1:570666/93‑1 (MQ=255)
                                      cAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAc                                                   <  2:3576128/93‑1 (MQ=255)
                                        ggCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACgt                                                 >  2:2082185/1‑92 (MQ=255)
                                        ggCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                                                 >  1:3451174/1‑93 (MQ=255)
                                        ggCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGa                                                 <  2:336379/93‑1 (MQ=255)
                                          cgcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATa                                               <  1:2554281/93‑1 (MQ=255)
                                          cgcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATa                                               <  2:4653906/93‑1 (MQ=255)
                                           gcAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATaa                                              <  2:3766298/93‑1 (MQ=255)
                                             aaTGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAc                                            >  1:363425/1‑93 (MQ=255)
                                              aTGCTTGGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                                           <  2:3501433/93‑1 (MQ=255)
                                              aTGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACa                                           <  1:1093040/93‑1 (MQ=255)
                                                  ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       <  2:924807/93‑1 (MQ=255)
                                                  ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       >  2:1616752/1‑93 (MQ=255)
                                                  ttCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGccc                                       >  2:4272484/1‑93 (MQ=255)
                                                   tCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCt                                      <  2:4675640/93‑1 (MQ=255)
                                                      tGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCg                                   >  2:1541824/1‑93 (MQ=255)
                                                       gACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt                                  <  2:1123545/93‑1 (MQ=255)
                                                       gACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt                                  <  2:1414623/93‑1 (MQ=255)
                                                       gACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGt                                  >  1:4303372/1‑93 (MQ=255)
                                                          gCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtat                               >  2:853683/1‑93 (MQ=255)
                                                          ccACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtat                               <  1:1432310/92‑1 (MQ=255)
                                                           cACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtata                              <  2:1973609/93‑1 (MQ=255)
                                                           cACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtata                              <  1:912135/93‑1 (MQ=255)
                                                           cACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTtata                              <  1:4241193/93‑1 (MQ=255)
                                                            aCCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTatat                             <  1:4548037/93‑1 (MQ=255)
                                                             ccAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATg                            <  1:926725/93‑1 (MQ=255)
                                                              cAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGa                           <  2:397175/93‑1 (MQ=255)
                                                                gCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGt                         <  2:667947/93‑1 (MQ=255)
                                                                 cGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGtt                        >  1:300450/1‑93 (MQ=255)
                                                                  gggACATAACGCTGCCACAGCGAGCGTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTa                       <  1:4394493/93‑1 (MQ=255)
                                                                  gggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTa                       >  1:704881/1‑93 (MQ=255)
                                                                   ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      >  2:1103496/1‑93 (MQ=255)
                                                                   ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      >  1:1286226/1‑93 (MQ=255)
                                                                   ggACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAt                      >  1:139136/1‑93 (MQ=255)
                                                                    gACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATc                     <  2:4254964/93‑1 (MQ=255)
                                                                    gACAAACCGCGGCCCGAGCGAGCGTTTATCCATTACACCTGCAGCGGGATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATc                     <  1:1238084/93‑1 (MQ=255)
                                                                     aCATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCg                    <  2:3086361/93‑1 (MQ=255)
                                                                     aCATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCg                    <  1:3838985/93‑1 (MQ=255)
                                                                        tAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCa                 <  1:337392/93‑1 (MQ=255)
                                                                        tAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCa                 <  1:1012174/93‑1 (MQ=255)
                                                                          aCGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATg               <  1:1361663/93‑1 (MQ=255)
                                                                           cGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGa              >  2:2017223/1‑93 (MQ=255)
                                                                            gCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGAt             <  2:2186400/93‑1 (MQ=255)
                                                                            gCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGAt             <  1:771474/93‑1 (MQ=255)
                                                                               gCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTAt          <  2:3887231/93‑1 (MQ=255)
                                                                                ccACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATc         <  1:1411657/93‑1 (MQ=255)
                                                                                ccACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATc         <  1:4123393/93‑1 (MQ=255)
                                                                                  acaGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATTCg       >  1:1039636/1‑93 (MQ=255)
                                                                                  acaGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCg       <  2:4120368/93‑1 (MQ=255)
                                                                                   caGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGt      <  1:2464938/93‑1 (MQ=255)
                                                                                    agcgagcgTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGtt     <  1:2016715/93‑1 (MQ=255)
                                                                                     gCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGttt    <  2:2675104/93‑1 (MQ=255)
                                                                                     gCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGttt    >  1:875925/1‑93 (MQ=255)
                                                                                       gAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCt  <  1:3291224/93‑1 (MQ=255)
                                                                                       gAGCGTTTATCCAGTACTCCTTCAGCGGTATAGATTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCTTTTCt  <  2:1403510/93‑1 (MQ=255)
                                                                                          |                                                                                         
AGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCT  >  NC_000913/1999682‑1999861

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGC  >  NC_000913/1999681‑1999863
                                                                                           |                                                                                           
CAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGT                                                                                           <  SRR4302535.1596992/94‑1 (MQ=60)
   ATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTA                                                                                        <  SRR4302535.4003520/94‑1 (MQ=60)
    TCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTAT                                                                                       >  SRR4302535.1554625/1‑94 (MQ=60)
     CGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATC                                                                                      >  SRR4302535.1384400/1‑94 (MQ=60)
      GTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCC                                                                                     <  SRR4302535.3153216/94‑1 (MQ=60)
               CACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCT                                                                            <  SRR4302535.68219/94‑1 (MQ=60)
                 CGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTC                                                                          <  SRR4302535.1693671/94‑1 (MQ=60)
                 CGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTC                                                                          <  SRR4302535.2061615/94‑1 (MQ=60)
                 CGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTC                                                                          >  SRR4302535.3239036/1‑94 (MQ=60)
                   CTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAG                                                                        <  SRR4302535.838946/94‑1 (MQ=60)
                     CGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCG                                                                      <  SRR4302535.4276479/94‑1 (MQ=60)
                      GCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGG                                                                     >  SRR4302535.2424882/1‑94 (MQ=60)
                       CGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGT                                                                    >  SRR4302535.1002854/1‑94 (MQ=60)
                       CGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGT                                                                    >  SRR4302535.4534076/1‑94 (MQ=60)
                       CGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGA                                                                    >  SRR4302535.355870/1‑94 (MQ=60)
                        GGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTA                                                                   <  SRR4302535.3438713/94‑1 (MQ=60)
                         GGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTAT                                                                  <  SRR4302535.2160577/94‑1 (MQ=60)
                           CAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAG                                                                >  SRR4302535.1927216/1‑94 (MQ=60)
                           CAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAG                                                                <  SRR4302535.1479974/94‑1 (MQ=60)
                           CAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAG                                                                <  SRR4302535.4480369/94‑1 (MQ=60)
                           CAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAG                                                                <  SRR4302535.405722/94‑1 (MQ=60)
                             GTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAG                                                              <  SRR4302535.2790596/94‑1 (MQ=60)
                             GTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAG                                                              <  SRR4302535.937367/94‑1 (MQ=60)
                              TCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGT                                                             >  SRR4302535.980205/1‑94 (MQ=60)
                                 AACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATCACGCTGCCACAGCGAGCGTTTATACATTACACCTTCAGCGGTATAGAGTGAA                                                          >  SRR4302535.671432/1‑94 (MQ=60)
                                 AACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAA                                                          >  SRR4302535.3628447/1‑94 (MQ=60)
                                    CTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTC                                                       >  SRR4302535.350261/1‑94 (MQ=60)
                                     TGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCA                                                      <  SRR4302535.716256/94‑1 (MQ=60)
                                      GCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAC                                                     <  SRR4302535.2628396/94‑1 (MQ=60)
                                      GCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCAC                                                     <  SRR4302535.2651392/94‑1 (MQ=60)
                                       CAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACG                                                    <  SRR4302535.574673/94‑1 (MQ=60)
                                        AGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGA                                                   >  SRR4302535.3475162/1‑94 (MQ=60)
                                         GGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGAT                                                  <  SRR4302535.2683462/94‑1 (MQ=60)
                                         GGCGCAATGCTTCGTGACGCACCAGCGGGACAGAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGAT                                                  <  SRR4302535.3755887/94‑1 (MQ=60)
                                          GCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATA                                                 <  SRR4302535.894485/94‑1 (MQ=60)
                                           CGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAA                                                <  SRR4302535.2572670/94‑1 (MQ=60)
                                             CAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAAC                                              >  SRR4302535.365929/1‑94 (MQ=60)
                                               ATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAG                                            <  SRR4302535.1100823/94‑1 (MQ=60)
                                                       TGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCGTGCGT                                    <  SRR4302535.4301672/94‑1 (MQ=60)
                                                       TGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGT                                    >  SRR4302535.4332815/1‑94 (MQ=60)
                                                       TGACGCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGT                                    >  SRR4302535.2018784/1‑94 (MQ=60)
                                                           CCACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATA                                <  SRR4302535.1442536/94‑1 (MQ=60)
                                                            CACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATAT                               <  SRR4302535.918622/94‑1 (MQ=60)
                                                            CACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATAT                               <  SRR4302535.4270222/94‑1 (MQ=60)
                                                             ACCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATG                              <  SRR4302535.4579133/94‑1 (MQ=60)
                                                              CCAGCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGA                             <  SRR4302535.933289/94‑1 (MQ=60)
                                                                 GCGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTT                          >  SRR4302535.302543/1‑94 (MQ=60)
                                                                  CGGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTA                         >  SRR4302535.709830/1‑94 (MQ=60)
                                                                   GGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACCCCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAT                        <  SRR4302535.4424565/94‑1 (MQ=60)
                                                                   GGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAT                        >  SRR4302535.140111/1‑94 (MQ=60)
                                                                   GGGACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTAT                        >  SRR4302535.1295393/1‑94 (MQ=60)
                                                                     gacaaaccgcggcccgAGCGAGCGTTTATCCATTACACCTGCAGCGGGATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCG                      <  SRR4302535.1246926/78‑1 (MQ=60)
                                                                     GACATACCGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCG                      <  SRR4302535.4084507/94‑1 (MQ=60)
                                                                     GACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATTGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCG                      <  SRR4302535.2536305/94‑1 (MQ=60)
                                                                     GACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCG                      <  SRR4302535.3398063/94‑1 (MQ=60)
                                                                     GACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCG                      >  SRR4302535.1226869/1‑94 (MQ=60)
                                                                      ACATAACGCTGCCACAGCGCGCGTTTATCCATTACACCGTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGG                     <  SRR4302535.1519062/94‑1 (MQ=60)
                                                                      ACATAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGG                     <  SRR4302535.3865446/94‑1 (MQ=60)
                                                                         TAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCAT                  <  SRR4302535.1019369/94‑1 (MQ=60)
                                                                         TAACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCAT                  <  SRR4302535.339728/94‑1 (MQ=60)
                                                                           ACGCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGA                <  SRR4302535.1371402/94‑1 (MQ=60)
                                                                             GCTGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATT              <  SRR4302535.776904/94‑1 (MQ=60)
                                                                               TGCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTAT            >  SRR4302535.3524647/1‑94 (MQ=60)
                                                                                GCCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATC           <  SRR4302535.593874/94‑1 (MQ=60)
                                                                                 CCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCC          <  SRR4302535.4151611/94‑1 (MQ=60)
                                                                                 CCACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCC          <  SRR4302535.1421734/94‑1 (MQ=60)
                                                                                  CACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATTCG         >  SRR4302535.1047012/1‑94 (MQ=60)
                                                                                  CACAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCG         <  SRR4302535.2076454/94‑1 (MQ=60)
                                                                                    CAGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTT       <  SRR4302535.2482712/94‑1 (MQ=60)
                                                                                     AGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT      <  SRR4302535.2031241/94‑1 (MQ=60)
                                                                                     AGCGAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTT      >  SRR4302535.882130/1‑94 (MQ=60)
                                                                                        GAGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTG   <  SRR4302535.3314108/94‑1 (MQ=60)
                                                                                         AGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGC  >  SRR4302535.2041833/1‑94 (MQ=60)
                                                                                         AGCGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGC  >  SRR4302535.4075227/1‑94 (MQ=60)
                                                                                           |                                                                                           
CAGATCGTCAAGTTCCACGCTCGCGGGCAGTCGAACCTGCAGGCGCAATGCTTCGTGACGCACCAGCGGGACATAACGCTGCCACAGCGAGTGTTTATCCATTACACCTTCAGCGGTATAGAGTGAATTCACGATAAACAGCCCTGCGTTATATGAGTTATCGGCATGATTATCCGTTTCTGC  >  NC_000913/1999681‑1999863

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 0 ≤ ATCG/ATCG < 2 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑