Sample Resequencing Stats

Note: The mutation counts shown below represent unfiltered mutation sets.
ALE, Flask, Isolate Predicted Mutations Mean Coverage Total Reads Percent Mapped Mapped Reads Average Read Length
A1 F1 I2 R1 593 264.2 13311096 97.8% 13018251 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,786,585 G→A G63D (GGC→GAC)  b2658 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,786,5850GA84.2% 276.0 / 36.2 101G63D (GGC→GAC) b2658hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (8/8);  new base A (28/57);  total (36/65)
Fisher's exact test for biased strand distribution p-value = 2.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.42e-01

Only 100 of 103 total aligned reads displayed.

ATACAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAA  >  NC_000913/2786496‑2786675
                                                                                         |                                                                                          
ataCAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACgg                                                                                         <  2:829173/93‑1 (MQ=255)
 taCAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGa                                                                                        >  1:2792497/1‑93 (MQ=255)
  aCAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGaa                                                                                       <  2:2207508/93‑1 (MQ=255)
   cAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAAc                                                                                      <  1:4962873/93‑1 (MQ=255)
     aaaaCATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAActc                                                                                    <  1:3376539/93‑1 (MQ=255)
      aaaCATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAActct                                                                                   >  1:1016758/1‑93 (MQ=255)
         cATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTaaa                                                                                <  2:6305761/93‑1 (MQ=255)
         cATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTaaa                                                                                <  2:3826563/93‑1 (MQ=255)
             gggAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATc                                                                            <  1:5416756/93‑1 (MQ=255)
              ggAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCa                                                                           >  2:4973618/1‑93 (MQ=255)
              ggAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCa                                                                           <  2:2685012/93‑1 (MQ=255)
               gAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCaa                                                                          <  1:3192686/93‑1 (MQ=255)
               gAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCaa                                                                          <  2:3056182/93‑1 (MQ=255)
               gAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCaa                                                                          <  1:6335828/93‑1 (MQ=255)
                 gTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAata                                                                        <  1:129797/93‑1 (MQ=255)
                  tAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAatat                                                                       >  2:1058426/1‑93 (MQ=255)
                   aaTAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATAtt                                                                      <  1:4666991/93‑1 (MQ=255)
                    aTAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTc                                                                     <  2:4693781/93‑1 (MQ=255)
                    aTAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTc                                                                     <  1:6428626/93‑1 (MQ=255)
                    aTAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTc                                                                     >  2:2440304/1‑93 (MQ=255)
                     tAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTcc                                                                    <  1:1275065/93‑1 (MQ=255)
                     tAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTcc                                                                    <  1:4616717/93‑1 (MQ=255)
                      aGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCt                                                                   >  1:5740563/1‑93 (MQ=255)
                       gCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTg                                                                  <  2:2901959/93‑1 (MQ=255)
                       gCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTg                                                                  <  1:3538097/93‑1 (MQ=255)
                        cGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTgg                                                                 <  1:80587/93‑1 (MQ=255)
                        cGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTgg                                                                 <  1:1236867/93‑1 (MQ=255)
                         gTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTggg                                                                <  2:3118884/93‑1 (MQ=255)
                          ttAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGc                                                               >  2:2353728/1‑93 (MQ=255)
                           tAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCg                                                              <  2:1760457/93‑1 (MQ=255)
                           tAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCg                                                              <  2:3970871/93‑1 (MQ=255)
                           tAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCg                                                              <  1:2371923/93‑1 (MQ=255)
                            aaCATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGt                                                             >  1:3700239/1‑93 (MQ=255)
                               attattattAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAgc                                                          >  2:733998/1‑93 (MQ=255)
                                 tattattaATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAgcgc                                                        <  2:5830105/93‑1 (MQ=255)
                                   ttattaATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGctc                                                      <  2:202690/93‑1 (MQ=255)
                                    tattaATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGctct                                                     <  1:3406002/93‑1 (MQ=255)
                                    tattaATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGctct                                                     <  1:4359606/93‑1 (MQ=255)
                                     attaATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCtctt                                                    >  2:5720994/1‑93 (MQ=255)
                                     attaATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCtctt                                                    <  2:1039128/93‑1 (MQ=255)
                                      ttaGTAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCtcttc                                                   <  1:2395596/93‑1 (MQ=255)
                                      ttaATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCtcttc                                                   >  2:1540117/1‑93 (MQ=255)
                                          taataCCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCtttt                                               <  2:1441406/93‑1 (MQ=255)
                                          taataCCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCtttt                                               >  2:5850554/1‑93 (MQ=255)
                                          taataCCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCtttt                                               >  2:798621/1‑93 (MQ=255)
                                            ataCCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAg                                             <  2:6005304/93‑1 (MQ=255)
                                             taCCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAgg                                            <  2:3678109/93‑1 (MQ=255)
                                              aCCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAggg                                           >  1:2902763/1‑93 (MQ=255)
                                               cccGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGa                                          >  1:4463451/1‑93 (MQ=255)
                                                ccGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGaa                                         >  2:6296304/1‑93 (MQ=255)
                                                 cGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGaaa                                        <  1:188550/93‑1 (MQ=255)
                                                  gAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGaaaa                                       <  2:3038768/93‑1 (MQ=255)
                                                  gAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGaaaa                                       <  1:6276913/93‑1 (MQ=255)
                                                   aaaaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAAc                                      <  2:4804722/93‑1 (MQ=255)
                                                     aaaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAg                                    <  2:6537518/93‑1 (MQ=255)
                                                      aaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAgg                                   <  1:5206845/93‑1 (MQ=255)
                                                      aaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAgg                                   >  2:5847361/1‑93 (MQ=255)
                                                      aaTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAgg                                   <  1:864644/93‑1 (MQ=255)
                                                       aTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGga                                  <  1:782145/93‑1 (MQ=255)
                                                       aTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGga                                  >  2:4485574/1‑93 (MQ=255)
                                                       aTTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGga                                  <  1:2678588/93‑1 (MQ=255)
                                                          aTCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGaa                               <  1:1896012/93‑1 (MQ=255)
                                                              ccccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgc                           <  2:3227095/93‑1 (MQ=255)
                                                              ccccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgc                           <  2:5977698/93‑1 (MQ=255)
                                                              ccccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgc                           >  1:1759998/1‑93 (MQ=255)
                                                              ccccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgc                           >  1:3953805/1‑93 (MQ=255)
                                                               cccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgcc                          <  1:1531123/93‑1 (MQ=255)
                                                               cccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgcc                          <  2:3458898/93‑1 (MQ=255)
                                                               cccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgcc                          >  1:3920628/1‑93 (MQ=255)
                                                               cccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgcc                          >  1:1988472/1‑93 (MQ=255)
                                                               cccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgcc                          <  1:3517091/93‑1 (MQ=255)
                                                               cccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgcc                          <  2:2230031/93‑1 (MQ=255)
                                                                ccTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgccg                         >  2:5180430/1‑93 (MQ=255)
                                                                ccTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAccgccg                         >  1:4500694/1‑93 (MQ=255)
                                                                 cTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGa                        >  1:3439094/1‑93 (MQ=255)
                                                                 cTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGa                        >  2:5882283/1‑93 (MQ=255)
                                                                  tGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGaa                       <  2:3169824/93‑1 (MQ=255)
                                                                  tGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGaa                       >  2:4600943/1‑93 (MQ=255)
                                                                     tAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTa                    >  1:3214688/1‑93 (MQ=255)
                                                                     tAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTa                    <  1:2471767/93‑1 (MQ=255)
                                                                     tAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTa                    >  1:2887848/1‑93 (MQ=255)
                                                                      aaaGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAt                   >  1:2905048/1‑93 (MQ=255)
                                                                       aaGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtc                  >  1:5113243/1‑93 (MQ=255)
                                                                       aaGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtc                  <  2:2736376/93‑1 (MQ=255)
                                                                       aaGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtc                  >  1:391570/1‑93 (MQ=255)
                                                                        aGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtct                 <  1:3515637/93‑1 (MQ=255)
                                                                        aGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtct                 <  1:5319505/93‑1 (MQ=255)
                                                                        aGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtct                 <  1:2617580/93‑1 (MQ=255)
                                                                        aGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtct                 <  2:5233154/93‑1 (MQ=255)
                                                                         gTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtctc                >  2:1172534/1‑93 (MQ=255)
                                                                         gTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtctc                <  2:2475424/93‑1 (MQ=255)
                                                                         gTATTTACGCCTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAtctc                <  2:1012551/93‑1 (MQ=255)
                                                                           aTTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTa              <  1:4324509/93‑1 (MQ=255)
                                                                             ttACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTAtt            <  2:5205130/93‑1 (MQ=255)
                                                                                 gCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTtata        <  1:5616393/93‑1 (MQ=255)
                                                                                   ttAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAAt      <  1:6111732/93‑1 (MQ=255)
                                                                                    tAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAAtt     <  1:733998/93‑1 (MQ=255)
                                                                                      gagaCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTa   >  2:466936/1‑93 (MQ=255)
                                                                                      gagaCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTa   >  2:3837096/1‑93 (MQ=255)
                                                                                       agGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTaa  >  1:1835013/1‑93 (MQ=255)
                                                                                         |                                                                                          
ATACAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAA  >  NC_000913/2786496‑2786675

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
GAATATACAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAA  >  NC_000913/2786492‑2786675
                                                                                             |                                                                                          
GAATATACAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGA                                                                                            <  SRR4302547.912742/94‑1 (MQ=60)
 AATATACAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGAC                                                                                           <  SRR4302547.1862984/94‑1 (MQ=60)
    ATACAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGA                                                                                        >  SRR4302547.2848415/1‑94 (MQ=60)
       CAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACT                                                                                     <  SRR4302547.5065573/94‑1 (MQ=60)
         AAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCT                                                                                   >  SRR4302547.1036667/1‑94 (MQ=60)
         AAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCT                                                                                   <  SRR4302547.3445135/94‑1 (MQ=60)
                 GGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCA                                                                           <  SRR4302547.5528813/94‑1 (MQ=60)
                   GAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAAT                                                                         <  SRR4302547.3256303/94‑1 (MQ=60)
                   GAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAAT                                                                         <  SRR4302547.6467749/94‑1 (MQ=60)
                     GTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATAT                                                                       <  SRR4302547.132144/94‑1 (MQ=60)
                       AATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTC                                                                     <  SRR4302547.4763641/94‑1 (MQ=60)
                        ATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCC                                                                    <  SRR4302547.6562641/94‑1 (MQ=60)
                         TAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCT                                                                   <  SRR4302547.1300748/94‑1 (MQ=60)
                         TAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCT                                                                   >  SRR4302547.5859192/1‑94 (MQ=60)
                         TAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCT                                                                   <  SRR4302547.4712313/94‑1 (MQ=60)
                           GCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGG                                                                 <  SRR4302547.3610324/94‑1 (MQ=60)
                            CGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGG                                                                <  SRR4302547.1261833/94‑1 (MQ=60)
                            CGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGG                                                                <  SRR4302547.82091/94‑1 (MQ=60)
                               TAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGT                                                             <  SRR4302547.2419997/94‑1 (MQ=60)
                               TAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGT                                                             >  SRR4302547.3775752/1‑94 (MQ=60)
                                        TATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTT                                                    <  SRR4302547.3475302/94‑1 (MQ=60)
                                        TATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTT                                                    <  SRR4302547.4449149/94‑1 (MQ=60)
                                          TTAGTAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCT                                                  <  SRR4302547.2444121/94‑1 (MQ=60)
                                                 TACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGG                                           >  SRR4302547.2960645/1‑94 (MQ=60)
                                                  ACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGA                                          >  SRR4302547.4555688/1‑94 (MQ=60)
                                                     CGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAA                                       <  SRR4302547.192104/94‑1 (MQ=60)
                                                      GAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAAC                                      <  SRR4302547.6407319/94‑1 (MQ=60)
                                                          AATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGA                                  <  SRR4302547.5313989/94‑1 (MQ=60)
                                                          AATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGA                                  <  SRR4302547.881188/94‑1 (MQ=60)
                                                           ATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAG                                 <  SRR4302547.796933/94‑1 (MQ=60)
                                                           ATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAG                                 <  SRR4302547.2732378/94‑1 (MQ=60)
                                                              ATCCCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAAC                              <  SRR4302547.1933506/94‑1 (MQ=60)
                                                                 CCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGC                           >  SRR4302547.1794875/1‑94 (MQ=60)
                                                                 CCCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGC                           >  SRR4302547.4034220/1‑94 (MQ=60)
                                                                  CCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCC                          >  SRR4302547.4000349/1‑94 (MQ=60)
                                                                  CCCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCC                          >  SRR4302547.2027954/1‑94 (MQ=60)
                                                                   CCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCG                         <  SRR4302547.1561597/94‑1 (MQ=60)
                                                                   CCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCG                         <  SRR4302547.3588885/94‑1 (MQ=60)
                                                                   CCCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCG                         >  SRR4302547.4593656/1‑94 (MQ=60)
                                                                    CCTGGTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGA                        >  SRR4302547.3509218/1‑94 (MQ=60)
                                                                        GTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTA                    >  SRR4302547.2945464/1‑94 (MQ=60)
                                                                        GTAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTA                    >  SRR4302547.3278819/1‑94 (MQ=60)
                                                                         TAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAT                   <  SRR4302547.2521730/94‑1 (MQ=60)
                                                                         TAAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTAT                   >  SRR4302547.2962978/1‑94 (MQ=60)
                                                                          AAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATC                  >  SRR4302547.5218615/1‑94 (MQ=60)
                                                                          AAAGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATC                  >  SRR4302547.398842/1‑94 (MQ=60)
                                                                            AGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTC                <  SRR4302547.5429187/94‑1 (MQ=60)
                                                                            AGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTC                >  SRR4302547.4560041/1‑94 (MQ=60)
                                                                            AGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTC                <  SRR4302547.3587407/94‑1 (MQ=60)
                                                                            AGTATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTC                <  SRR4302547.2670318/94‑1 (MQ=60)
                                                                               ATTTACGCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTAT             <  SRR4302547.4413530/94‑1 (MQ=60)
                                                                                     GCTTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAA       <  SRR4302547.5732294/94‑1 (MQ=60)
                                                                                       TTAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATT     <  SRR4302547.6238512/94‑1 (MQ=60)
                                                                                        TAGAGACGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTT    <  SRR4302547.747711/94‑1 (MQ=60)
                                                                                          GAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAA  >  SRR4302547.1871290/1‑94 (MQ=60)
                                                                                             |                                                                                          
GAATATACAAAAACATGGGGAGTAATAGCGTTAACATTATTATTAATAATACCCGAAAAATTATCCCCCCTGGTAAAGTATTTACGCTTAGAGGCGGAACTCTAAATATCAATATTCCTGGGCGTAGCGCTCTTCTTTTAGGGAAAACAGGAGAACCGCCGAACTATCTCTATTTATAATTTAA  >  NC_000913/2786492‑2786675

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