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 4,225,139 T→C intergenic (+49/‑171) metH → / → yjbB B12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,225,1390TC92.6% 253.4 / 9.0 95intergenic (+49/‑171)metH/yjbBB12‑dependent homocysteine‑N5‑methyltetrahydrofolate transmethylase/putative alpha helix protein
Reads supporting (aligned to +/- strand):  ref base T (5/2);  new base C (46/42);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 4.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

Only 100 of 113 total aligned reads displayed.

GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225049‑4225226
                                                                                          |                                                                                        
gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgcgc                                                                                        >  1:1036598/3‑93 (MQ=255)
gggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAgcgc                                                                                        >  1:1619385/3‑93 (MQ=255)
 ggcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGcc                                                                                       <  2:69503/92‑1 (MQ=255)
  gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGccc                                                                                      <  1:3501655/93‑1 (MQ=255)
  gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGccc                                                                                      <  1:2591310/93‑1 (MQ=255)
  gcgcTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGccc                                                                                      <  2:4554733/93‑1 (MQ=255)
      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  1:4694870/93‑1 (MQ=255)
      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  1:4005787/93‑1 (MQ=255)
      tggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtc                                                                                  <  2:1617809/93‑1 (MQ=255)
       ggctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGtct                                                                                 >  1:2311024/1‑93 (MQ=255)
        gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCTATGAGCGCCCTGtctc                                                                                >  1:4219409/1‑93 (MQ=255)
        gctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCCCCCTGtctc                                                                                <  2:1637749/93‑1 (MQ=255)
         ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTctct                                                                               >  2:1105845/1‑93 (MQ=255)
         ctggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTctct                                                                               >  2:1228678/1‑93 (MQ=255)
          tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtt                                                                              <  2:1673729/93‑1 (MQ=255)
          tggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtt                                                                              <  1:2499506/93‑1 (MQ=255)
           ggcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCttt                                                                             >  2:2103401/1‑93 (MQ=255)
           ggcACCGAAGCTGGGGTGTGACGCGGACTGATTCACAAATCTGTCATTTTGCATTACAACAAACAGGGCGTTCAATGAGCGCCCTGTCTCttt                                                                             <  2:1700663/93‑1 (MQ=255)
            gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTtta                                                                            <  1:740660/93‑1 (MQ=255)
            gcACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCtttt                                                                            >  2:2747760/1‑92 (MQ=255)
              aCCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTttatt                                                                          <  1:1539735/93‑1 (MQ=255)
                  aaTCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGCGCTCTGTCTCTTTATTtata                                                                      >  1:4220022/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  1:1245056/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  2:3284300/1‑93 (MQ=255)
                      tGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaa                                                                  >  2:1528494/1‑93 (MQ=255)
                       ggggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGaaa                                                                 <  2:3762880/93‑1 (MQ=255)
                        gggTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAac                                                                <  1:3311950/93‑1 (MQ=255)
                          gTATGGCGGGAACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCTCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAacac                                                              <  2:2374154/93‑1 (MQ=255)
                           tATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACt                                                             <  1:2302347/93‑1 (MQ=255)
                            aTGACGAGCACTTATTCAAAACTCTAACACTTTACTTAACAACAACCGGGGCCCGCAAGGGGCGCCCTGTCTCTTTATTAATATGAAACACtt                                                            <  2:3591747/93‑1 (MQ=255)
                             tGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGACCTGTCTCTTTATTAATATGAAACACTta                                                           >  1:1574009/1‑93 (MQ=255)
                               aCGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTtata                                                         <  1:3301162/93‑1 (MQ=255)
                                 gcgGACTGATCCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACt                                                       <  2:390043/93‑1 (MQ=255)
                                  cgGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTg                                                      <  1:549841/93‑1 (MQ=255)
                                   ggACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTgg                                                     >  2:4460082/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:3911562/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGa                                                    >  1:3553357/1‑93 (MQ=255)
                                    gACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCCCAAAGAGCGCCCTGGCCCCCTTTTAATAAGAAAAACATacacaaga                                                    >  1:4204617/1‑89 (MQ=255)
                                      cTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCCGTCTCTTTTTTAATATGAAACACTTATACTGGaaa                                                  >  2:142943/1‑93 (MQ=255)
                                       tGATTCACAAATCTGTCACTTTCCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAc                                                 <  2:3000715/93‑1 (MQ=255)
                                        gATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACa                                                >  1:4572993/1‑93 (MQ=255)
                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                              >  1:2198596/1‑93 (MQ=255)
                                          ttCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAgg                                              >  1:2623059/1‑93 (MQ=255)
                                           tCAAAAATCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGc                                             <  1:1154871/93‑1 (MQ=255)
                                             acaAATCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                          <  1:1336055/93‑1 (MQ=255)
                                             acaAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTg                                          <  2:4554091/93‑1 (MQ=255)
                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTgg                                         <  1:636111/93‑1 (MQ=255)
                                              caAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAAAAGGCTgg                                         >  2:1271649/1‑93 (MQ=255)
                                               aaaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTCATATGAAACACTTATACTGGAAACAGGCTGGa                                        >  1:2967095/1‑93 (MQ=255)
                                                aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                       >  1:314942/1‑93 (MQ=255)
                                                aaTCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaa                                       <  2:4047678/93‑1 (MQ=255)
                                                aaTCCGCCACTTTTCCCTACAAAAAACGGGGCGCTCAACAATCCCCCTGTCTCTTTATTAATATGAAACCCTTATACTGGAAACAGGCTGGaa                                       <  1:1279865/93‑1 (MQ=255)
                                                 agcTGTCCTTTTTCCTAACAACAAAAGGGGCATGAAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaat                                      <  2:1303997/91‑1 (MQ=255)
                                                  tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                     >  2:2459906/1‑93 (MQ=255)
                                                  tCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaata                                     >  1:1156548/1‑93 (MQ=255)
                                                    tGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGaataat                                   >  1:489683/1‑93 (MQ=255)
                                                     gTCACTTTTCCTTACAACAAAGAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                  <  2:441335/93‑1 (MQ=255)
                                                     gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                  >  1:1661322/1‑93 (MQ=255)
                                                     gTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATc                                  >  2:2816590/1‑93 (MQ=255)
                                                      tCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                 <  2:1562804/93‑1 (MQ=255)
                                                      tCACTGTTCCTTACAACAAACAGGGCGCTCAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCt                                 <  1:746581/93‑1 (MQ=255)
                                                       cACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCtt                                >  2:2768917/1‑93 (MQ=255)
                                                        ccTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTa                               <  2:1609870/92‑1 (MQ=255)
                                                          ttttCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             <  2:2004668/93‑1 (MQ=255)
                                                          ttttCCTAACAACAAACAGGGCGCTAAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGc                             <  1:2035410/93‑1 (MQ=255)
                                                             tCCTTACAACAAACAGGGGGCTCCAAGAGGGCCCCGTCTCTTTTTTTATATTTAAAACTTATAATCGGAACAGGCTGGGAAAAACTTTAGCCg                           >  2:1876739/1‑93 (MQ=255)
                                                             tCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                          <  1:3971618/93‑1 (MQ=255)
                                                             tCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCgg                          <  2:3851146/93‑1 (MQ=255)
                                                              ccTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAAGCTGGAATAATATTAGCCGGa                         >  2:458164/1‑93 (MQ=255)
                                                               cTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaa                        >  1:1975538/1‑93 (MQ=255)
                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGCCGTaaa                       >  1:527562/1‑93 (MQ=255)
                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                       >  2:3142556/1‑93 (MQ=255)
                                                                ttACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaa                       <  2:2788553/93‑1 (MQ=255)
                                                                  acaacaAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGaaaaa                     >  1:4085365/1‑93 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGCAAAACg                   >  2:235014/1‑93 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   <  1:4424948/93‑1 (MQ=255)
                                                                    aacaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   >  1:4330243/1‑93 (MQ=255)
                                                                    aaaaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACg                   <  1:3012164/90‑1 (MQ=255)
                                                                     aaaaacgGGGCGCTAAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                  <  2:3447524/91‑1 (MQ=255)
                                                                     aaaaacgGGGCGCTAAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGa                  <  1:4172137/91‑1 (MQ=255)
                                                                      caaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                 >  1:2389164/1‑93 (MQ=255)
                                                                      aaaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGaa                 <  1:2233497/92‑1 (MQ=255)
                                                                       aaacaGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAt                >  1:3123141/1‑93 (MQ=255)
                                                                        agcggggCACTCAGGGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGAAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATg               <  2:4215167/89‑1 (MQ=255)
                                                                        aacaGGGCGCTCAATGAGCGCCCTGTCTCTTTTTTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAAGAATg               >  1:1959440/1‑93 (MQ=255)
                                                                         acggggCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGt              <  1:309447/90‑1 (MQ=255)
                                                                          caGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATATc             >  2:4006359/1‑93 (MQ=255)
                                                                          caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTc             <  1:2694478/93‑1 (MQ=255)
                                                                          caGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAAGAATGTc             >  1:4545064/1‑93 (MQ=255)
                                                                           aGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCa            >  1:2909787/1‑93 (MQ=255)
                                                                            gggCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaa           >  2:1052547/1‑93 (MQ=255)
                                                                              gcgcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCaaaa         <  2:2224376/93‑1 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  1:1038196/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  2:3801448/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  1:1320567/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       >  2:1592270/1‑93 (MQ=255)
                                                                                gcTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAgg       <  1:1512010/93‑1 (MQ=255)
                                                                                  tCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAggga     >  1:945388/1‑91 (MQ=255)
                                                                                     aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg  >  2:1115223/1‑93 (MQ=255)
                                                                                     aTGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAg  >  2:2916311/1‑93 (MQ=255)
                                                                                          |                                                                                        
GAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAG  >  NC_000913/4225049‑4225226

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGA  >  NC_000913/4225046‑4225232
                                                                                             |                                                                                             
GTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGC                                                                                               <  SRR4302544.2472985/94‑1 (MQ=60)
GTTGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGC                                                                                               <  SRR4302544.4290140/94‑1 (MQ=60)
  TGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGC                                                                                             >  SRR4302544.1046904/1‑94 (MQ=60)
  TGGGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGC                                                                                             >  SRR4302544.1635338/1‑94 (MQ=60)
     GCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCT                                                                                          <  SRR4302544.2616234/94‑1 (MQ=60)
     GCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCT                                                                                          <  SRR4302544.3208580/94‑1 (MQ=60)
     GCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCT                                                                                          <  SRR4302544.3535061/94‑1 (MQ=60)
      CGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCATTACAACAAACAGGGCGCTCAATGAGCGCCCTG                                                                                         <  SRR4302544.94162/94‑1 (MQ=60)
         TGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCT                                                                                      <  SRR4302544.4741228/94‑1 (MQ=60)
         TGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCT                                                                                      <  SRR4302544.4045649/94‑1 (MQ=60)
         TGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCT                                                                                      >  SRR4302544.1387695/1‑94 (MQ=60)
         TGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCT                                                                                      >  SRR4302544.2332875/1‑94 (MQ=60)
          GGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCTATGAGCGCCCTGTCTC                                                                                     >  SRR4302544.4262130/1‑94 (MQ=60)
             TGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTT                                                                                  <  SRR4302544.2523467/94‑1 (MQ=60)
               GCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTAT                                                                                <  SRR4302544.747912/94‑1 (MQ=60)
                 ACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTA                                                                              <  SRR4302544.1555034/94‑1 (MQ=60)
                    GAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGCGCTCTGTCTCTTTATTTATA                                                                           >  SRR4302544.4262760/1‑94 (MQ=60)
                        CTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAA                                                                       >  SRR4302544.1257327/1‑94 (MQ=60)
                           GGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACA                                                                    <  SRR4302544.3343489/94‑1 (MQ=60)
                             GTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACT                                                                  >  SRR4302544.412825/1‑94 (MQ=60)
                              TATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTT                                                                 <  SRR4302544.2324090/94‑1 (MQ=60)
                               ATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGACCTGTCTCTTTATTAATATGAAACACTTA                                                                >  SRR4302544.1589569/1‑94 (MQ=60)
                                TGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTAT                                                               <  SRR4302544.4669221/94‑1 (MQ=60)
                                  ACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATAC                                                             <  SRR4302544.3332604/94‑1 (MQ=60)
                                     CGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGG                                                          <  SRR4302544.555230/94‑1 (MQ=60)
                                      GGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGA                                                         >  SRR4302544.3950420/1‑94 (MQ=60)
                                      GGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGA                                                         >  SRR4302544.3587397/1‑94 (MQ=60)
                                          TGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACA                                                     >  SRR4302544.4618561/1‑94 (MQ=60)
                                          TGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACA                                                     <  SRR4302544.2107651/94‑1 (MQ=60)
                                           GATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCCATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAG                                                    >  SRR4302544.2705867/1‑94 (MQ=60)
                                            ATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGG                                                   >  SRR4302544.2648212/1‑94 (MQ=60)
                                            ATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGG                                                   >  SRR4302544.2219318/1‑94 (MQ=60)
                                              TCAAAAATCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCT                                                 <  SRR4302544.1166359/94‑1 (MQ=60)
                                                ACAAATCTGTCACTTTTCCTTACAACAAACGGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGG                                               <  SRR4302544.1349228/94‑1 (MQ=60)
                                                 CAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTCATATGAAACACTTATACTGGAAACAGGCTGGA                                              >  SRR4302544.2995559/1‑94 (MQ=60)
                                                 CAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGA                                              <  SRR4302544.642383/94‑1 (MQ=60)
                                                  AAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAA                                             >  SRR4302544.318214/1‑94 (MQ=60)
                                                    ATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATA                                           >  SRR4302544.1168059/1‑94 (MQ=60)
                                                      CTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAAT                                         >  SRR4302544.494487/1‑94 (MQ=60)
                                                       TGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATC                                        >  SRR4302544.1677682/1‑94 (MQ=60)
                                                         TCACTGTTCCTTACAACAAACAGGGCGCTCAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTT                                      <  SRR4302544.753897/94‑1 (MQ=60)
                                                             TTTTCCTAACAACAAACAGGGCGCTAAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCC                                  <  SRR4302544.2055083/94‑1 (MQ=60)
                                                                TCCTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGA                               <  SRR4302544.4011125/94‑1 (MQ=60)
                                                                 CCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAA                              >  SRR4302544.1994736/1‑94 (MQ=60)
                                                                  CTTACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGGATAATCTTAGCCGTAAA                             >  SRR4302544.532798/1‑94 (MQ=60)
                                                                    TACAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAA                           >  SRR4302544.4125997/1‑94 (MQ=60)
                                                                      CAACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACG                         >  SRR4302544.4374337/1‑94 (MQ=60)
                                                                       AACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGA                        <  SRR4302544.4469588/94‑1 (MQ=60)
                                                                       AAAAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGA                        <  SRR4302544.3041041/94‑1 (MQ=60)
                                                                        ACAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAA                       >  SRR4302544.2411788/1‑94 (MQ=60)
                                                                        AAAAACGGGGCGCTAAATGAGCTCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAA                       <  SRR4302544.4213463/94‑1 (MQ=60)
                                                                         CAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAT                      >  SRR4302544.3152825/1‑94 (MQ=60)
                                                                         AAAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAAT                      <  SRR4302544.2254503/94‑1 (MQ=60)
                                                                          AAACAGGGCGCTCAATGAGCGCCCTGTCTCTTTTTTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAAGAATG                     >  SRR4302544.1978452/1‑94 (MQ=60)
                                                                            ACGGGGCGCTCAATGAGCCCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTC                   <  SRR4302544.312667/94‑1 (MQ=60)
                                                                            ACAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAAAGAATGTC                   >  SRR4302544.4590469/1‑94 (MQ=60)
                                                                             CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCA                  >  SRR4302544.2937781/1‑94 (MQ=60)
                                                                             CAGGGCGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCA                  <  SRR4302544.2720344/94‑1 (MQ=60)
                                                                                  CGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGG             >  SRR4302544.1048512/1‑94 (MQ=60)
                                                                                  CGCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGG             >  SRR4302544.1333556/1‑94 (MQ=60)
                                                                                   GCTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGC            <  SRR4302544.1526996/94‑1 (MQ=60)
                                                                                    CTCAATGAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGGA           >  SRR4302544.954699/1‑94 (MQ=60)
                                                                                          GAGCGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAA     >  SRR4302544.288703/1‑94 (MQ=60)
                                                                                             CGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGA  >  SRR4302544.4619386/1‑94 (MQ=60)
                                                                                             |                                                                                             
GTTGAGCGCTGGCTGGCACCGAATCTGGGGTATGACGCGGACTGATTCACAAATCTGTCACTTTTCCTTACAACAAACAGGGCGCTCAATGAGTGCCCTGTCTCTTTATTAATATGAAACACTTATACTGGAAACAGGCTGGAATAATCTTAGCCGGAAAAACGAATGTCAAAAGGCACAGTAAAGA  >  NC_000913/4225046‑4225232

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