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 1873 183.1 9376148 98.2% 9207377 93.0

Breseq alignment

BRESEQ :: Evidence
Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,995,916 T→C Y196C (TAC→TGC)  yecS ← putative transport system permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,995,9160TC99.1% 365.4 / ‑4.6 116Y196C (TAC→TGC) yecSputative transport system permease protein
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (73/42);  total (74/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.80e-01

Only 100 of 121 total aligned reads displayed.

AATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGG  >  NC_000913/1995826‑1996005
                                                                                          |                                                                                         
aaTGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCcagcag                                                                                         >  1:513176/1‑93 (MQ=255)
 aTGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGa                                                                                        >  2:3733237/1‑93 (MQ=255)
   ggCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATg                                                                                      <  2:1152257/93‑1 (MQ=255)
     cACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAg                                                                                    >  2:123947/1‑93 (MQ=255)
     cACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAg                                                                                    <  2:3080549/93‑1 (MQ=255)
      aCTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGc                                                                                   >  2:3903263/1‑93 (MQ=255)
      aCTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGc                                                                                   >  1:1125315/1‑93 (MQ=255)
        tCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGa                                                                                 >  2:4380365/1‑93 (MQ=255)
        tCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGa                                                                                 <  2:720838/93‑1 (MQ=255)
         cATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGaa                                                                                <  2:1859735/93‑1 (MQ=255)
         cATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGaa                                                                                <  2:1594449/93‑1 (MQ=255)
          aTTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGaac                                                                               >  1:295124/1‑92 (MQ=255)
           ttttGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgc                                                                              >  1:1450026/1‑93 (MQ=255)
           ttttGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgc                                                                              >  2:4585814/1‑93 (MQ=255)
           ttttGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgc                                                                              <  2:3118863/93‑1 (MQ=255)
           ttttGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgc                                                                              >  2:2939517/1‑93 (MQ=255)
           ttttGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgc                                                                              <  2:1025774/93‑1 (MQ=255)
            tttGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgcc                                                                             <  1:634977/93‑1 (MQ=255)
               ggTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgccgcc                                                                          >  1:405479/1‑93 (MQ=255)
               ggTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAgccgcc                                                                          >  2:3954681/1‑93 (MQ=255)
                 ttCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAg                                                                        <  2:281371/93‑1 (MQ=255)
                 ttCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGACATGATCCAGCAGATGAGCGAAGCCGCCAg                                                                        >  2:4339122/1‑93 (MQ=255)
                  tctctctcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGa                                                                       <  2:2594634/93‑1 (MQ=255)
                  tctctctcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGa                                                                       <  1:456341/93‑1 (MQ=255)
                   ctctctcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGAt                                                                      >  2:2764016/1‑93 (MQ=255)
                    tctctcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATa                                                                     >  1:114310/1‑93 (MQ=255)
                      tctcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACa                                                                   >  1:2500674/1‑93 (MQ=255)
                       ctcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACAt                                                                  >  1:4598568/1‑93 (MQ=255)
                        tcCTGGCTAGTAAGTTGATTCTCAAAATGGTTCTGCAACGGCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATg                                                                 <  1:3853132/93‑1 (MQ=255)
                        tcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATg                                                                 <  1:2033667/93‑1 (MQ=255)
                        tcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATg                                                                 >  1:3689958/1‑93 (MQ=255)
                        tcCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATg                                                                 >  1:3314296/1‑93 (MQ=255)
                          cTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATggc                                                               >  2:1033803/1‑92 (MQ=255)
                          cTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGt                                                               >  1:4484107/1‑93 (MQ=255)
                           tGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTg                                                              >  1:3460257/1‑93 (MQ=255)
                           tGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTg                                                              >  2:1816889/1‑93 (MQ=255)
                           tGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTg                                                              >  2:1746082/1‑93 (MQ=255)
                              cGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAac                                                           <  2:455878/93‑1 (MQ=255)
                               gATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAaca                                                          <  1:4359158/93‑1 (MQ=255)
                               gATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAaca                                                          >  1:1669373/1‑93 (MQ=255)
                               gATTAAGGTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAaca                                                          >  2:3804248/1‑93 (MQ=255)
                                  tAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCt                                                       <  1:1868246/93‑1 (MQ=255)
                                  tAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCt                                                       <  1:894197/93‑1 (MQ=255)
                                   aaGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTc                                                      >  2:1985092/1‑93 (MQ=255)
                                   aaGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTc                                                      >  1:3071509/1‑93 (MQ=255)
                                   aaGTTGATTCTCAAAATGGATCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAAAACCTc                                                      >  2:1160741/1‑93 (MQ=255)
                                    aGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACcccc                                                     >  1:1284438/1‑93 (MQ=255)
                                    aGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTcc                                                     >  2:1572743/1‑93 (MQ=255)
                                     gTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCa                                                    >  1:279761/1‑93 (MQ=255)
                                     gTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCa                                                    >  1:1922609/1‑93 (MQ=255)
                                     gTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCa                                                    >  1:2243278/1‑93 (MQ=255)
                                     gTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCa                                                    >  2:2175660/1‑93 (MQ=255)
                                      ttGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAg                                                   >  2:4062086/1‑93 (MQ=255)
                                         aTTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCga                                                >  1:2708236/1‑92 (MQ=255)
                                         aTTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGt                                                >  2:4182119/1‑93 (MQ=255)
                                            ctcAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTAcg                                             >  2:63032/1‑93 (MQ=255)
                                             tcAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTAcgc                                            >  1:1711033/1‑93 (MQ=255)
                                               aaaaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGa                                          >  2:218771/1‑93 (MQ=255)
                                                aaaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAg                                         >  1:2606542/1‑93 (MQ=255)
                                                aaaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAg                                         <  2:3132684/93‑1 (MQ=255)
                                                aaaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAg                                         >  2:243286/1‑93 (MQ=255)
                                                aaaTGGTTCTGCAACGTCGATAAAACAGTCGCCATGATCCAGCAGATGAGCGACGACGCCAGATACATGGGGAAGACCTCCAGCGTATGCGAg                                         >  2:461693/1‑93 (MQ=255)
                                                aaaTGCTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAg                                         <  1:2623628/93‑1 (MQ=255)
                                                 aaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAgg                                        >  1:4178777/1‑93 (MQ=255)
                                                 aaTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAgg                                        >  2:1057427/1‑93 (MQ=255)
                                                  aTGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGt                                       >  2:2200234/1‑93 (MQ=255)
                                                   tGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTa                                      >  1:4243055/1‑93 (MQ=255)
                                                     gTTCTGCAACGTCTATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAAt                                    <  2:1470935/93‑1 (MQ=255)
                                                      ttCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATc                                   <  2:274800/93‑1 (MQ=255)
                                                        cTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCaa                                 >  2:159379/1‑93 (MQ=255)
                                                         tGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAAc                                <  1:4345690/93‑1 (MQ=255)
                                                          gCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACt                               <  1:2035428/93‑1 (MQ=255)
                                                          gCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACt                               <  1:1536722/93‑1 (MQ=255)
                                                            aaCGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTgc                             <  2:3246989/93‑1 (MQ=255)
                                                             aCGTCGATAACAGAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTgcg                            >  2:3758794/1‑93 (MQ=255)
                                                             aCGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTgcg                            >  1:2331542/1‑93 (MQ=255)
                                                              cGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTgcgc                           <  2:1492639/93‑1 (MQ=255)
                                                               gTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGcc                          >  1:3406345/1‑93 (MQ=255)
                                                                 cGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGCGAACACCTCCAGCGTACGCGAGGCAATCAACTGCGCCTg                        >  2:4144756/1‑93 (MQ=255)
                                                                  gATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGa                       >  2:1948284/1‑93 (MQ=255)
                                                                  gATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGa                       <  2:2939949/93‑1 (MQ=255)
                                                                  gATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGa                       >  2:2747432/1‑93 (MQ=255)
                                                                   aTAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTGCGCGAGGTAATCAACTGCGCCTGAc                      >  1:1860744/1‑93 (MQ=255)
                                                                   aTAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGAc                      >  1:4094790/1‑93 (MQ=255)
                                                                    tAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACg                     >  2:560124/1‑93 (MQ=255)
                                                                      acacAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGa                   >  2:2940786/1‑93 (MQ=255)
                                                                          aGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAg               >  2:2236434/1‑93 (MQ=255)
                                                                              gCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTcc           <  1:1301861/93‑1 (MQ=255)
                                                                              gCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTcc           <  2:3658756/93‑1 (MQ=255)
                                                                              gCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTcc           <  1:1174248/93‑1 (MQ=255)
                                                                               ccATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACCCCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCg          <  1:3673578/93‑1 (MQ=255)
                                                                                cATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCgg         >  2:3463019/1‑93 (MQ=255)
                                                                                 aTGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGc        >  2:2742154/1‑93 (MQ=255)
                                                                                  tGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGcc       >  2:728276/1‑92 (MQ=255)
                                                                                   gATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCAc      <  1:1487397/93‑1 (MQ=255)
                                                                                   gATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCAc      >  2:1717339/1‑93 (MQ=255)
                                                                                    aTCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCAcc     <  2:2343112/93‑1 (MQ=255)
                                                                                       cagcagATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTgg  >  1:963193/1‑93 (MQ=255)
                                                                                       cagcagATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTgg  <  1:605086/93‑1 (MQ=255)
                                                                                       cagcagATGAGCGAAGCCGCCAGATACATGGGGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTgg  >  1:377786/1‑93 (MQ=255)
                                                                                          |                                                                                         
AATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGG  >  NC_000913/1995826‑1996005

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ 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
CAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTG  >  NC_000913/1995825‑1996009
                                                                                           |                                                                                             
CAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAG                                                                                             >  SRR4302535.516762/1‑94 (MQ=60)
CAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAG                                                                                             <  SRR4302535.3357340/94‑1 (MQ=60)
 AATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGA                                                                                            >  SRR4302535.1824011/1‑94 (MQ=60)
      CACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGC                                                                                       >  SRR4302535.1133349/1‑94 (MQ=60)
          CATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAG                                                                                   >  SRR4302535.3727622/1‑94 (MQ=60)
          CATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAC                                                                                   >  SRR4302535.297184/1‑94 (MQ=60)
           ATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGC                                                                                  >  SRR4302535.1460397/1‑94 (MQ=60)
             TTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCG                                                                                <  SRR4302535.639400/94‑1 (MQ=60)
               TGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCC                                                                              >  SRR4302535.408319/1‑94 (MQ=60)
                   TCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGAT                                                                          <  SRR4302535.459533/94‑1 (MQ=60)
                    CTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATA                                                                         >  SRR4302535.115113/1‑94 (MQ=60)
                      CTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACA                                                                       >  SRR4302535.2518707/1‑94 (MQ=60)
                       TCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACAT                                                                      >  SRR4302535.4629850/1‑94 (MQ=60)
                        CTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATG                                                                     >  SRR4302535.3337354/1‑94 (MQ=60)
                        CTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATG                                                                     >  SRR4302535.3715497/1‑94 (MQ=60)
                         TCCTGGCTAGTAAGTTGATTCTCAAAATGGTTCTGCAACGGCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGG                                                                    <  SRR4302535.3879638/94‑1 (MQ=60)
                         TCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGG                                                                    <  SRR4302535.2048320/94‑1 (MQ=60)
                          CCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGT                                                                   >  SRR4302535.4514890/1‑94 (MQ=60)
                           CTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTG                                                                  <  SRR4302535.2181616/94‑1 (MQ=60)
                           CTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTG                                                                  >  SRR4302535.3484314/1‑94 (MQ=60)
                               CGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACA                                                              >  SRR4302535.1681451/1‑94 (MQ=60)
                                GATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACAC                                                             <  SRR4302535.4388997/94‑1 (MQ=60)
                                  ataaagttaatcttaaaaagggtcctgcaccgtgaaaacccaagccgCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCT                                                           <  SRR4302535.2507283/47‑1 (MQ=60)
                                   TAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTC                                                          >  SRR4302535.3092955/1‑94 (MQ=60)
                                   TAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTC                                                          <  SRR4302535.900537/94‑1 (MQ=60)
                                   TAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTC                                                          <  SRR4302535.1881712/94‑1 (MQ=60)
                                    AAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCC                                                         <  SRR4302535.4226466/94‑1 (MQ=60)
                                    AAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCCCC                                                         >  SRR4302535.1293592/1‑94 (MQ=60)
                                     AGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCA                                                        >  SRR4302535.2259232/1‑94 (MQ=60)
                                     AGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCA                                                        >  SRR4302535.281712/1‑94 (MQ=60)
                                     AGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCA                                                        >  SRR4302535.1936469/1‑94 (MQ=60)
                                      GTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAG                                                       >  SRR4302535.812785/1‑94 (MQ=60)
                                         GATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGA                                                    >  SRR4302535.2727588/1‑94 (MQ=60)
                                             CTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGC                                                >  SRR4302535.1723400/1‑94 (MQ=60)
                                                AAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAG                                             >  SRR4302535.2625247/1‑94 (MQ=60)
                                                 AAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGG                                            <  SRR4302535.2642431/94‑1 (MQ=60)
                                                 AAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGG                                            >  SRR4302535.135001/1‑94 (MQ=60)
                                                 AAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGG                                            >  SRR4302535.4207380/1‑94 (MQ=60)
                                                 AAATGCTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGG                                            <  SRR4302535.2642430/94‑1 (MQ=60)
                                                   ATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTA                                          >  SRR4302535.4272095/1‑94 (MQ=60)
                                                       TTCTGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCA                                      >  SRR4302535.4599194/1‑94 (MQ=60)
                                                          TGCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACT                                   <  SRR4302535.4375425/94‑1 (MQ=60)
                                                           GCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTG                                  <  SRR4302535.1547795/94‑1 (MQ=60)
                                                           GCAACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTG                                  <  SRR4302535.2050091/94‑1 (MQ=60)
                                                             AACGTCGATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCG                                >  SRR4302535.2348220/1‑94 (MQ=60)
                                                               CGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCC                              >  SRR4302535.3430053/1‑94 (MQ=60)
                                                                   GATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTGCGCGAGGTAATCAACTGCGCCTGAC                          >  SRR4302535.1874160/1‑94 (MQ=60)
                                                                   GATAACACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGAC                          >  SRR4302535.4122810/1‑94 (MQ=60)
                                                                        CACAGTCGCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAAC                     <  SRR4302535.3102451/94‑1 (MQ=60)
                                                                         acataCGCCATGATCCATCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACA                    <  SRR4302535.2332106/89‑1 (MQ=60)
                                                                               GCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCG              <  SRR4302535.1311160/94‑1 (MQ=60)
                                                                               GCCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCG              <  SRR4302535.1182630/94‑1 (MQ=60)
                                                                                CCATGATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACCCCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGG             <  SRR4302535.3699015/94‑1 (MQ=60)
                                                                                    GATCCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACC         <  SRR4302535.1498050/94‑1 (MQ=60)
                                                                                       CCAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGG      >  SRR4302535.970032/1‑94 (MQ=60)
                                                                                       CCAGCAGATGAGCGAAGCCGCCAGATACATGGGGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGG      >  SRR4302535.380405/1‑94 (MQ=60)
                                                                                        CAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGA     >  SRR4302535.1914217/1‑94 (MQ=60)
                                                                                        CAGCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGA     <  SRR4302535.609318/94‑1 (MQ=60)
                                                                                          GCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATT   <  SRR4302535.4268544/94‑1 (MQ=60)
                                                                                          GCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATT   >  SRR4302535.1053023/1‑94 (MQ=60)
                                                                                          GCAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATT   <  SRR4302535.1991727/94‑1 (MQ=60)
                                                                                           CAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTG  >  SRR4302535.782527/1‑94 (MQ=60)
                                                                                           |                                                                                             
CAATGGCACTCATTTTGGTTCTCTCTCCTGGCGATTAAGTTGATTCTCAAAATGGTTCTGCAACGTCGATAACACAGTCGCCATGATCCAGTAGATGAGCGAAGCCGCCAGATACATGGTGAACACCTCCAGCGTACGCGAGGTAATCAACTGCGCCTGACGGAACAGCTCCGGCACCTGGATTG  >  NC_000913/1995825‑1996009

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