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 647 270.1 21258240 98.4% 20918108 94.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,041,126 A→G E43G (GAG→GGG)  dsbA → protein disulfide isomerase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,041,1260AG98.9% 2182.7 / 1.5 642E43G (GAG→GGG) dsbAprotein disulfide isomerase I
Reads supporting (aligned to +/- strand):  ref base A (0/7);  new base G (302/332);  total (302/340)
Fisher's exact test for biased strand distribution p-value = 1.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 655 total aligned reads displayed.

TTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGT  >  NC_000913/4041035‑4041216
                                                                                           |                                                                                          
ttagcgtttagcgCATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGt                                                                                          <  1:6565676/94‑1 (MQ=255)
    cgtttagcgCATAGGCGGCGCAGTATGGAGATGGTAAACAGTACACTACCCTGGGAAAACCGGAAGTTTGCGCGCCGCCAGGGCTGGGGttttt                                                                                      >  1:9096448/1‑94 (MQ=255)
       ttagcgCATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTtctt                                                                                   >  2:2222337/1‑93 (MQ=255)
       ttagcgCATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTtctc                                                                                   >  2:10016182/1‑94 (MQ=255)
          gcgcATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCtcttt                                                                                >  2:5020640/1‑94 (MQ=255)
            gcATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCtctttct                                                                              <  2:2017590/94‑1 (MQ=255)
            gcATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCtctttct                                                                              >  1:3250877/1‑94 (MQ=255)
             cATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCtctttctt                                                                             <  2:1257186/94‑1 (MQ=255)
              aTCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGGTGGCGCGCCGCAAGTGCTGGGGCTTTTCTCTttcttc                                                                            >  2:3717207/1‑94 (MQ=255)
              aTCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTttcttc                                                                            >  1:9210708/1‑94 (MQ=255)
               tCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTtcttct                                                                           <  1:8845648/94‑1 (MQ=255)
                cggcggCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTg                                                                          <  2:4155772/94‑1 (MQ=255)
                cggcggCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTg                                                                          <  2:5370495/94‑1 (MQ=255)
                 ggcggcGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGc                                                                         <  1:2304760/94‑1 (MQ=255)
                   cggcgCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGccc                                                                       <  1:9635928/94‑1 (MQ=255)
                   cggcgCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGccc                                                                       <  1:399218/94‑1 (MQ=255)
                    ggcgCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCg                                                                      >  1:1828375/1‑94 (MQ=255)
                       gcAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCAc                                                                   >  2:548988/1‑94 (MQ=255)
                       ccAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCAc                                                                   <  1:4893752/93‑1 (MQ=255)
                         aGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTg                                                                 >  2:7619656/1‑94 (MQ=255)
                         aGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTg                                                                 >  1:6218030/1‑94 (MQ=255)
                         aGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTg                                                                 >  2:6389199/1‑94 (MQ=255)
                          gTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGc                                                                >  1:5099792/1‑94 (MQ=255)
                           tATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCt                                                               <  2:9455570/94‑1 (MQ=255)
                           tATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCt                                                               <  1:8214612/94‑1 (MQ=255)
                            aTGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTa                                                              >  1:10161131/1‑94 (MQ=255)
                             tGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTAt                                                             <  2:7379604/94‑1 (MQ=255)
                              gAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATc                                                            >  1:7199274/1‑94 (MQ=255)
                               aaGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCa                                                           >  1:6901074/1‑94 (MQ=255)
                                aGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAg                                                          <  1:5541410/94‑1 (MQ=255)
                                 aaTGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCGGGGTTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGt                                                         <  1:4634301/93‑1 (MQ=255)
                                  aTGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGtt                                                        <  2:6319959/94‑1 (MQ=255)
                                      gaaaCAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTgaa                                                    <  2:6988030/93‑1 (MQ=255)
                                       aaaCAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTgaag                                                   >  2:613248/1‑94 (MQ=255)
                                        aaCAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTgaaga                                                  <  1:849213/94‑1 (MQ=255)
                                        aaCAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTgaaga                                                  <  2:8673068/94‑1 (MQ=255)
                                        aaCAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTgaaga                                                  >  1:7310180/1‑94 (MQ=255)
                                           aGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGt                                               <  2:4293163/94‑1 (MQ=255)
                                           aGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGt                                               >  2:3284673/1‑94 (MQ=255)
                                             tACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTc                                             >  1:5002377/1‑94 (MQ=255)
                                              acacTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCt                                            <  2:9127745/94‑1 (MQ=255)
                                               cacTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTg                                           <  2:965687/94‑1 (MQ=255)
                                               cacTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTg                                           <  1:10444574/94‑1 (MQ=255)
                                               cacTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTg                                           >  2:8865138/1‑94 (MQ=255)
                                                   aCCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTCTCAGTTTGCAGAAGTTCTGCata                                       >  1:8082264/1‑94 (MQ=255)
                                                    cccTGGAAAAACCGGTAGCTGGCTCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCatat                                      <  1:7553732/94‑1 (MQ=255)
                                                    cccTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCatat                                      <  2:1424607/94‑1 (MQ=255)
                                                    cccTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCatat                                      <  2:2462066/94‑1 (MQ=255)
                                                      cTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATAttt                                    <  2:5912951/94‑1 (MQ=255)
                                                        ggAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCt                                  >  1:8387155/1‑94 (MQ=255)
                                                         gAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTg                                 <  2:1822162/94‑1 (MQ=255)
                                                         gAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTg                                 <  1:6066634/94‑1 (MQ=255)
                                                         gAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTg                                 <  2:5765248/94‑1 (MQ=255)
                                                            aaaCCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATa                              >  2:3861295/1‑94 (MQ=255)
                                                            aaaCCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATa                              <  1:4328305/94‑1 (MQ=255)
                                                              aCCGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAAt                            <  2:7922755/94‑1 (MQ=255)
                                                               ccGGTCACTGGCGCGCAGCTAGTGCTGGGGTTTTTCTCTTTCTTCTGACCGCACTGCTATCAGTTTGAAGAAGTTCAGCATTTTTCTGATAAtg                           <  1:9402227/94‑1 (MQ=255)
                                                                cGGTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAAtgt                          >  2:8551631/1‑94 (MQ=255)
                                                                 ggTAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAAtgtg                         <  2:7137258/94‑1 (MQ=255)
                                                                   tAGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTgaa                       <  1:4789159/94‑1 (MQ=255)
                                                                    aGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTgaag                      >  1:10314548/1‑94 (MQ=255)
                                                                    aGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTgaag                      <  2:9226345/94‑1 (MQ=255)
                                                                    aGCTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTgaag                      >  1:2105837/1‑94 (MQ=255)
                                                                      cTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTgaagaa                    <  1:3715505/94‑1 (MQ=255)
                                                                      cTGGCGCGCCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTgaagaa                    <  2:4924009/94‑1 (MQ=255)
                                                                         gcgcgcCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGaaaaa                 <  1:5333450/94‑1 (MQ=255)
                                                                          cgcgcCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGaaaaaa                <  2:4709717/94‑1 (MQ=255)
                                                                           gcgcCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAAc               >  2:8251163/1‑94 (MQ=255)
                                                                           gcgcCGCAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAAc               <  1:5942530/94‑1 (MQ=255)
                                                                            cgccgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACt              >  1:2506004/1‑94 (MQ=255)
                                                                            cgccgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACt              <  2:307790/94‑1 (MQ=255)
                                                                             gccgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTg             >  1:5395842/1‑94 (MQ=255)
                                                                             gccgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTg             >  1:3493349/1‑94 (MQ=255)
                                                                              ccgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGc            <  1:9836513/94‑1 (MQ=255)
                                                                              ccgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACGGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGGAGAAAAAACTGc            <  1:1734812/94‑1 (MQ=255)
                                                                               cgcAAGTTCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGcc           <  2:10510833/94‑1 (MQ=255)
                                                                               cgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGcg           >  2:3351752/1‑93 (MQ=255)
                                                                               cgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGcc           >  2:9847850/1‑94 (MQ=255)
                                                                               cgcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGcc           >  1:4167916/1‑94 (MQ=255)
                                                                                gcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCg          <  2:6628110/94‑1 (MQ=255)
                                                                                gcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCg          >  2:7976728/1‑94 (MQ=255)
                                                                                gcAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCg          >  1:6431777/1‑94 (MQ=255)
                                                                                 cAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCgg         <  1:3045804/94‑1 (MQ=255)
                                                                                 cAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCgg         <  1:1620866/94‑1 (MQ=255)
                                                                                 cAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCgg         >  1:7815619/1‑94 (MQ=255)
                                                                                 cAAGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCgg         <  1:3874756/94‑1 (MQ=255)
                                                                                  aaGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGa        >  2:452834/1‑94 (MQ=255)
                                                                                  aaGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGa        >  1:2776896/1‑94 (MQ=255)
                                                                                   aGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGaa       >  1:2697735/1‑94 (MQ=255)
                                                                                   aGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGaa       >  2:6820716/1‑94 (MQ=255)
                                                                                   aGTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGaa       <  1:10546312/94‑1 (MQ=255)
                                                                                    gTGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAg      <  1:1005943/94‑1 (MQ=255)
                                                                                     tGCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAgg     <  2:2847187/94‑1 (MQ=255)
                                                                                      gCTGGTGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGc    <  2:8753059/94‑1 (MQ=255)
                                                                                      gCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGc    >  1:7741845/1‑94 (MQ=255)
                                                                                      gCTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGc    <  1:1453240/94‑1 (MQ=255)
                                                                                       cTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCg   >  2:7460679/1‑94 (MQ=255)
                                                                                       cTGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCg   <  1:1209051/94‑1 (MQ=255)
                                                                                        tGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTTTCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGt  <  2:10217623/94‑1 (MQ=255)
                                                                                        tGGGGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGt  >  1:8609045/1‑94 (MQ=255)
                                                                                           |                                                                                          
