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 I1 R1 1357 179.3 9500450 96.1% 9129932 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,448,171 T→C N147S (AAC→AGC)  b2333 ← putative fimbrial‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,448,1710TC97.0% 411.9 / 1.9 135N147S (AAC→AGC) b2333putative fimbrial‑like protein
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (63/68);  total (66/69)
Fisher's exact test for biased strand distribution p-value = 3.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-01

Only 100 of 138 total aligned reads displayed.

TCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACT  >  NC_000913/2448081‑2448260
                                                                                          |                                                                                         
tCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTa                                                                                         >  2:2732970/1‑93 (MQ=255)
 cAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTAt                                                                                        >  2:2710257/1‑93 (MQ=255)
 cAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTAt                                                                                        <  2:4336353/93‑1 (MQ=255)
  aaaTGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTAtt                                                                                       <  2:1577418/93‑1 (MQ=255)
  aaaTGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTAtt                                                                                       >  1:2072958/1‑93 (MQ=255)
  aaaTGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTAtt                                                                                       <  2:4630621/93‑1 (MQ=255)
  aaaTGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTAtt                                                                                       >  1:4037431/1‑93 (MQ=255)
    aggTCTCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGc                                                                                     <  1:3091895/91‑1 (MQ=255)
    aTGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGc                                                                                     >  1:2834965/1‑93 (MQ=255)
    aTGTCGCAGTTAAAGAGGCGGTAAATTCACAGATCGTCAGATCACGGCCGCTTCTCGCCTGTCACCGGCATCTGAAATTGATATCGCGATTGc                                                                                     >  2:3851359/1‑93 (MQ=255)
      gTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCa                                                                                   >  1:4553315/1‑93 (MQ=255)
      gTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCa                                                                                   <  2:2853605/93‑1 (MQ=255)
      gTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCa                                                                                   >  1:652715/1‑93 (MQ=255)
       tCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCAt                                                                                  >  1:4695355/1‑93 (MQ=255)
       tCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCAt                                                                                  >  2:2880094/1‑93 (MQ=255)
       tCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCAt                                                                                  <  1:3241860/93‑1 (MQ=255)
        cGCAGTTAAAGAGGCGGTAAATTCTCCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCAtt                                                                                 <  1:677984/93‑1 (MQ=255)
         gCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTGGCCAGGCTCTGAAATTGATATCGCTATTGCCATTg                                                                                >  1:4097122/1‑93 (MQ=255)
         gCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTg                                                                                >  1:1712375/1‑93 (MQ=255)
         gCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTg                                                                                <  1:999785/93‑1 (MQ=255)
           aGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGt                                                                              <  2:4733928/93‑1 (MQ=255)
            gTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGtt                                                                             >  1:3725777/1‑93 (MQ=255)
             ttAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTa                                                                            <  2:3888567/93‑1 (MQ=255)
             ttAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTa                                                                            <  1:3678546/93‑1 (MQ=255)
             ttAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTa                                                                            <  1:3511868/93‑1 (MQ=255)
             ttAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTa                                                                            <  2:1900841/93‑1 (MQ=255)
              tAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAg                                                                           >  2:80157/1‑93 (MQ=255)
              tAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAg                                                                           >  2:619842/1‑93 (MQ=255)
              tAAAGAGCAGGTCAATTCACCGATCGTCACATCACGGCCGCTTGTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAg                                                                           <  2:241920/93‑1 (MQ=255)
               aaaGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAgc                                                                          <  2:4380888/93‑1 (MQ=255)
                aaGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAgcg                                                                         >  2:499290/1‑93 (MQ=255)
                 agagGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAgcgc                                                                        <  1:2301813/93‑1 (MQ=255)
                 agagGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATGGCCATTGGTTAgcgc                                                                        <  1:2301814/93‑1 (MQ=255)
                     gCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAg                                                                    >  1:3373101/1‑93 (MQ=255)
                     gCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAg                                                                    <  2:3812696/93‑1 (MQ=255)
                       ggTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGt                                                                  <  2:3879179/93‑1 (MQ=255)
                       ggTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGt                                                                  >  1:3090591/1‑93 (MQ=255)
                       ggTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGt                                                                  <  1:2202908/93‑1 (MQ=255)
                        gTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTc                                                                 >  2:4668722/1‑93 (MQ=255)
                         tAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCg                                                                <  2:4591155/93‑1 (MQ=255)
                         tAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCg                                                                >  2:2437445/1‑93 (MQ=255)
                            aTTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGccc                                                             <  1:2960446/93‑1 (MQ=255)
                               cACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCAt                                                          <  2:1918896/93‑1 (MQ=255)
                                 ccGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCAttt                                                        <  2:1653598/93‑1 (MQ=255)
                                  cGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTg                                                       >  1:251185/1‑93 (MQ=255)
                                   gATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGt                                                      >  1:936087/1‑93 (MQ=255)
                                    aTCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTa                                                     <  1:2981818/93‑1 (MQ=255)
                                      cGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAtt                                                   >  1:4495777/1‑93 (MQ=255)
                                      cGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAtt                                                   <  1:925516/93‑1 (MQ=255)
                                       gTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttg                                                  >  2:3714987/1‑93 (MQ=255)
                                       gTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttg                                                  >  2:1788391/1‑93 (MQ=255)
                                        tCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttgt                                                 >  1:1552403/1‑93 (MQ=255)
                                        tCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttgt                                                 >  2:2754862/1‑93 (MQ=255)
                                         cacaTCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttgtt                                                <  1:2453861/93‑1 (MQ=255)
                                         cacaTCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttgtt                                                <  2:211858/93‑1 (MQ=255)
                                          acaTCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttgttg                                               >  2:3324259/1‑93 (MQ=255)
                                          acaTCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAttgttg                                               >  2:1688068/1‑93 (MQ=255)
                                           caTCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGa                                              <  2:2839202/93‑1 (MQ=255)
                                           caTCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGa                                              <  1:1156148/93‑1 (MQ=255)
                                           caTCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGa                                              <  1:3288834/93‑1 (MQ=255)
                                                 ggCCGCTTTTCGCCTGTAGCCAGGCTCCGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCaa                                        <  2:4105896/93‑1 (MQ=255)
                                                  gCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAAc                                       >  2:3081249/1‑93 (MQ=255)
                                                    cGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCt                                     >  2:344404/1‑93 (MQ=255)
                                                     gCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTg                                    >  2:4072768/1‑93 (MQ=255)
                                                     gCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTg                                    >  1:3365123/1‑93 (MQ=255)
                                                     gCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTg                                    <  2:3433454/93‑1 (MQ=255)
                                                      cTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGc                                   <  2:771141/93‑1 (MQ=255)
                                                      cTTTTCGACTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGc                                   >  1:2706891/1‑93 (MQ=255)
                                                       ttttCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCa                                  >  2:1655928/1‑93 (MQ=255)
                                                       ttttCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCa                                  <  2:525471/93‑1 (MQ=255)
                                                         ttCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCAtt                                <  2:2019753/93‑1 (MQ=255)
                                                           cGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTcc                              <  1:1350833/93‑1 (MQ=255)
                                                           cGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTcc                              <  1:3511380/93‑1 (MQ=255)
                                                            gCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTATATTGTTGATGGCAACCTGCATTCCg                             <  1:2730950/93‑1 (MQ=255)
                                                             ccTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGt                            <  2:1867729/93‑1 (MQ=255)
                                                             ccTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGt                            >  2:1223325/1‑93 (MQ=255)
                                                               tGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCg                          >  2:4697057/1‑93 (MQ=255)
                                                               tGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCg                          <  2:4743729/93‑1 (MQ=255)
                                                               tGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCg                          >  2:2209319/1‑93 (MQ=255)
                                                                gTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCgg                         <  1:2994145/93‑1 (MQ=255)
                                                                   gCCAGGCTCTGAAATTGATATCGCTATCGCCATGGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGttt                      <  2:3190337/93‑1 (MQ=255)
                                                                      aGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGtttttt                   >  1:1799295/1‑93 (MQ=255)
                                                                       ggCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTc                  >  1:3580937/1‑93 (MQ=255)
                                                                       ggCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTc                  <  2:1556656/93‑1 (MQ=255)
                                                                        gCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTcc                 <  1:4390469/93‑1 (MQ=255)
                                                                         ctctGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCa                <  1:3263941/93‑1 (MQ=255)
                                                                          tctGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCAt               >  2:1722684/1‑93 (MQ=255)
                                                                          tctGAAATGGATATCGCTATTGCCATTGGCTAGCGCAAAGGTCGCGCCATTTGTATTGTTGATTTCAACCTGCATTCCGTCGGTTTTTTCCAt               >  2:3543477/1‑93 (MQ=255)
                                                                            tGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGc             >  2:2482066/1‑93 (MQ=255)
                                                                                aTTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAAt         >  1:937592/1‑93 (MQ=255)
                                                                                 ttGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAAtt        >  1:107388/1‑93 (MQ=255)
                                                                                   gATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTcc      <  1:2298971/93‑1 (MQ=255)
                                                                                    atatCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTccc     >  1:1550007/1‑93 (MQ=255)
                                                                                    atatCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTccc     <  2:1080029/93‑1 (MQ=255)
                                                                                     tatCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCa    <  1:4481937/93‑1 (MQ=255)
                                                                                      atCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCAc   >  2:2922466/1‑93 (MQ=255)
                                                                                      atCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCAc   <  1:2117811/93‑1 (MQ=255)
                                                                                       tCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACt  >  1:3738481/1‑93 (MQ=255)
                                                                                       tCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACt  >  1:4632366/1‑93 (MQ=255)
                                                                                       tCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACt  <  1:1506279/93‑1 (MQ=255)
                                                                                          |                                                                                         
TCAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACT  >  NC_000913/2448081‑2448260

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ 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
CAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACTCCCT  >  NC_000913/2448082‑2448264
                                                                                         |                                                                                             
CAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATT                                                                                           >  SRR4302544.2092874/1‑94 (MQ=60)
CAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATT                                                                                           >  SRR4302544.4077636/1‑94 (MQ=60)
  AATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGC                                                                                         >  SRR4302544.2862354/1‑94 (MQ=60)
   AGGTCTCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCC                                                                                        <  SRR4302544.3121346/94‑1 (MQ=60)
    TGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCA                                                                                       >  SRR4302544.4598765/1‑94 (MQ=60)
    TGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCA                                                                                       >  SRR4302544.659178/1‑94 (MQ=60)
     GTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCAT                                                                                      >  SRR4302544.4741719/1‑94 (MQ=60)
      TCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATT                                                                                     <  SRR4302544.3272890/94‑1 (MQ=60)
       CGCAGTTAAAGAGGCGGTAAATTCTCCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTG                                                                                    <  SRR4302544.684706/94‑1 (MQ=60)
       CGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTGGCCAGGCTCTGAAATTGATATCGCTATTGCCATTG                                                                                    >  SRR4302544.4137896/1‑94 (MQ=60)
       CGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTG                                                                                    >  SRR4302544.1729228/1‑94 (MQ=60)
        GCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGG                                                                                   <  SRR4302544.1009687/94‑1 (MQ=60)
          AGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTT                                                                                 >  SRR4302544.3762828/1‑94 (MQ=60)
            TTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAG                                                                               <  SRR4302544.3715237/94‑1 (MQ=60)
            TTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAG                                                                               <  SRR4302544.3545372/94‑1 (MQ=60)
              AAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCG                                                                             <  SRR4302544.532684/94‑1 (MQ=60)
                AGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCA                                                                           <  SRR4302544.2323548/94‑1 (MQ=60)
                AGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATGGCCATTGGTTAGCGCA                                                                           <  SRR4302544.2323549/94‑1 (MQ=60)
                   GGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAG                                                                        >  SRR4302544.3405306/1‑94 (MQ=60)
                     CGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGT                                                                      >  SRR4302544.3120036/1‑94 (MQ=60)
                      GGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTC                                                                     <  SRR4302544.2691231/94‑1 (MQ=60)
                      GGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTC                                                                     <  SRR4302544.1819817/94‑1 (MQ=60)
                      GGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTC                                                                     <  SRR4302544.2223659/94‑1 (MQ=60)
                       GTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGACATTGGTTAGCGCAAAGGACG                                                                    >  SRR4302544.3139561/1‑94 (MQ=60)
                           ATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCC                                                                <  SRR4302544.2988864/94‑1 (MQ=60)
                              CACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATT                                                             >  SRR4302544.3470673/1‑94 (MQ=60)
                              CACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATT                                                             >  SRR4302544.697639/1‑94 (MQ=60)
                                CCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTG                                                           >  SRR4302544.253826/1‑94 (MQ=60)
                                 CGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGT                                                          >  SRR4302544.945311/1‑94 (MQ=60)
                                   ATCGTCACATCACGGCCTCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAT                                                        <  SRR4302544.1378405/94‑1 (MQ=60)
                                   ATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTAT                                                        <  SRR4302544.3010415/94‑1 (MQ=60)
                                    TCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATT                                                       >  SRR4302544.4540900/1‑94 (MQ=60)
                                     CGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTG                                                      <  SRR4302544.934637/94‑1 (MQ=60)
                                      GTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGT                                                     >  SRR4302544.1567856/1‑94 (MQ=60)
                                        CACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTG                                                   >  SRR4302544.2880787/1‑94 (MQ=60)
                                        CACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTG                                                   <  SRR4302544.2477268/94‑1 (MQ=60)
                                          CATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGAT                                                 >  SRR4302544.1601529/1‑94 (MQ=60)
                                          CATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGAT                                                 <  SRR4302544.3320172/94‑1 (MQ=60)
                                          CATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGAT                                                 <  SRR4302544.1167653/94‑1 (MQ=60)
                                           ATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATG                                                <  SRR4302544.3882631/94‑1 (MQ=60)
                                            TCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGG                                               <  SRR4302544.977615/94‑1 (MQ=60)
                                             CACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGC                                              <  SRR4302544.4166835/94‑1 (MQ=60)
                                             CACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGC                                              <  SRR4302544.3052839/94‑1 (MQ=60)
                                                   CGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTG                                        >  SRR4302544.3397245/1‑94 (MQ=60)
                                                   CGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTG                                        >  SRR4302544.3397246/1‑94 (MQ=60)
                                                    GCTTTTCGACTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGC                                       >  SRR4302544.2732914/1‑94 (MQ=60)
                                                     CTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCA                                      >  SRR4302544.3636265/1‑94 (MQ=60)
                                                          CGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCG                                 <  SRR4302544.3544880/94‑1 (MQ=60)
                                                          CGCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCG                                 <  SRR4302544.1364222/94‑1 (MQ=60)
                                                           GCCTGTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTATATTGTTGATGGCAACCTGCATTCCGT                                <  SRR4302544.2757272/94‑1 (MQ=60)
                                                               GTAGCCAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGT                            <  SRR4302544.3022851/94‑1 (MQ=60)
                                                                    CAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTT                       >  SRR4302544.4563172/1‑94 (MQ=60)
                                                                    CAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTT                       >  SRR4302544.1816882/1‑94 (MQ=60)
                                                                    CAGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTT                       <  SRR4302544.2320352/94‑1 (MQ=60)
                                                                     AGGCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTC                      >  SRR4302544.3615290/1‑94 (MQ=60)
                                                                       GCTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCA                    <  SRR4302544.4434897/94‑1 (MQ=60)
                                                                        CTCTGAAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCAT                   <  SRR4302544.3295109/94‑1 (MQ=60)
                                                                             AAATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGACAA              >  SRR4302544.2338023/1‑94 (MQ=60)
                                                                              AATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAAT             >  SRR4302544.946828/1‑94 (MQ=60)
                                                                               ATTGATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATT            >  SRR4302544.108370/1‑94 (MQ=60)
                                                                                  GATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCC         <  SRR4302544.2320667/94‑1 (MQ=60)
                                                                                  GATATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCC         >  SRR4302544.1565445/1‑94 (MQ=60)
                                                                                    TATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCAC       <  SRR4302544.4526983/94‑1 (MQ=60)
                                                                                     ATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACT      >  SRR4302544.4678280/1‑94 (MQ=60)
                                                                                     ATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACT      <  SRR4302544.2138149/94‑1 (MQ=60)
                                                                                     ATCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACT      >  SRR4302544.3775640/1‑94 (MQ=60)
                                                                                      TCGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACTC     <  SRR4302544.1521205/94‑1 (MQ=60)
                                                                                       CGCTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACTCC    >  SRR4302544.3630846/1‑94 (MQ=60)
                                                                                         CTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACTCCCT  >  SRR4302544.2069895/1‑94 (MQ=60)
                                                                                         |                                                                                             
CAAATGTCGCAGTTAAAGAGGCGGTAAATTCACCGATCGTCACATCACGGCCGCTTTTCGCCTGTAGCCAGGCTCTGAAATTGATATCGTTATTGCCATTGGTTAGCGCAAAGGTCGCCCCATTTGTATTGTTGATGGCAACCTGCATTCCGTCGGTTTTTTCCATGCCAATTCCCACTCCCT  >  NC_000913/2448082‑2448264

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