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,890,225 C→T R345R (CGG→CGA b1808 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,890,2250CT100.0% 403.9 / NA 112R345R (CGG→CGAb1808hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (58/54);  total (58/54)

Only 100 of 116 total aligned reads displayed.

AGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTT  >  NC_000913/1890136‑1890313
                                                                                         |                                                                                        
aGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgt                                                                                       >  2:2287537/1‑93 (MQ=255)
 gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                      >  1:3171274/1‑93 (MQ=255)
 gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                      <  1:3632242/93‑1 (MQ=255)
 gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                      >  1:2209084/1‑93 (MQ=255)
 gcggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtg                                                                                      >  2:1894088/1‑93 (MQ=255)
  cggcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCgtgt                                                                                     <  1:1664379/93‑1 (MQ=255)
    gcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGc                                                                                   >  1:4670628/1‑93 (MQ=255)
    gcgTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGc                                                                                   >  1:4580538/1‑93 (MQ=255)
      gTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGa                                                                                 <  1:1738752/93‑1 (MQ=255)
      gTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGa                                                                                 <  1:4316217/93‑1 (MQ=255)
       tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                >  2:1329190/1‑93 (MQ=255)
       tGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAt                                                                                <  1:2311019/93‑1 (MQ=255)
        gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                               <  1:2112998/93‑1 (MQ=255)
        gagaGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATa                                                                               <  2:3804156/93‑1 (MQ=255)
            gCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAAc                                                                           >  2:663275/1‑93 (MQ=255)
             ccggagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACg                                                                          >  1:1179444/4‑93 (MQ=255)
              cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGt                                                                         >  2:4339933/1‑93 (MQ=255)
              cAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGt                                                                         >  2:653998/1‑93 (MQ=255)
               agagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTg                                                                        <  1:3204121/93‑1 (MQ=255)
                gagTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTgg                                                                       >  1:1468494/1‑93 (MQ=255)
                 agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                      <  1:550963/93‑1 (MQ=255)
                 agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                      <  2:1162368/93‑1 (MQ=255)
                 agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                      <  2:2449129/93‑1 (MQ=255)
                 agTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgc                                                                      >  1:2141808/1‑93 (MQ=255)
                  gTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTTTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcg                                                                     >  1:144208/1‑93 (MQ=255)
                   tAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgc                                                                    >  2:3634963/1‑93 (MQ=255)
                   tAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgc                                                                    <  2:3789690/93‑1 (MQ=255)
                   tAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgc                                                                    >  1:2467595/1‑93 (MQ=255)
                    aaaaTGGCCCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgcg                                                                   >  2:4584868/1‑93 (MQ=255)
                     aaaTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgcgc                                                                  <  2:2177831/93‑1 (MQ=255)
                      aaTGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGgcgcgtt                                                                 <  2:2960098/93‑3 (MQ=255)
                          gcgcGGAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                             <  1:4487660/93‑1 (MQ=255)
                          gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                             <  1:2331765/93‑1 (MQ=255)
                          gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                             <  1:2255625/93‑1 (MQ=255)
                          gcgcGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCaaa                                                             <  1:2544587/93‑1 (MQ=255)
                           cgcgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAt                                                            >  1:4123790/1‑93 (MQ=255)
                            gcgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATg                                                           >  2:3159994/1‑93 (MQ=255)
                            gcgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATg                                                           <  2:1773827/93‑1 (MQ=255)
                             cgAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGc                                                          <  2:3961448/93‑1 (MQ=255)
                              gAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGcc                                                         >  2:2140063/1‑93 (MQ=255)
                                aGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGc                                                       >  2:4580200/1‑93 (MQ=255)
                                 gTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCa                                                      >  1:3210564/1‑93 (MQ=255)
                                    cATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATcc                                                   >  2:211794/1‑93 (MQ=255)
                                     aTTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCa                                                  >  1:1062476/1‑93 (MQ=255)
                                       tCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGc                                                <  1:1045582/93‑1 (MQ=255)
                                       tCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGc                                                >  2:3076639/1‑93 (MQ=255)
                                        cGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCa                                               <  1:3768920/93‑1 (MQ=255)
                                         gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                              <  1:962926/93‑1 (MQ=255)
                                         gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                              <  1:2956988/93‑1 (MQ=255)
                                         gTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCaa                                              <  2:2927294/93‑1 (MQ=255)
                                          tACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAg                                             >  2:44699/1‑93 (MQ=255)
                                            ccAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGc                                           <  2:1537019/93‑1 (MQ=255)
                                                tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                       <  2:4545732/93‑1 (MQ=255)
                                                tAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGaa                                       <  1:2522564/93‑1 (MQ=255)
                                                  gCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACg                                     >  2:74543/1‑93 (MQ=255)
                                                   cTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGc                                    >  1:119937/1‑93 (MQ=255)
                                                       aGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCa                                >  2:3536249/1‑93 (MQ=255)
                                                        gCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAg                               >  2:3056150/1‑93 (MQ=255)
                                                         cccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGt                              <  2:4275601/93‑1 (MQ=255)
                                                          ccGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTg                             <  1:2230479/93‑1 (MQ=255)
                                                           cGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGa                            >  1:3092900/1‑93 (MQ=255)
                                                             gCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCCGTGACa                          >  2:1624270/1‑93 (MQ=255)
                                                              cTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAg                         <  2:2043996/93‑1 (MQ=255)
                                                                gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGGCGCATCCAGCAAAGCGGAACGGCCCAGTGagaggt                       >  2:2509781/1‑91 (MQ=255)
                                                                gattgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGtt                       >  1:1480191/1‑93 (MQ=255)
                                                                 attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                      >  2:3855902/1‑93 (MQ=255)
                                                                 attgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGttt                      <  2:4596635/93‑1 (MQ=255)
                                                                  ttgatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTc                     <  2:1664133/93‑1 (MQ=255)
                                                                    gatCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGc                   >  2:2560253/1‑93 (MQ=255)
                                                                     atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                  <  1:943832/93‑1 (MQ=255)
                                                                     atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                  >  2:2248668/1‑93 (MQ=255)
                                                                     atCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGcc                  >  2:3569264/1‑93 (MQ=255)
                                                                      tctcTTTTAGCCGCTTCAGTCGTGTGCGATACCACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCa                 >  2:2133020/1‑93 (MQ=255)
                                                                      tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGTCCCAGTGACAGTTTTGctc                 >  2:786685/1‑91 (MQ=255)
                                                                      tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCa                 >  2:1542163/1‑93 (MQ=255)
                                                                      tctcTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTCACAGCTTCGCCa                 >  2:4526735/1‑93 (MQ=255)
                                                                        tctTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACg               >  1:4421539/1‑93 (MQ=255)
                                                                          ttttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTc             >  2:3107301/1‑93 (MQ=255)
                                                                          ttttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTc             <  1:1826284/93‑1 (MQ=255)
                                                                           tttAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCa            >  2:4342975/1‑93 (MQ=255)
                                                                             tAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATa          <  1:1847472/93‑1 (MQ=255)
                                                                             tAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATa          <  2:1781405/93‑1 (MQ=255)
                                                                              ggCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCGGTGACAGTTTCGCCACGTCATaa         <  2:2969418/92‑1 (MQ=255)
                                                                              aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa         <  2:994981/93‑1 (MQ=255)
                                                                              aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa         <  2:1992628/93‑1 (MQ=255)
                                                                              aGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATaa         <  1:1855771/93‑1 (MQ=255)
                                                                               gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc        <  1:3475613/93‑1 (MQ=255)
                                                                               gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc        >  2:3438412/1‑93 (MQ=255)
                                                                               gCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAc        <  1:296529/93‑1 (MQ=255)
                                                                                 cGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACaa      <  1:3680507/93‑1 (MQ=255)
                                                                                  gCTTCAGTCGTGTGCGATACACCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACaaa     <  1:4252140/93‑1 (MQ=255)
                                                                                   cTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAg    >  2:891549/1‑93 (MQ=255)
                                                                                   cTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAg    <  2:2451388/93‑1 (MQ=255)
                                                                                    ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt   >  2:3992820/1‑93 (MQ=255)
                                                                                    ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt   >  2:1825807/1‑93 (MQ=255)
                                                                                    ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt   >  1:70530/1‑93 (MQ=255)
                                                                                    ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt   >  2:1024868/1‑93 (MQ=255)
                                                                                    ttCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt   >  2:256312/1‑93 (MQ=255)
                                                                                    tgcAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGt   <  2:2693921/91‑1 (MQ=255)
                                                                                     tCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGtt  <  2:1821244/93‑1 (MQ=255)
                                                                                         |                                                                                        
AGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTT  >  NC_000913/1890136‑1890313

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
CTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGG  >  NC_000913/1890133‑1890318
                                                                                            |                                                                                             
CTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTC                                                                                              <  SRR4302535.3025960/94‑1 (MQ=60)
CTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTC                                                                                              >  SRR4302535.559635/1‑94 (MQ=60)
CTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTC                                                                                              >  SRR4302535.4335415/1‑94 (MQ=60)
 TGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCG                                                                                             <  SRR4302535.1049421/94‑1 (MQ=60)
   AGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTG                                                                                           >  SRR4302535.3193232/1‑94 (MQ=60)
   AGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTG                                                                                           >  SRR4302535.2224819/1‑94 (MQ=60)
    GCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGT                                                                                          <  SRR4302535.3657416/94‑1 (MQ=60)
     CGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTG                                                                                         <  SRR4302535.1676407/94‑1 (MQ=60)
      GGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGC                                                                                        >  SRR4302535.4611735/1‑94 (MQ=60)
      GGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGC                                                                                        >  SRR4302535.4702175/1‑94 (MQ=60)
         GTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAT                                                                                     <  SRR4302535.1751318/94‑1 (MQ=60)
         GTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGAT                                                                                     <  SRR4302535.4345743/94‑1 (MQ=60)
          TGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATA                                                                                    <  SRR4302535.2327557/94‑1 (MQ=60)
           GAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATAC                                                                                   <  SRR4302535.2128248/94‑1 (MQ=60)
               GCCGGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACG                                                                               >  SRR4302535.1187860/1‑94 (MQ=60)
                  AGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGG                                                                            <  SRR4302535.3226311/94‑1 (MQ=60)
                  AGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGG                                                                            >  SRR4302535.1479015/1‑94 (MQ=60)
                   GAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGC                                                                           >  SRR4302535.2157221/1‑94 (MQ=60)
                   GAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGC                                                                           >  SRR4302535.1616609/1‑94 (MQ=60)
                    AGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTTTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCG                                                                          >  SRR4302535.145213/1‑94 (MQ=60)
                    AGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCG                                                                          <  SRR4302535.554819/94‑1 (MQ=60)
                     GTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGC                                                                         >  SRR4302535.2485384/1‑94 (MQ=60)
                             GCGCGGAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 <  SRR4302535.4518468/94‑1 (MQ=60)
                             GCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 <  SRR4302535.2271684/94‑1 (MQ=60)
                             GCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 <  SRR4302535.2348443/94‑1 (MQ=60)
                             GCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 >  SRR4302535.4152009/1‑94 (MQ=60)
                             GCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAAT                                                                 <  SRR4302535.2562907/94‑1 (MQ=60)
                                   AGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCA                                                           >  SRR4302535.3232802/1‑94 (MQ=60)
                                       CATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCA                                                       >  SRR4302535.1070038/1‑94 (MQ=60)
                                          TCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCA                                                    <  SRR4302535.1052990/94‑1 (MQ=60)
                                           CGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAA                                                   <  SRR4302535.3794934/94‑1 (MQ=60)
                                            GTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAG                                                  <  SRR4302535.2977756/94‑1 (MQ=60)
                                            GTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAG                                                  <  SRR4302535.969765/94‑1 (MQ=60)
                                                   TAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAAC                                           <  SRR4302535.2540741/94‑1 (MQ=60)
                                                     GCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGC                                         >  SRR4302535.120777/1‑94 (MQ=60)
                                                            CCCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTG                                  <  SRR4302535.3181920/94‑1 (MQ=60)
                                                             CCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGA                                 >  SRR4302535.3114439/1‑94 (MQ=60)
                                                             CCGGCTGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGA                                 <  SRR4302535.2246338/94‑1 (MQ=60)
                                                                  TGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTT                            >  SRR4302535.1490780/1‑94 (MQ=60)
                                                                  TGATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTT                            >  SRR4302535.2437857/1‑94 (MQ=60)
                                                                   GATTGATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTT                           <  SRR4302535.4653413/94‑1 (MQ=60)
                                                                        ATCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCA                      <  SRR4302535.950510/94‑1 (MQ=60)
                                                                         TCTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCAC                     <  SRR4302535.2130316/94‑1 (MQ=60)
                                                                          CTCTTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACG                    >  SRR4302535.4451818/1‑94 (MQ=60)
                                                                             TTTTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCA                 <  SRR4302535.1839466/94‑1 (MQ=60)
                                                                               TTAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATA               <  SRR4302535.4147748/94‑1 (MQ=60)
                                                                                TAGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAA              <  SRR4302535.1860785/94‑1 (MQ=60)
                                                                                 AGCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAAC             <  SRR4302535.1869139/94‑1 (MQ=60)
                                                                                  GCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACA            <  SRR4302535.298594/94‑1 (MQ=60)
                                                                                  GCCGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACA            <  SRR4302535.3499766/94‑1 (MQ=60)
                                                                                    CGCTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAA          <  SRR4302535.3705979/94‑1 (MQ=60)
                                                                                     GCTTCAGTCGTGTGCGATACACCGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAG         <  SRR4302535.4281277/94‑1 (MQ=60)
                                                                                      CTTCAGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGT        >  SRR4302535.71046/1‑94 (MQ=60)
                                                                                          AGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCA    >  SRR4302535.4553668/1‑94 (MQ=60)
                                                                                          AGTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCA    >  SRR4302535.1263937/1‑94 (MQ=60)
                                                                                           GTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAG   >  SRR4302535.1481809/1‑94 (MQ=60)
                                                                                           GTCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAG   >  SRR4302535.999457/1‑94 (MQ=60)
                                                                                            TCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGG  <  SRR4302535.3845371/94‑1 (MQ=60)
                                                                                            |                                                                                             
CTGAGCGGCGTGAGAGCCAGAGTAAAATGGCGCGAAGTGCATTCGTACCAGTAGCTGTAGCCCGGCTGATTGATCTCTTTTAGCCGCTTCAGCCGTGTGCGATACAACGTGGCGCGCTCAAATGCCGCATCCAGCAAGGCGGAACGCCCCAGTGACAGTTTCGCCACGTCATAACAAAGTTCCAGG  >  NC_000913/1890133‑1890318

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