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,863,065 G→A A199V (GCC→GTC)  b1781 ← putative aldehyde reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,863,0650GA99.2% 415.8 / ‑5.5 121A199V (GCC→GTC) b1781putative aldehyde reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (55/64);  total (57/64)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

Only 100 of 129 total aligned reads displayed.

GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGA  >  NC_000913/1862975‑1863154
                                                                                          |                                                                                         
gCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGAcc                                                                                         <  1:3634959/93‑1 (MQ=255)
 ccAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCt                                                                                        >  1:3871157/1‑93 (MQ=255)
  caacaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTg                                                                                       >  1:602350/1‑93 (MQ=255)
    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTggg                                                                                     >  1:4171495/1‑93 (MQ=255)
    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTggg                                                                                     >  1:1093339/1‑93 (MQ=255)
    acaaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTggg                                                                                     >  1:2732965/1‑93 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  2:946280/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  1:3123741/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  1:2857776/93‑1 (MQ=255)
     caaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGc                                                                                    <  2:1362615/93‑1 (MQ=255)
      aaTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCt                                                                                   >  2:4399561/1‑93 (MQ=255)
       aTACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTa                                                                                  >  2:3508934/1‑93 (MQ=255)
        tACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTaa                                                                                 <  2:3710776/93‑1 (MQ=255)
         aCTTGTGCCGCGCTTATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAAc                                                                                <  2:912934/93‑1 (MQ=255)
         aCTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAAc                                                                                <  1:2322225/93‑1 (MQ=255)
           ttGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACgg                                                                              >  2:742536/1‑93 (MQ=255)
            tgtgCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGa                                                                             >  1:3235632/1‑93 (MQ=255)
               gCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTg                                                                          >  2:4119276/1‑93 (MQ=255)
                ccGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGt                                                                         >  2:3643949/1‑93 (MQ=255)
                ccGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGt                                                                         <  1:3449333/93‑1 (MQ=255)
                ccGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGt                                                                         >  2:4594798/1‑93 (MQ=255)
                   cgcTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAg                                                                      <  2:4427260/93‑1 (MQ=255)
                     cTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGcc                                                                    <  2:4185893/93‑1 (MQ=255)
                      tGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCa                                                                   >  1:53272/1‑93 (MQ=255)
                       gATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCAt                                                                  >  2:378606/1‑93 (MQ=255)
                       gATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCAt                                                                  <  1:126868/93‑1 (MQ=255)
                        atatTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATc                                                                 >  2:1205168/1‑93 (MQ=255)
                         tatTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCa                                                                <  1:1427916/93‑1 (MQ=255)
                         tatTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCa                                                                >  1:4411342/1‑93 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               >  1:71497/1‑93 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               >  1:510317/1‑93 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               <  1:1841395/93‑1 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               <  2:150229/93‑1 (MQ=255)
                          atTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAc                                                               <  1:1016322/93‑1 (MQ=255)
                           ttGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAcc                                                              >  2:3330646/1‑93 (MQ=255)
                           ttGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAcc                                                              >  1:1054583/1‑93 (MQ=255)
                            tgtgAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCg                                                             <  2:1098302/93‑1 (MQ=255)
                             gtgAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCgg                                                            <  1:1472975/93‑1 (MQ=255)
                                aGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCAt                                                         <  1:2043194/93‑1 (MQ=255)
                                   aTGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctg                                                      >  2:3703363/1‑93 (MQ=255)
                                    tgtgCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctgc                                                     <  2:3932248/93‑1 (MQ=255)
                                      tgCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATctgctg                                                   >  2:3976940/1‑93 (MQ=255)
                                       gCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGt                                                  <  2:3356191/93‑1 (MQ=255)
                                           tttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCt                                              <  1:1908594/93‑1 (MQ=255)
                                           tttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCt                                              <  2:484252/93‑1 (MQ=255)
                                            ttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTg                                             <  2:3219798/93‑1 (MQ=255)
                                            ttCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTg                                             >  2:1984661/1‑93 (MQ=255)
                                             tCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTgg                                            <  1:2944119/93‑1 (MQ=255)
                                              cGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGc                                           >  2:1592263/1‑93 (MQ=255)
                                               gTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCa                                          >  2:470960/1‑93 (MQ=255)
                                               gTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCa                                          >  2:2317385/1‑93 (MQ=255)
                                                  gACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCa                                       <  1:920079/93‑1 (MQ=255)
                                                   aCTACCGCGTTTTTTAACAGTCCATTGCGCACCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAg                                      <  2:610729/93‑1 (MQ=255)
                                                   aCTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAg                                      >  2:1670630/1‑93 (MQ=255)
                                                   aCTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAg                                      <  1:3929818/93‑1 (MQ=255)
                                                     tACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAggg                                    <  1:3805571/93‑1 (MQ=255)
                                                     tACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAggg                                    <  1:3206689/93‑1 (MQ=255)
                                                       ccGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGa                                  <  2:2156036/93‑1 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 >  1:1739060/1‑93 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 <  2:3669081/93‑1 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 >  2:2180478/1‑93 (MQ=255)
                                                        cgcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAg                                 <  1:1140662/93‑1 (MQ=255)
                                                         gcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGt                                >  1:2004028/1‑93 (MQ=255)
                                                         gcgTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGt                                <  1:3074879/93‑1 (MQ=255)
                                                           gTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGCGTAg                              >  2:204252/1‑93 (MQ=255)
                                                            gtttttAACAGTCCATTGCGCAACCGCCCGACGTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGa                             <  1:1506034/92‑1 (MQ=255)
                                                             tttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGAt                            <  1:3160584/93‑1 (MQ=255)
                                                             tttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGAt                            <  2:71679/93‑1 (MQ=255)
                                                             tttttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGCGAGGAGAt                            >  1:1121468/1‑93 (MQ=255)
                                                               tttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGGGTAGATCg                          <  2:2496426/93‑1 (MQ=255)
                                                                ttAACAGTCCATTGCGCAACCGCCCGTCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCGGGCACCAGGGGAGTAGATCGt                         <  2:2051019/93‑1 (MQ=255)
                                                                ttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATATGCTGTTGCTGGCACCAGGGGAGTAGATCGt                         >  1:202179/1‑93 (MQ=255)
                                                                ttAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATAACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGCTCGt                         >  2:2022047/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        >  1:38271/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        >  1:1134614/1‑93 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        <  1:2550426/93‑1 (MQ=255)
                                                                 tAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTa                        <  2:3158799/93‑1 (MQ=255)
                                                                   aCAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACt                      >  1:945130/1‑93 (MQ=255)
                                                                   aCAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACt                      >  2:1436309/1‑93 (MQ=255)
                                                                       tCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAAt                  <  1:678369/93‑1 (MQ=255)
                                                                         cATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTc                >  2:1715843/1‑93 (MQ=255)
                                                                         cATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTc                >  2:4306542/1‑93 (MQ=255)
                                                                          aTTGCCCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTcc               <  1:2461304/93‑1 (MQ=255)
                                                                           ttGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTCCTCAATTCCt              >  1:324707/1‑93 (MQ=255)
                                                                           ttGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCt              <  2:2161699/93‑1 (MQ=255)
                                                                            tGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTc             <  1:3593763/93‑1 (MQ=255)
                                                                             gcgcAACCGCCGGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCg            <  2:4051661/93‑1 (MQ=255)
                                                                             gcgcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCg            <  1:2829444/93‑1 (MQ=255)
                                                                               gcAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTg          >  2:4033830/1‑93 (MQ=255)
                                                                                cAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGa         >  2:4613743/1‑93 (MQ=255)
                                                                                cAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGa         >  1:3956039/1‑93 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       >  1:1500752/1‑93 (MQ=255)
                                                                                  aCCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAc       <  2:2313168/93‑1 (MQ=255)
                                                                                   ccGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAcc      <  1:2794315/93‑1 (MQ=255)
                                                                                     gCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCga    <  2:302419/93‑1 (MQ=255)
                                                                                     gCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCga    <  1:4596262/93‑1 (MQ=255)
                                                                                      cccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCACTTCCTCCTGAACCgag   >  2:714079/1‑93 (MQ=255)
                                                                                      cccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgag   <  2:1805083/93‑1 (MQ=255)
                                                                                       ccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgaga  >  1:3114398/1‑93 (MQ=255)
                                                                                       ccGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCgaga  <  2:3154295/93‑1 (MQ=255)
                                                                                          |                                                                                         
GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGA  >  NC_000913/1862975‑1863154

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

GATK/CNVnator alignment

BRESEQ :: bam2aln output
CGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGATGG  >  NC_000913/1862974‑1863157
                                                                                           |                                                                                            
CGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACC                                                                                            <  SRR4302535.2286373/94‑1 (MQ=60)
 GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCT                                                                                           <  SRR4302535.3660148/94‑1 (MQ=60)
 GCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCT                                                                                           >  SRR4302535.3897759/1‑94 (MQ=60)
  CCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTG                                                                                          >  SRR4302535.606566/1‑94 (MQ=60)
    AACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGG                                                                                        >  SRR4302535.4200053/1‑94 (MQ=60)
    AACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGG                                                                                        >  SRR4302535.2752471/1‑94 (MQ=60)
    AACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGG                                                                                        >  SRR4302535.1101125/1‑94 (MQ=60)
      CAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCT                                                                                      <  SRR4302535.2878022/94‑1 (MQ=60)
      CAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCT                                                                                      <  SRR4302535.3145496/94‑1 (MQ=60)
          ACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACG                                                                                  <  SRR4302535.1396325/94‑1 (MQ=60)
          ACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACG                                                                                  <  SRR4302535.2338842/94‑1 (MQ=60)
          ACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACG                                                                                  <  SRR4302535.2776093/94‑1 (MQ=60)
            TTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGA                                                                                >  SRR4302535.3258093/1‑94 (MQ=60)
                 CCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTA                                                                           <  SRR4302535.3473314/94‑1 (MQ=60)
                      CTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCA                                                                      >  SRR4302535.53663/1‑94 (MQ=60)
                        GATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATC                                                                    <  SRR4302535.127762/94‑1 (MQ=60)
                         ATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCA                                                                   >  SRR4302535.4441533/1‑94 (MQ=60)
                          TATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGCCCTGGGCTAACGGACTGTAAGCCATCAC                                                                  >  SRR4302535.3162299/1‑94 (MQ=60)
                          TATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAC                                                                  >  SRR4302535.72019/1‑94 (MQ=60)
                          TATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAC                                                                  >  SRR4302535.513886/1‑94 (MQ=60)
                          TATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCAC                                                                  <  SRR4302535.1438107/94‑1 (MQ=60)
                           ATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACC                                                                 <  SRR4302535.1023545/94‑1 (MQ=60)
                           ATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACC                                                                 <  SRR4302535.1854675/94‑1 (MQ=60)
                           ATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACC                                                                 >  SRR4302535.1062083/1‑94 (MQ=60)
                              GTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGC                                                              <  SRR4302535.1483520/94‑1 (MQ=60)
                                 AGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATC                                                           <  SRR4302535.2057913/94‑1 (MQ=60)
                                     TGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCTCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCT                                                       <  SRR4302535.1135832/94‑1 (MQ=60)
                                           ATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCT                                                 <  SRR4302535.2355575/94‑1 (MQ=60)
                                            TTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTG                                                <  SRR4302535.1922353/94‑1 (MQ=60)
                                              TCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGC                                              <  SRR4302535.2964815/94‑1 (MQ=60)
                                                   GACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAG                                         <  SRR4302535.926604/94‑1 (MQ=60)
                                                    ACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGG                                        <  SRR4302535.3956657/94‑1 (MQ=60)
                                                    ACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGG                                        <  SRR4302535.159147/94‑1 (MQ=60)
                                                      TACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGG                                      <  SRR4302535.3831802/94‑1 (MQ=60)
                                                      TACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGG                                      <  SRR4302535.3228899/94‑1 (MQ=60)
                                                        CCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAG                                    >  SRR4302535.1751627/1‑94 (MQ=60)
                                                        CCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAG                                    >  SRR4302535.304453/1‑94 (MQ=60)
                                                         CGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGT                                   <  SRR4302535.1148789/94‑1 (MQ=60)
                                                         CGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGT                                   >  SRR4302535.2018470/1‑94 (MQ=60)
                                                          GCGTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTA                                  <  SRR4302535.3096338/94‑1 (MQ=60)
                                                            GTTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGA                                >  SRR4302535.4574771/1‑94 (MQ=60)
                                                             TTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGAT                               <  SRR4302535.3394189/94‑1 (MQ=60)
                                                             TTTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGCGAGGAGAT                               >  SRR4302535.1129467/1‑94 (MQ=60)
                                                             GTTTTTAACAGTCCATTGCGCAACCGCCCGACGTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGAT                               <  SRR4302535.1516861/94‑1 (MQ=60)
                                                              TTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATC                              <  SRR4302535.3182511/94‑1 (MQ=60)
                                                              TTTTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATC                              <  SRR4302535.2004274/94‑1 (MQ=60)
                                                                TTTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATATGCTGTTGCTGGCACCAGGGGAGTAGATCGT                            >  SRR4302535.203560/1‑94 (MQ=60)
                                                                 TTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTA                           <  SRR4302535.566027/94‑1 (MQ=60)
                                                                 TTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTA                           >  SRR4302535.1142710/1‑94 (MQ=60)
                                                                 TTAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTA                           >  SRR4302535.38567/1‑94 (MQ=60)
                                                                  TAACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTAC                          <  SRR4302535.2568788/94‑1 (MQ=60)
                                                                   AACAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACT                         >  SRR4302535.951823/1‑94 (MQ=60)
                                                                     CAGTCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCA                       <  SRR4302535.3913383/94‑1 (MQ=60)
                                                                        TCCATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATT                    <  SRR4302535.683113/94‑1 (MQ=60)
                                                                           ATTGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTCCTCAATTCCT                 >  SRR4302535.326954/1‑94 (MQ=60)
                                                                           ATTGCCCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCT                 <  SRR4302535.2479048/94‑1 (MQ=60)
                                                                             TGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCG               <  SRR4302535.3684054/94‑1 (MQ=60)
                                                                             TGCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCG               <  SRR4302535.3618704/94‑1 (MQ=60)
                                                                              GCGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGT              <  SRR4302535.2849528/94‑1 (MQ=60)
                                                                               CGCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTG             >  SRR4302535.939792/1‑94 (MQ=60)
                                                                                GCAACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGA            >  SRR4302535.3983053/1‑94 (MQ=60)
                                                                                  AACCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAAC          >  SRR4302535.1511548/1‑94 (MQ=60)
                                                                                    CCGCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCG        <  SRR4302535.2814195/94‑1 (MQ=60)
                                                                                      GCCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAG      <  SRR4302535.4627535/94‑1 (MQ=60)
                                                                                       CCCGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGA     >  SRR4302535.3136080/1‑94 (MQ=60)
                                                                                         CGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGATG   <  SRR4302535.1929441/94‑1 (MQ=60)
                                                                                         CGACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGATG   <  SRR4302535.2348830/94‑1 (MQ=60)
                                                                                          GACCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGATGG  >  SRR4302535.2456343/1‑94 (MQ=60)
                                                                                           |                                                                                            
CGCCAACAATACTTGTGCCGCGCTGATATTGTGAGCATGTGCAATTTCGTTGACTACCGCGTTTTTTAACAGTCCATTGCGCAACCGCCCGGCCTGGGCTAACGGACTGTAAGCCATCACCGGCATCTGCTGTTGCTGGCACCAGGGGAGTAGATCGTACTCAATTCCTCGTGAACCGAGATGG  >  NC_000913/1862974‑1863157

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