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,349,989 T→C M136V (ATG→GTG)  glpT ← sn‑glycerol‑3‑phosphate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,349,9890TC96.7% 273.1 / ‑1.4 92M136V (ATG→GTG) glpTsn‑glycerol‑3‑phosphate permease
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (41/48);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 102 total aligned reads displayed.

GTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349900‑2350078
                                                                                         |                                                                                         
gtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCaccccacccc                                                                                        <  1:245215/93‑1 (MQ=255)
  gcgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCt                                                                                      <  1:2403267/93‑1 (MQ=255)
   cgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTg                                                                                     <  1:568290/93‑1 (MQ=255)
    gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTgg                                                                                    <  2:4012100/93‑1 (MQ=255)
    gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTgg                                                                                    <  1:701989/93‑1 (MQ=255)
    cctcaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTgg                                                                                    <  2:4012101/90‑1 (MQ=255)
     cacaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATCGTACGACCACACGGCGGCCACCCCACCCCCTGGa                                                                                   >  2:3173645/1‑93 (MQ=255)
     cacaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGa                                                                                   >  2:140399/1‑93 (MQ=255)
      acaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaa                                                                                  >  2:3966244/1‑93 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATATTACGACCACACGGCGGCCACCCCACCCCCTGGaccc                                                                                <  2:4724705/93‑1 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacc                                                                                <  1:3555160/93‑1 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACAACCACACGGCGGCCACCCCACCCCCTGTaacc                                                                                >  2:1609351/1‑93 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTGTCTGCGACCACCAGTGCACCATAGTAAGACCACACGGCGGCCACCCCACCCCCTCGaaca                                                                                >  2:1609350/1‑92 (MQ=255)
         gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacca                                                                               >  1:2846031/1‑93 (MQ=255)
          ttCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaa                                                                              >  1:2786213/1‑93 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaaccaac                                                                             <  1:4505283/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:163341/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:2076314/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:4167550/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:3657040/93‑1 (MQ=255)
            ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaacc                                                                            >  2:3147025/1‑93 (MQ=255)
            ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaacc                                                                            >  2:4445890/1‑93 (MQ=255)
             cacaCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCg                                                                           >  2:2785098/1‑93 (MQ=255)
             cacaCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCg                                                                           >  2:2785097/1‑93 (MQ=255)
             cacaCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCg                                                                           <  1:2717054/93‑1 (MQ=255)
                acTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCag                                                                        <  2:2211268/93‑1 (MQ=255)
                 cTGACACAATGCCGCCTCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGCCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       <  2:2052414/93‑1 (MQ=255)
                 cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       >  2:4451438/1‑93 (MQ=255)
                  tGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCagag                                                                      <  2:3482719/93‑1 (MQ=255)
                  cgACACACTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCagag                                                                      <  1:2001376/92‑1 (MQ=255)
                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAgg                                                                     >  2:2918886/1‑93 (MQ=255)
                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAgg                                                                     <  2:4468791/93‑1 (MQ=255)
                    acacAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGTCCACACGGCGGGCACCCCACCCCCTGGAGCCAACCGCAGAGGa                                                                    >  2:3999286/1‑93 (MQ=255)
                       caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaca                                                                 <  1:4055752/93‑1 (MQ=255)
                        aaTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaa                                                                <  2:3119481/93‑1 (MQ=255)
                        aaTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaa                                                                <  1:1191258/93‑1 (MQ=255)
                         aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               >  1:702/1‑93 (MQ=255)
                          tGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAacaaca                                                              <  2:3361950/93‑1 (MQ=255)
                          agccgccCCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAacaaca                                                              <  1:2533412/92‑1 (MQ=255)
                           gccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAg                                                             >  2:2693309/1‑93 (MQ=255)
                           gccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAg                                                             <  1:4123762/93‑1 (MQ=255)
                            ccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGt                                                            <  1:2398869/93‑1 (MQ=255)
                             cgccCCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           <  2:2088840/93‑1 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           <  1:3496573/93‑1 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           <  1:797612/93‑1 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           >  1:4576877/1‑93 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           >  1:794222/1‑93 (MQ=255)
                              gccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                                                          >  1:3328362/1‑93 (MQ=255)
                               ccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACa                                                         >  1:605596/1‑93 (MQ=255)
                                   gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     <  2:2884869/93‑1 (MQ=255)
                                   gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     <  1:3642239/93‑1 (MQ=255)
                                    ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACac                                                    >  1:1213642/1‑92 (MQ=255)
                                    ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt                                                    >  2:1782126/1‑93 (MQ=255)
                                      ctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCa                                                  <  2:1762331/93‑1 (MQ=255)
                                        ttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACGGTACAAACATCAcc                                                <  2:1576562/93‑1 (MQ=255)
                                         tctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCg                                               >  1:4034362/1‑93 (MQ=255)
                                           tGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCa                                             >  2:2993084/1‑93 (MQ=255)
                                             cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           <  2:3777638/93‑1 (MQ=255)
                                              gACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATg                                          >  1:3037497/1‑93 (MQ=255)
                                              gACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATg                                          <  2:3600278/93‑1 (MQ=255)
                                               accaccAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGc                                         >  2:1634956/1‑93 (MQ=255)
                                                 caccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc                                       <  2:4043672/93‑1 (MQ=255)
                                                 caccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc                                       >  1:3016525/1‑93 (MQ=255)
                                                  cccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      <  2:1331511/92‑1 (MQ=255)
                                                  accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      >  2:4449900/1‑93 (MQ=255)
                                                   ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACGGTACAAACATCACCGCAATGCTCGa                                     <  1:3989620/93‑1 (MQ=255)
                                                   ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     >  1:617919/1‑93 (MQ=255)
                                                    caGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGAc                                    <  1:3720404/93‑1 (MQ=255)
                                                     aGTGCCCCATAGTACGACCTCACGGCGGCCCCCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACg                                   <  1:1053662/93‑1 (MQ=255)
                                                      gTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACgc                                  >  2:756399/1‑92 (MQ=255)
                                                      gTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGt                                  <  1:1806028/93‑1 (MQ=255)
                                                       ggCACCATAGTACGACCACACGGCGGTCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTc                                 <  1:2566135/92‑1 (MQ=255)
                                                        gCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCg                                >  1:1330551/1‑93 (MQ=255)
                                                          aCCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGcc                              <  1:2037548/93‑1 (MQ=255)
                                                            cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  2:3585953/93‑1 (MQ=255)
                                                            cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  1:2621701/93‑1 (MQ=255)
                                                            cATAGTACCACCCCACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  1:2621702/93‑1 (MQ=255)
                                                             aTAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAt                           >  1:94299/1‑93 (MQ=255)
                                                                gTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGc                        >  1:2016470/1‑93 (MQ=255)
                                                                   cGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGcaca                     >  2:1947514/1‑93 (MQ=255)
                                                                   cGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGcaaa                     >  2:3584943/1‑91 (MQ=255)
                                                                     aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   <  1:2638661/93‑1 (MQ=255)
                                                                     aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   >  2:2841016/1‑93 (MQ=255)
                                                                       cacaCGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGc                 >  1:106464/1‑93 (MQ=255)
                                                                        acacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGcc                >  1:1345527/1‑93 (MQ=255)
                                                                         cacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGccc               <  1:2908388/93‑1 (MQ=255)
                                                                         cacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGccc               <  1:1138387/93‑1 (MQ=255)
                                                                          acGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              >  2:168438/1‑93 (MQ=255)
                                                                            ggcggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATa            <  2:909352/93‑1 (MQ=255)
                                                                            ggcggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATa            <  2:1728443/93‑1 (MQ=255)
                                                                             gcggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATaa           <  1:259638/93‑1 (MQ=255)
                                                                                  caccccacccCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaac      <  2:428050/93‑1 (MQ=255)
                                                                                  cACCCCATCCGCTGGAACCAACCGCAAAAGAACAACAGTACAAACATCACCCCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaac      <  2:1578314/93‑1 (MQ=255)
                                                                                   accccaccccCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     >  2:1283587/1‑93 (MQ=255)
                                                                                    ccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACAt    <  1:4743703/93‑1 (MQ=255)
                                                                                     cccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   >  2:2892827/1‑93 (MQ=255)
                                                                                     cccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   >  2:2388388/1‑93 (MQ=255)
                                                                                     cccACCCCCTGGAACCAACCGCAGAGGAACAACAGGACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   >  2:3700819/1‑93 (MQ=255)
                                                                                      ccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCa  >  2:3473540/1‑93 (MQ=255)
                                                                                      ccACCCCATGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATCCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCa  <  2:2705779/93‑1 (MQ=255)
                                                                                         |                                                                                         
GTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349900‑2350078

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCACTGC  >  NC_000913/2349896‑2350082
                                                                                             |                                                                                             
CGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCAC                                                                                               >  SRR4302544.2222462/1‑94 (MQ=60)
 GTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACC                                                                                              <  SRR4302544.2245635/94‑1 (MQ=60)
  TTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCC                                                                                             <  SRR4302544.2240868/94‑1 (MQ=60)
    GTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCC                                                                                           <  SRR4302544.247779/94‑1 (MQ=60)
      GCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTG                                                                                         <  SRR4302544.2426075/94‑1 (MQ=60)
       CGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGG                                                                                        <  SRR4302544.573841/94‑1 (MQ=60)
        GCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGA                                                                                       <  SRR4302544.708978/94‑1 (MQ=60)
            AGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCA                                                                                   <  SRR4302544.3589229/94‑1 (MQ=60)
            AGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCA                                                                                   >  SRR4302544.2873521/1‑94 (MQ=60)
             GTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAA                                                                                  >  SRR4302544.2813026/1‑94 (MQ=60)
               TCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACC                                                                                <  SRR4302544.4550466/94‑1 (MQ=60)
               TCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACC                                                                                <  SRR4302544.3693525/94‑1 (MQ=60)
               TCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACC                                                                                <  SRR4302544.4208839/94‑1 (MQ=60)
               TCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACC                                                                                <  SRR4302544.2096267/94‑1 (MQ=60)
               TCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACC                                                                                <  SRR4302544.164937/94‑1 (MQ=60)
                 CACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGC                                                                              <  SRR4302544.2743204/94‑1 (MQ=60)
                      CGACACACTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGG                                                                         <  SRR4302544.2020810/94‑1 (MQ=60)
                           CAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACAGGGCGGCAACCCCACCCCCTGGAACCAccagacgagagaaaa                                                                    >  SRR4302544.1092672/1‑79 (MQ=60)
                           CAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAA                                                                    <  SRR4302544.4096147/94‑1 (MQ=60)
                            AATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAAC                                                                   >  SRR4302544.710/1‑94 (MQ=60)
                            AATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAAC                                                                   <  SRR4302544.1203078/94‑1 (MQ=60)
                              AGCCGCCCCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAG                                                                 <  SRR4302544.2557749/94‑1 (MQ=60)
                               GCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGT                                                                <  SRR4302544.4164679/94‑1 (MQ=60)
                                CCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTA                                                               <  SRR4302544.2421618/94‑1 (MQ=60)
                                CCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTA                                                               >  SRR4302544.4622481/1‑94 (MQ=60)
                                CCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTA                                                               >  SRR4302544.801907/1‑94 (MQ=60)
                                 CGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTAC                                                              <  SRR4302544.3529913/94‑1 (MQ=60)
                                 CGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAC                                                              >  SRR4302544.3360059/1‑94 (MQ=60)
                                 CGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAC                                                              <  SRR4302544.805351/94‑1 (MQ=60)
                                  GCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACA                                                             >  SRR4302544.611540/1‑94 (MQ=60)
                                       GTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAT                                                        <  SRR4302544.3678578/94‑1 (MQ=60)
                                       GTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAC                                                        >  SRR4302544.1225656/1‑94 (MQ=60)
                                            TTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCG                                                   >  SRR4302544.4074549/1‑94 (MQ=60)
                                                 CGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATG                                              >  SRR4302544.3066571/1‑94 (MQ=60)
                                                    CCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC                                           >  SRR4302544.3045445/1‑94 (MQ=60)
                                                      ACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGA                                         >  SRR4302544.623987/1‑94 (MQ=60)
                                                       CCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACGGTACAAACATCACCGCAATGCTCGAC                                        <  SRR4302544.4029335/94‑1 (MQ=60)
                                                        CAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACG                                       <  SRR4302544.3757409/94‑1 (MQ=60)
                                                         AGTGCCCCATAGTACGACCTCACGGCGGCCCCCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGT                                      <  SRR4302544.1064106/94‑1 (MQ=60)
                                                          GTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTC                                     <  SRR4302544.1823678/94‑1 (MQ=60)
                                                           TGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCG                                    >  SRR4302544.1343644/1‑94 (MQ=60)
                                                           GGCACCATAGTACGACCACACGGCGGTCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCG                                    <  SRR4302544.2590789/94‑1 (MQ=60)
                                                              ACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCC                                 <  SRR4302544.2057230/94‑1 (MQ=60)
                                                                CATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT                               >  SRR4302544.95169/1‑94 (MQ=60)
                                                                CATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT                               <  SRR4302544.2646845/94‑1 (MQ=60)
                                                                CATAGTACCACCCCACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAT                               <  SRR4302544.2646846/94‑1 (MQ=60)
                                                                   AGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGC                            >  SRR4302544.2036057/1‑94 (MQ=60)
                                                                         ACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAG                      <  SRR4302544.2663954/94‑1 (MQ=60)
                                                                          CCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGC                     >  SRR4302544.107438/1‑94 (MQ=60)
                                                                           CACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCC                    >  SRR4302544.1358843/1‑94 (MQ=60)
                                                                             CACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCA                  <  SRR4302544.2936373/94‑1 (MQ=60)
                                                                             CACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCA                  <  SRR4302544.1149707/94‑1 (MQ=60)
                                                                                 GCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAA              <  SRR4302544.262385/94‑1 (MQ=60)
                                                                                        CCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATC       <  SRR4302544.4790451/94‑1 (MQ=60)
                                                                                           CACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCACT    <  SRR4302544.3961589/94‑1 (MQ=60)
                                                                                             CCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCACTGC  >  SRR4302544.4514520/1‑94 (MQ=60)
                                                                                             CCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCACTGC  <  SRR4302544.1083790/94‑1 (MQ=60)
                                                                                             |                                                                                             
CGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCACTGC  >  NC_000913/2349896‑2350082

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