TTAGCGTTTAGCGCATCGGCGGCGCAGTATGAAGATGGTAAACAGTACACTACCCTGGAAAAACCGGTAGCTGGCGCGCCGCAAGTGCTGGAGTTTTTCTCTTTCTTCTGCCCGCACTGCTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGT  >  NC_000913/4041035‑4041216

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ 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
CAATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCT  >  NC_000913/4041562‑4041717
                                                                                            |                                                                
CAATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAAT                                                                 <  SRR4302545.2253155/94‑1 (MQ=60)
CAATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAAT                                                                 <  SRR4302545.870259/94‑1 (MQ=60)
  ATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATAT                                                               <  SRR4302545.10939550/94‑1 (MQ=60)
  ATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATAT                                                               <  SRR4302545.6407443/94‑1 (MQ=60)
   TATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATT                                                              <  SRR4302545.10601991/94‑1 (MQ=60)
   TATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATT                                                              <  SRR4302545.10849014/94‑1 (MQ=60)
   TATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATT                                                              <  SRR4302545.7923141/94‑1 (MQ=60)
    ATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTC                                                             <  SRR4302545.7157618/94‑1 (MQ=60)
     TGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCA                                                            <  SRR4302545.1533676/94‑1 (MQ=60)
     TGGATGTTTTTGTTAAGCAGCAAGCTGATACAGTGAAATATCTGTCTAAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCA                                                            <  SRR4302545.9070622/94‑1 (MQ=60)
      GGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAC                                                           <  SRR4302545.5518289/94‑1 (MQ=60)
      GGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAC                                                           <  SRR4302545.9727456/94‑1 (MQ=60)
      GGATGTTTTTGTTCAGCAGGATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCAC                                                           <  SRR4302545.10237613/94‑1 (MQ=60)
       GATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACG                                                          <  SRR4302545.10154058/94‑1 (MQ=60)
        ATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGG                                                         <  SRR4302545.8412118/94‑1 (MQ=60)
         TGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGG                                                        <  SRR4302545.6807456/94‑1 (MQ=60)
          GTTTTTGTTGAGCAGTATTCTGATTCAGTTAAATCGCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTTAATATTCACGGGC                                                       <  SRR4302545.3455019/94‑1 (MQ=60)
          GTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC                                                       <  SRR4302545.9027602/94‑1 (MQ=60)
          GTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGC                                                       <  SRR4302545.7515293/94‑1 (MQ=60)
           TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCT                                                     <  SRR4302545.926533/94‑1 (MQ=60)
           TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCT                                                     <  SRR4302545.1142021/94‑1 (MQ=60)
           TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCT                                                     <  SRR4302545.5420022/94‑1 (MQ=60)
           TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCT                                                     <  SRR4302545.4860765/94‑1 (MQ=60)
           TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCT                                                     <  SRR4302545.2952472/94‑1 (MQ=60)
           TTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCT                                                     <  SRR4302545.9423174/94‑1 (MQ=60)
           TTTTAGATCAACAGTGTGCTGATACAGTGAAATATCTGTCAGAAAAAGAAAAATTCAATGTAAATTACATAAAGCACGTGAATATTCACGAGCT                                                     <  SRR4302545.9653784/94‑1 (MQ=42)
            TTTTGTTCAGCATTGTGTTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCCCGGGCTT                                                    <  SRR4302545.7965996/94‑1 (MQ=60)
            TTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTT                                                    <  SRR4302545.155804/94‑1 (MQ=60)
            TTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTT                                                    <  SRR4302545.1452068/94‑1 (MQ=60)
            TTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTT                                                    <  SRR4302545.8973715/94‑1 (MQ=60)
              TTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTT                                                  <  SRR4302545.8877230/94‑1 (MQ=60)
              TTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTT                                                  <  SRR4302545.7173411/94‑1 (MQ=60)
              TTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTT                                                  <  SRR4302545.7338569/94‑1 (MQ=60)
               TGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTT                                                 <  SRR4302545.10656456/94‑1 (MQ=60)
                TTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTT                                                <  SRR4302545.2480110/94‑1 (MQ=60)
                GTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTT                                                <  SRR4302545.132036/94‑1 (MQ=60)
                GTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTT                                                <  SRR4302545.3505300/94‑1 (MQ=60)
                 TTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTT                                               <  SRR4302545.490825/94‑1 (MQ=60)
                 TTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTT                                               <  SRR4302545.4874680/94‑1 (MQ=60)
                 TTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTT                                               <  SRR4302545.1318420/94‑1 (MQ=60)
                  TCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTta                                              <  SRR4302545.4773796/94‑3 (MQ=60)
                   CAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtat                                             <  SRR4302545.5167667/94‑4 (MQ=60)
                   CAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTtat                                             <  SRR4302545.6242919/94‑4 (MQ=60)
                    AGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTT‑‑TTATT                                          <  SRR4302545.6415249/94‑1 (MQ=60)
                     GCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTA                                            <  SRR4302545.3965613/94‑1 (MQ=60)
                     GCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTA                                            <  SRR4302545.4926285/94‑1 (MQ=60)
                      CAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTAT                                           <  SRR4302545.3076604/94‑1 (MQ=60)
                       AGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATT                                          <  SRR4302545.7674792/94‑1 (MQ=60)
                       AGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATT                                          <  SRR4302545.10356139/94‑1 (MQ=60)
                        GTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTT                                         <  SRR4302545.4447859/94‑1 (MQ=60)
                        GTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTT                                         <  SRR4302545.9785766/94‑1 (MQ=60)
                         TATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTA                                        <  SRR4302545.8645829/94‑1 (MQ=60)
                         TATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTA                                        <  SRR4302545.7676280/94‑1 (MQ=60)
                           TGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAAT                                      <  SRR4302545.2156480/94‑1 (MQ=60)
                            GCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATA                                     <  SRR4302545.3581678/94‑1 (MQ=60)
                             CTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAA                                    <  SRR4302545.6489831/94‑1 (MQ=60)
                             CTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAA                                    <  SRR4302545.10803610/94‑1 (MQ=60)
                                       GAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATAC                          >  SRR4302545.7380862/1‑94 (MQ=60)
                                                     AGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCA            <  SRR4302545.9039177/94‑1 (MQ=60)
                                                               ATTCAATGTAAATTACATAAAGCCCGTGAATATTCACGGGCTTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCT  >  SRR4302545.2500086/1‑94 (MQ=60)
                                                                                            |                                                                
CAATATGGATGTTTTTGTTCAGCAGTATGCTGATACAGTGAAATATCTGTCCGAGAAAAAATAATTCAATGTAAATTACATAAAGCCCGTGAGTATTCACGGGCTTTTTTATTATTTAATAAATATAAATACATTCTGATAATGCATCCTGCCGCT  >  NC_000913/4041562‑4041717

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